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.     

Natural Selection Shapes Variation in Genome-wide Recombination Rate in Drosophila pseudoobscura

1. Download : Download high-res image (210KB)
2. Download : Download full-size image

     

Specificity of RNA Binding by CPEB: Requirement for RNA Recognition Motifs and a Novel Zinc Finger

CPEB is an RNA binding protein that interacts with the maturation-type cytoplasmic polyadenylation element (CPE) (consensus UUUUUAU) to promote polyadenylation and translational activation of maternal mRNAs in Xenopus laevis. CPEB, which is conserved from mammals to invertebrates, is composed of three regions: an amino-terminal portion with no obvious functional motif, two RNA recognition motifs (RRMs), and a cysteine-histidine region that is reminiscent of a zinc finger. In this study, we investigated the physical properties of CPEB required for RNA binding. CPEB can interact with RNA as a monomer, and phosphorylation, which modifies the protein during oocyte maturation, has little effect on RNA binding. Deletion mutations of CPEB have been overexpressed in Escherichia coli and used in a series of RNA gel shift experiments. Although a full-length and a truncated CPEB that lacks 139 amino-terminal amino acids bind CPE-containing RNA avidly, proteins that have had either RRM deleted bind RNA much less efficiently. CPEB that has had the cysteine-histidine region deleted has no detectable capacity to bind RNA. Single alanine substitutions of specific cysteine or histidine residues within this region also abolish RNA binding, pointing to the importance of this highly conserved domain of the protein. Chelation of metal ions by 1,10-phenanthroline inhibits the ability of CPEB to bind RNA; however, RNA binding is restored if the reaction is supplemented with zinc. CPEB also binds other metals such as cobalt and cadmium, but these destroy RNA binding. These data indicate that the RRMs and a zinc finger region of CPEB are essential for RNA binding.     

Genetic and Evolutionary Correlates of Fine-Scale Recombination Rate Variation in Drosophila persimilis

Recombination is fundamental to meiosis in many species and generates variation on which natural selection can act, yet fine-scale linkage maps are cumbersome to construct. We generated a fine-scale map of recombination rates across two major chromosomes in Drosophila persimilis using 181 SNP markers spanning two of five major chromosome arms. Using this map, we report significant fine-scale heterogeneity of local recombination rates. However, we also observed “recombinational neighborhoods,” where adjacent intervals had similar recombination rates after excluding regions near the centromere and telomere. We further found significant positive associations of fine-scale recombination rate with repetitive element abundance and a 13-bp sequence motif known to associate with human recombination rates. We noted strong crossover interference extending 5–7 Mb from the initial crossover event. Further, we observed that fine-scale recombination rates in D. persimilis are strongly correlated with those obtained from a comparable study of its sister species, D. pseudoobscura. We documented a significant relationship between recombination rates and intron nucleotide sequence diversity within species, but no relationship between recombination rate and intron divergence between species. These results are consistent with selection models (hitchhiking and background selection) rather than mutagenic recombination models for explaining the relationship of recombination with nucleotide diversity within species. Finally, we found significant correlations between recombination rate and GC content, supporting both GC-biased gene conversion (BGC) models and selection-driven codon bias models. Overall, this genome-enabled map of fine-scale recombination rates allowed us to confirm findings of broader-scale studies and identify multiple novel features that merit further investigation.     

Selective Sensitization of Zinc Finger Protein Oxidation by Reactive Oxygen Species through Arsenic Binding

Cysteine oxidation induced by reactive oxygen species (ROS) on redox-sensitive targets such as zinc finger proteins plays a critical role in redox signaling and subsequent biological outcomes. We found that arsenic exposure led to oxidation of certain zinc finger proteins based on arsenic interaction with zinc finger motifs. Analysis of zinc finger proteins isolated from arsenic-exposed cells and zinc finger peptides by mass spectrometry demonstrated preferential oxidation of C3H1 and C4 zinc finger configurations. C2H2 zinc finger proteins that do not bind arsenic were not oxidized by arsenic-generated ROS in the cellular environment. The findings suggest that selectivity in arsenic binding to zinc fingers with three or more cysteines defines the target proteins for oxidation by ROS. This represents a novel mechanism of selective protein oxidation and demonstrates how an environmental factor may sensitize certain target proteins for oxidation, thus altering the oxidation profile and redox regulation.     

Water-Induced Finger Wrinkles Do Not Affect Touch Acuity or Dexterity in Handling Wet Objects

Human non-hairy (glabrous) skin of the fingers, palms and soles wrinkles after prolonged exposure to water. Wrinkling is a sympathetic nervous system-dependent process but little is known about the physiology and potential functions of water-induced skin wrinkling. Here we investigated the idea that wrinkling might improve handling of wet objects by measuring the performance of a large cohort of human subjects (n = 40) in a manual dexterity task. We also tested the idea that skin wrinkling has an impact on tactile acuity or vibrotactile sensation using two independent sensory tasks. We found that skin wrinkling did not improve dexterity in handling wet objects nor did it affect any aspect of touch sensitivity measured. Thus water-induced wrinkling appears to have no significant impact on tactile driven performance or dexterity in handling wet or dry objects.     

Disparity and Evolutionary Rate Do Not Explain Diversity Patterns in Muroid Rodents (Rodentia: Muroidea)

A positive correlation between diversity and disparity/evolutionary rate is predicted by multiple evolutionary theories. However, recent empirical studies in various taxa do not always find such an association. Similarly, we find no correlation between these two levels of variation, based on cranial morphometric data and molecular phylogenetic data from 317 muroid rodent species and dipodoid outgroups, analyzed using three-dimensional geometric morphometrics. This disassociation was found using both phylogenetic and non-phylogenetic approaches, indicating that an increase in clade richness is not necessarily followed by an increase in morphological divergence and vice versa. Furthermore, the distribution of muroid families in morphospace is highly overlapping suggesting greater variation within than between clades. Taken together with the observation that families with the most distinctive cranial morphologies (nesomyids, dipodids, and spalacids) are the least diverse, indicates that evolution of new cranial morphologies may not play an important role in the diversification of muroid rodents.     

Structure of Genetic Variation within and between Populations of Mycophagous Drosophila

Patterns of genetic variation within and between populations of five species of mycophagous Drosophila were examined by gel electrophoresis of several polymorphic loci. Populations of the five species could not be shown to be subdivided into sympatric host-adapted races. Statistically significant, but small, between-host differences in gene frequencies were observed at three of 15 loci. Mean gene frequencies at all loci were similar in New York and Tennessee, and, with one exception, relatively little genetic differentiation was observed among study sites within those two regions. Gene frequencies generally were stable over several years of collecting as well. The unpredictable nature of the fungal hosts may preclude the site fidelity and continuity of diversifying selection necessary for adaptive divergence of populations.     

Natural Selection and Recombination Rate Variation Shape Nucleotide Polymorphism Across the Genomes of Three Related Populus Species

A central aim of evolutionary genomics is to identify the relative roles that various evolutionary forces have played in generating and shaping genetic variation within and among species. Here we use whole-genome resequencing data to characterize and compare genome-wide patterns of nucleotide polymorphism, site frequency spectrum, and population-scaled recombination rates in three species of Populus: Populus tremula, P. tremuloides, and P. trichocarpa. We find that P. tremuloides has the highest level of genome-wide variation, skewed allele frequencies, and population-scaled recombination rates, whereas P. trichocarpa harbors the lowest. Our findings highlight multiple lines of evidence suggesting that natural selection, due to both purifying and positive selection, has widely shaped patterns of nucleotide polymorphism at linked neutral sites in all three species. Differences in effective population sizes and rates of recombination largely explain the disparate magnitudes and signatures of linked selection that we observe among species. The present work provides the first phylogenetic comparative study on a genome-wide scale in forest trees. This information will also improve our ability to understand how various evolutionary forces have interacted to influence genome evolution among related species.     

Variation of body size within and between wild populations of Drosophila buzzatii

Four populations of the cactophilous species D. buzzatii have been compared with respect to the phenotypic variation of thorax and wing length of wild versus laboratory reared flies. Three of the strains were intercrossed to provide parent, F₁ and F₂ comparisons as a test of co-adaptation. The genetic contribution to phenotypic variation of laboratory reared flies was estimated from the correlation between sibs derived from random pair mating and reared individually in separate cultures. The average natural temperature during development was estimated from the relations between the wing/thorax ratio and temperature in laboratory tests.The variance of thorax and wing length of wild flies was several times greater than that of laboratory reared flies and the increase was attributed primarily to variation in larval food supply although temperature fluctuation is also important. There was no evidence of heterosis or F₂ break-down in the crosses. For two of the populations the heritability of thorax length was high, 60–70%, and substantially lower for the third. The average temperature estimated from the wing/thorax/temperature relationship differed between sites. The reduction of body size below the potential maximum averaged 30% for two and 20% for the other population, with a wide spread about these values. The evidence is discussed in relation to assessing the nature of ecological variation by comparing the variation of morphological traits in wild and laboratory reared flies.     

Phylogenetic Constraints Do Not Explain the Rarity of Nitrogen-Fixing Trees in Late-Successional Temperate Forests

Background

Symbiotic nitrogen (N)-fixing trees are rare in late-successional temperate forests, even though these forests are often N limited. Two hypotheses could explain this paradox. The ‘phylogenetic constraints hypothesis’ states that no late-successional tree taxa in temperate forests belong to clades that are predisposed to N fixation. Conversely, the ‘selective constraints hypothesis’ states that such taxa are present, but N-fixing symbioses would lower their fitness. Here we test the phylogenetic constraints hypothesis.

Methodology/Principal Findings

Using U.S. forest inventory data, we derived successional indices related to shade tolerance and stand age for N-fixing trees, non-fixing trees in the ‘potentially N-fixing clade’ (smallest angiosperm clade that includes all N fixers), and non-fixing trees outside this clade. We then used phylogenetically independent contrasts (PICs) to test for associations between these successional indices and N fixation. Four results stand out from our analysis of U.S. trees. First, N fixers are less shade-tolerant than non-fixers both inside and outside of the potentially N-fixing clade. Second, N fixers tend to occur in younger stands in a given geographical region than non-fixers both inside and outside of the potentially N-fixing clade. Third, the potentially N-fixing clade contains numerous late-successional non-fixers. Fourth, although the N fixation trait is evolutionarily conserved, the successional traits are relatively labile.

Conclusions/Significance

These results suggest that selective constraints, not phylogenetic constraints, explain the rarity of late-successional N-fixing trees in temperate forests. Because N-fixing trees could overcome N limitation to net primary production if they were abundant, this study helps to understand the maintenance of N limitation in temperate forests, and therefore the capacity of this biome to sequester carbon.     

Germinal Excisions of the Maize Transposon Activator Do Not Stimulate Meiotic Recombination or Homology-Dependent Repair at the Bz Locus

Double-strand breaks have been implicated both in the initiation of meiotic recombination in yeast and as intermediates in the transposition process of nonreplicative transposons. Some transposons of this class, notably P of Drosophila and Tc1 of Caenorhabditis elegans, promote a form of homology-dependent premeiotic gene conversion upon excision. In this work, we have looked for evidence of an interaction between Ac transposition and meiotic recombination at the bz locus in maize. We find that the frequency of meiotic recombination between homologues is not enhanced by the presence of Ac in one of the bz heteroalleles and, conversely, that the presence of a homologous sequence in either trans (homologous chromosome) or cis (tandem duplication) does not promote conversion of the Ac insertion site. However, a tandem duplication of the bz locus may be destabilized by the insertion of Ac. We discuss possible reasons for the lack of interaction between Ac excision and homologous meiotic recombination in maize.     

DNA Sequence Variation at the Period Locus within and among Species of the Drosophila Melanogaster Complex

A 1.9-kilobase region of the period locus was sequenced in six individuals of Drosophila melanogaster and from six individuals of each of three sibling species: Drosophila simulans, Drosophila sechellia and Drosophila mauritiana. Extensive genealogical analysis of 174 polymorphic sites reveals a complex history. It appears that D. simulans, as a large population still segregating very old lineages, gave rise to the island species D. mauritiana and D. sechellia. Rather than considering these speciation events as having produced ``sister' taxa, it seems more appropriate to consider D. simulans a parent species to D. sechellia and D. mauritiana. The order, in time, of these two phylogenetic events remains unclear. D. mauritiana supports a large number of polymorphisms, many of which are shared with D. simulans, and so appears to have begun and persisted as a large population. In contrast, D. sechellia has very little variation and seems to have experienced a severe population bottleneck. Alternatively, the low variation in D. sechellia could be due to recent directional selection and genetic hitchhiking at or near the per locus.     

Calories Do Not Explain Extension of Life Span by Dietary Restriction in Drosophila

Dietary restriction (DR) extends life span in diverse organisms, including mammals, and common mechanisms may be at work. DR is often known as calorie restriction, because it has been suggested that reduction of calories, rather than of particular nutrients in the diet, mediates extension of life span in rodents. We here demonstrate that extension of life span by DR in Drosophila is not attributable to the reduction in calorie intake. Reduction of either dietary yeast or sugar can reduce mortality and extend life span, but by an amount that is unrelated to the calorie content of the food, and with yeast having a much greater effect per calorie than does sugar. Calorie intake is therefore not the key factor in the reduction of mortality rate by DR in this species.     

Length Variation of Cag/Caa Trinucleotide Repeats in Natural Populations of Drosophila Melanogaster and Its Relation to the Recombination Rate

Eleven genes distributed along the Drosophila melanogaster chromosome 2 and showing exonic tandem repeats of glutamine codons (CAG or CAA) were surveyed for length variation in a sample of four European and African populations. Only one gene was monomorphic. Eight genes were polymorphic in all populations, with a total number of alleles varying between five and 12 for 120 chromosomes. The average heterozygozity per locus and population was 0.41. Selective neutrality in length variation could not be rejected under the assumptions of the infinite allele model. Significant population subdivision was found though no geographical pattern emerged, all populations being equally different. Significant linkage disequilibrium was found in four out of seven cases where the genetic distance between loci was <1 cM and was negligible when the distance was larger. There is evidence that these associations were established after the populations separated. An unexpected result was that variation at each locus was independent of the coefficient of exchange, although the latter ranged from zero to the relatively high value of 6.7%. This would indicate that background selection and selective hitchhiking, which are thought to affect levels of nucleotide substitution polymorphism, have no effect on trinucleotide repeat variation.     

Variation within and between Closely Related Species Uncovers High Intra-Specific Variability in Dispersal

Mounting evidence shows that contrasting selection pressures generate variability in dispersal patterns among individuals or populations of the same species, with potential impacts on both species dynamics and evolution. However, this variability is hardly considered in empirical works, where a single dispersal function is considered to adequately reflect the species-specific dispersal ability, suggesting thereby that within-species variation is negligible as regard to inter-specific differences in dispersal abilities. We propose here an original method to make the comparison of intra- and inter-specific variability in dispersal, by decomposing the diversity of that trait along a phylogeny of closely related species. We used as test group European butterflies that are classic study organisms in spatial ecology. We apply the analysis separately to eight metrics that reflect the dispersal propensity, the dispersal ability or the dispersal efficiency of populations and species. At the inter-specific level, only the dispersal ability showed the signature of a phylogenetic signal while neither the dispersal propensity nor the dispersal efficiency did. At the within-species level, the partitioning of dispersal diversity showed that dispersal was variable or highly variable among populations: intra-specific variability represented from 11% to 133% of inter-specific variability in dispersal metrics. This finding shows that dispersal variation is far from negligible in the wild. Understanding the processes behind this high within-species variation should allow us to properly account for dispersal in demographic models. Accordingly, to encompass the within species variability in life histories the use of more than one value per trait per species should be encouraged in the construction of databases aiming at being sources for modelling purposes.