Inferences of Demography and Selection in an African Population of Drosophila melanogaster

It remains a central problem in population genetics to infer the past action of natural selection, and these inferences pose a challenge because demographic events will also substantially affect patterns of polymorphism and divergence. Thus it is imperative to explicitly model the underlying demographic history of the population whenever making inferences about natural selection. In light of the considerable interest in adaptation in African populations of Drosophila melanogaster, which are considered ancestral to the species, we generated a large polymorphism data set representing 2.1 Mb from each of 20 individuals from a Ugandan population of D. melanogaster. In contrast to previous inferences of a simple population expansion in eastern Africa, our demographic modeling of this ancestral population reveals a strong signature of a population bottleneck followed by population expansion, which has significant implications for future demographic modeling of derived populations of this species. Taking this more complex underlying demographic history into account, we also estimate a mean X-linked region-wide rate of adaptation of 6 × 10⁻¹¹/site/generation and a mean selection coefficient of beneficial mutations of 0.0009. These inferences regarding the rate and strength of selection are largely consistent with most other estimates from D. melanogaster and indicate a relatively high rate of adaptation driven by weakly beneficial mutations.     

Chromosomal divergence and evolutionary inferences in Rhodniini based on the chromosomal location of ribosomal genes

In this study, we used fluorescence in situ hybridisation to determine the chromosomal location of 45S rDNA clusters in 10 species of the tribe Rhodniini (Hemiptera: Reduviidae: Triatominae). The results showed striking inter and intraspecific variability, with the location of the rDNA clusters restricted to sex chromosomes with two patterns: either on one (X chromosome) or both sex chromosomes (X and Y chromosomes). This variation occurs within a genus that has an unchanging diploid chromosome number (2n = 22, including 20 autosomes and 2 sex chromosomes) and a similar chromosome size and genomic DNA content, reflecting a genome dynamic not revealed by these chromosome traits. The rDNA variation in closely related species and the intraspecific polymorphism in Rhodnius ecuadoriensis suggested that the chromosomal position of rDNA clusters might be a useful marker to identify recently diverged species or populations. We discuss the ancestral position of ribosomal genes in the tribe Rhodniini and the possible mechanisms involved in the variation of the rDNA clusters, including the loss of rDNA loci on the Y chromosome, transposition and ectopic pairing. The last two processes involve chromosomal exchanges between both sex chromosomes, in contrast to the widely accepted idea that the achiasmatic sex chromosomes of Heteroptera do not interchange sequences.     

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.     

Inferring the demographic history and rate of adaptive substitution in Drosophila

An important goal of population genetics is to determine the forces that have shaped the pattern of genetic variation in natural populations. We developed a maximum likelihood method that allows us to infer demographic changes and detect recent positive selection (selective sweeps) in populations of varying size from DNA polymorphism data. Applying this approach to single nucleotide polymorphism data at more than 250 noncoding loci on the X chromosome of Drosophila melanogaster from an (ancestral) African population and a (derived) European, we found that the African population expanded about 60,000 y ago and that the European population split off from the African lineage about 15,800 y ago, thereby suffering a severe population size bottleneck. We estimated that about 160 beneficial mutations (with selection coefficients s between 0.05% and 0.5%) were fixed in the euchromatic portion of the X in the African population since population size expansion, and about 60 mutations (with s around 0.5%) in the diverging European lineage.     

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.     

A scan of molecular variation leads to the narrow localization of a selective sweep affecting both Afrotropical and cosmopolitan populations of Drosophila melanogaster

Drosophila melanogaster originated in tropical Africa but has achieved a cosmopolitan distribution in association with human habitation. Cosmopolitan populations of D. melanogaster are known to have reduced genetic variation, particularly on the X chromosome. However, the relative importance of population bottlenecks and selective sweeps in explaining this reduction is uncertain. We surveyed variation at 31 microsatellites across a 330-kb section of the X chromosome located between the white and kirre genes. Two linked clusters of loci were observed with reduced variation and a skew toward rare alleles in both an Ecuador and a Zimbabwe population sample. Examining Zimbabwe DNA sequence polymorphism within one of these regions allowed us to localize a selective sweep to a 361-bp window within the 5' regulatory region of the roughest gene, with one nucleotide substitution representing the best candidate for the target of selection. Estimates of sweep age suggested that this fixation event occurred prior to the expansion of D. melanogaster from sub-Saharan Africa. For both putative sweep regions in our data set, cosmopolitan populations showed wider footprints of selection compared to those in Zimbabwe. This pattern appears consistent with the demographic amplification of preexisting sweep signals due to one or more population bottlenecks.     

Demographic Inference Reveals African and European Admixture in the North American Drosophila melanogaster Population

Drosophila melanogaster spread from sub-Saharan Africa to the rest of the world colonizing new environments. Here, we modeled the joint demography of African (Zimbabwe), European (The Netherlands), and North American (North Carolina) populations using an approximate Bayesian computation (ABC) approach. By testing different models (including scenarios with continuous migration), we found that admixture between Africa and Europe most likely generated the North American population, with an estimated proportion of African ancestry of 15%. We also revisited the demography of the ancestral population (Africa) and found—in contrast to previous work—that a bottleneck fits the history of the population of Zimbabwe better than expansion. Finally, we compared the site-frequency spectrum of the ancestral population to analytical predictions under the estimated bottleneck model.     

Population Genomics of Sub-Saharan Drosophila melanogaster: African Diversity and Non-African Admixture

Drosophila melanogaster has played a pivotal role in the development of modern population genetics. However, many basic questions regarding the demographic and adaptive history of this species remain unresolved. We report the genome sequencing of 139 wild-derived strains of D. melanogaster, representing 22 population samples from the sub-Saharan ancestral range of this species, along with one European population. Most genomes were sequenced above 25X depth from haploid embryos. Results indicated a pervasive influence of non-African admixture in many African populations, motivating the development and application of a novel admixture detection method. Admixture proportions varied among populations, with greater admixture in urban locations. Admixture levels also varied across the genome, with localized peaks and valleys suggestive of a non-neutral introgression process. Genomes from the same location differed starkly in ancestry, suggesting that isolation mechanisms may exist within African populations. After removing putatively admixed genomic segments, the greatest genetic diversity was observed in southern Africa (e.g. Zambia), while diversity in other populations was largely consistent with a geographic expansion from this potentially ancestral region. The European population showed different levels of diversity reduction on each chromosome arm, and some African populations displayed chromosome arm-specific diversity reductions. Inversions in the European sample were associated with strong elevations in diversity across chromosome arms. Genomic scans were conducted to identify loci that may represent targets of positive selection within an African population, between African populations, and between European and African populations. A disproportionate number of candidate selective sweep regions were located near genes with varied roles in gene regulation. Outliers for Europe-Africa F_ST were found to be enriched in genomic regions of locally elevated cosmopolitan admixture, possibly reflecting a role for some of these loci in driving the introgression of non-African alleles into African populations.     

Further evidence of intraspecific heterochromatin variants in Drosophila melanogaster

The analysis of inter-strain heterochromatin polymorphism in mitotic chromosomes of Drosophila melanogaster was extended to some stocks characterized by chromosomal mutations. In particular, the present investigation aims to compare, in the same cell, the quinacrine banding of two different Y chromosomes of male hybrids derived from crosses using special stocks. A direct comparison of homologous heteromorphic chromosomes in F₁ hybrids provided additional evidence of differences in the fluorescence pattern of the Y chromosome, as well as in the length of the heterochromatin segment of the X chromosome.     

Molecular population genetics of the OBP83 genomic region in Drosophila subobscura and D. guanche: contrasting the effects of natural selection and gene arrangement expansion in the patterns of nucleotide variation

Chromosomal inversion polymorphism play a major role in the evolutionary dynamics of populations and species because of their effects on the patterns of genetic variability in the genomic regions within inversions. Though there is compelling evidence for the adaptive character of chromosomal polymorphisms, the mechanisms responsible for their maintenance in natural populations is not fully understood. For this type of analysis, Drosophila subobscura is a good model species as it has a rich and extensively studied chromosomal inversion polymorphism system. Here, we examine the patterns of DNA variation in two natural populations segregating for chromosomal arrangements that differentially affect the surveyed genomic region; in particular, we analyse both nucleotide substitutions and insertion/deletion variations in the genomic region encompassing the odorant-binding protein genes Obp83a and Obp83b (Obp83 region). We show that the two main gene arrangements are genetically differentiated, but are consistent with a monophyletic origin of inversions. Nevertheless, these arrangements interchange some genetic information, likely by gene conversion. We also find that the frequency spectrum-based tests indicate that the pattern of nucleotide variation is not at equilibrium; this feature probably reflects the rapid increase in the frequency of the new gene arrangement promoted by positive selection (that is an adaptive change). Furthermore, a comparative analysis of polymorphism and divergence patterns reveals a relaxation of the functional constraints at the Obp83b gene, which might be associated with particular ecological or demographic features of the Canary island endemic species D. guanche     

Fine-Scale Heterogeneity in Crossover Rate in the garnet-scalloped Region of the Drosophila melanogaster X Chromosome

Homologous recombination affects myriad aspects of genome evolution, from standing levels of nucleotide diversity to the efficacy of natural selection. Rates of crossing over show marked variability at all scales surveyed, including species-, population-, and individual-level differences. Even within genomes, crossovers are nonrandomly distributed in a wide diversity of taxa. Although intra- and intergenomic heterogeneities in crossover distribution have been documented in Drosophila, the scale and degree of crossover rate heterogeneity remain unclear. In addition, the genetic features mediating this heterogeneity are unknown. Here we quantify fine-scale heterogeneity in crossover distribution in a 2.1-Mb region of the Drosophila melanogaster X chromosome by localizing crossover breakpoints in 2500 individuals, each containing a single crossover in this specific X chromosome region. We show 90-fold variation in rates of crossing over at a 5-kb scale, place this variation in the context of several aspects of genome evolution, and identify several genetic features associated with crossover rates. Our results shed new light on the scale and magnitude of crossover rate heterogeneity in D. melanogaster and highlight potential features mediating this heterogeneity.     

Large-Scale Selective Sweep among Segregation Distorter Chromosomes in African Populations of Drosophila melanogaster

Segregation Distorter (SD) is a selfish, coadapted gene complex on chromosome 2 of Drosophila melanogaster that strongly distorts Mendelian transmission; heterozygous SD/SD⁺ males sire almost exclusively SD-bearing progeny. Fifty years of genetic, molecular, and theory work have made SD one of the best-characterized meiotic drive systems, but surprisingly the details of its evolutionary origins and population dynamics remain unclear. Earlier analyses suggested that the SD system arose recently in the Mediterranean basin and then spread to a low, stable equilibrium frequency (1–5%) in most natural populations worldwide. In this report, we show, first, that SD chromosomes occur in populations in sub-Saharan Africa, the ancestral range of D. melanogaster, at a similarly low frequency (~2%), providing evidence for the robustness of its equilibrium frequency but raising doubts about the Mediterranean-origins hypothesis. Second, our genetic analyses reveal two kinds of SD chromosomes in Africa: inversion-free SD chromosomes with little or no transmission advantage; and an African-endemic inversion-bearing SD chromosome, SD-Mal, with a perfect transmission advantage. Third, our population genetic analyses show that SD-Mal chromosomes swept across the African continent very recently, causing linkage disequilibrium and an absence of variability over 39% of the length of the second chromosome. Thus, despite a seemingly stable equilibrium frequency, SD chromosomes continue to evolve, to compete with one another, or evade suppressors in the genome.     

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.     

Variation in Y Chromosome Segregation in Natural Populations of Drosophila melanogaster

Functional variation among Y chromosomes in natural populations of Drosophila melanogaster was assayed by a segregation study. A total of 36 Y chromosomes was extracted and ten generations of replacement backcrossing yielded stocks with Y chromosomes in two different genetic backgrounds. Eleven of the Y chromosomes were from diverse geographic origins, and the remaining 25 were from locally captured flies. Segregation of sexes in adult offspring was scored for the four possible crosses among the two backgrounds with each Y chromosome. Although the design confounds meiotic drive and effects on viability, statistical partitioning of these effects reveals significant variation among lines in Y chromosome segregation. Results are discussed in regards to models of Y-linked segregation and viability effects, which suggest that Y-linked adaptive polymorphism is unlikely.     

Genomic Patterns of Geographic Differentiation in Drosophila simulans

Geographic patterns of genetic differentiation have long been used to understand population history and to learn about the biological mechanisms of adaptation. Here we present an examination of genomic patterns of differentiation between northern and southern populations of Australian and North American Drosophila simulans, with an emphasis on characterizing signals of parallel differentiation. We report on the genomic scale of differentiation and functional enrichment of outlier SNPs. While, overall, signals of shared differentiation are modest, we find the strongest support for parallel differentiation in genomic regions that are associated with regulation. Comparisons to Drosophila melanogaster yield potential candidate genes involved in local adaptation in both species, providing insight into common selective pressures and responses. In contrast to D. melanogaster, in D. simulans we observe patterns of variation that are inconsistent with a model of temperate adaptation out of a tropical ancestral range, highlighting potential differences in demographic and colonization histories of this cosmopolitan species pair.     

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.     

Identifying new sex-linked genes through BAC sequencing in the dioecious plant Silene latifolia

Background

Silene latifolia represents one of the best-studied plant sex chromosome systems. A new approach using RNA-seq data has recently identified hundreds of new sex-linked genes in this species. However, this approach is expected to miss genes that are either not expressed or are expressed at low levels in the tissue(s) used for RNA-seq. Therefore other independent approaches are needed to discover such sex-linked genes.

Results

Here we used 10 well-characterized S. latifolia sex-linked genes and their homologs in Silene vulgaris, a species without sex chromosomes, to screen BAC libraries of both species. We isolated and sequenced 4 Mb of BAC clones of S. latifolia X and Y and S. vulgaris genomic regions, which yielded 59 new sex-linked genes (with S. vulgaris homologs for some of them). We assembled sequences that we believe represent the tip of the Xq arm. These sequences are clearly not pseudoautosomal, so we infer that the S. latifolia X has a single pseudoautosomal region (PAR) on the Xp arm. The estimated mean gene density in X BACs is 2.2 times lower than that in S. vulgaris BACs, agreeing with the genome size difference between these species. Gene density was estimated to be extremely low in the Y BAC clones. We compared our BAC-located genes with the sex-linked genes identified in previous RNA-seq studies, and found that about half of them (those with low expression in flower buds) were not identified as sex-linked in previous RNA-seq studies. We compiled a set of ~70 validated X/Y genes and X-hemizygous genes (without Y copies) from the literature, and used these genes to show that X-hemizygous genes have a higher probability of being undetected by the RNA-seq approach, compared with X/Y genes; we used this to estimate that about 30 % of our BAC-located genes must be X-hemizygous. The estimate is similar when we use BAC-located genes that have S. vulgaris homologs, which excludes genes that were gained by the X chromosome.

Conclusions

Our BAC sequencing identified 59 new sex-linked genes, and our analysis of these BAC-located genes, in combination with RNA-seq data suggests that gene losses from the S. latifolia Y chromosome could be as high as 30 %, higher than previous estimates of 10-20 %.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1698-7) contains supplementary material, which is available to authorized users.     

Genetic perspectives on “Lion Conservation Units” in Eastern and Southern Africa

Current understanding of genetic variation in lions (Panthera leo) is inadequate to guide many management decisions necessary for conservation of the species. We studied sequence variation in the mitochondrial cytochrome-b (cyt-b) gene of 75 lions and nuclear variation at 11 microsatellite loci of 480 lions across 8 range states (Cameroon, Uganda, Kenya, Zambia, Zimbabwe, South Africa, Botswana, and Namibia) and 13 Lion Conservation Units (LCUs) plus two other unassigned sites (Cameroon and Zimbabwe). A total of 11 cyt-b haplotypes were found, whose variation follows an isolation-by-distance model. In combination with previously known sequences, the haplotypes document the close relationship, derived position, and limited variability of Asian and West and Central African lions relative to other extant lions. Both phylogenetic analyses and substitution networks identify two clades in Eastern and Southern Africa—one restricted to Namibia and South Africa and the other more widespread across the region. However, these analyses are equivocal on which of these is closest to the ancestor of modern lions. Microsatellite analyses showed high levels of variation within and among populations, subdivision among most LCUs, and evidence of isolation by distance. While rates of gene flow are generally low, admixture among lions in northern Botswana, Caprivi Strip (Namibia) and Zambia is apparent from STRUCTURE analyses. Conservation management plans should incorporate information on genetic variability and gene flow in delimiting management units and in guiding translocations of lions to minimize inbreeding and to control problem animals.     

Mapping of two suppressors of OVATE (sov) loci in tomato

Tomato fruit shape varies significantly in the cultivated germplasm. To a large extent, this variation can be explained by four genes including OVATE. While most varieties with the OVATE mutation bear elongated fruits, some accessions carry round fruit, suggesting the existence of suppressors of OVATE in the germplasm. We developed three intraspecific F₂ populations with parents that carried the OVATE mutation but differed in fruit shape. We used a bulk segregant analysis approach and genotyped the extreme classes using a high-throughput genotyping platform, the SolCAP Infinium Assay. The analyses revealed segregation at two quantitative trait loci (QTLs), sov1 and sov2. These loci were confirmed by genotyping and QTL analyses of the entire population. More precise location of those loci using progeny testing confirmed that sov1 on chromosome 10 controlled obovoid and elongated shape, whereas sov2 on chromosome 11 controlled mainly elongated fruit shape. Both loci were located in intervals of <2.4 Mb on their respective chromosomes.     

Natural variation of the Y chromosome suppresses sex ratio distortion and modulates testis-specific gene expression in Drosophila simulans

X-linked sex-ratio distorters that disrupt spermatogenesis can cause a deficiency in functional Y-bearing sperm and a female-biased sex ratio. Y-linked modifiers that restore a normal sex ratio might be abundant and favored when a X-linked distorter is present. Here we investigated natural variation of Y-linked suppressors of sex-ratio in the Winters systems and the ability of these chromosomes to modulate gene expression in Drosophila simulans. Seventy-eight Y chromosomes of worldwide origin were assayed for their resistance to the X-linked sex-ratio distorter gene Dox. Y chromosome diversity caused males to sire ∼63% to ∼98% female progeny. Genome-wide gene expression analysis revealed hundreds of genes differentially expressed between isogenic males with sensitive (high sex ratio) and resistant (low sex ratio) Y chromosomes from the same population. Although the expression of about 75% of all testis-specific genes remained unchanged across Y chromosomes, a subset of post-meiotic genes was upregulated by resistant Y chromosomes. Conversely, a set of accessory gland-specific genes and mitochondrial genes were downregulated in males with resistant Y chromosomes. The D. simulans Y chromosome also modulated gene expression in XXY females in which the Y-linked protein-coding genes are not transcribed. The data suggest that the Y chromosome might exert its regulatory functions through epigenetic mechanisms that do not require the expression of protein-coding genes. The gene network that modulates sex ratio distortion by the Y chromosome is poorly understood, other than that it might include interactions with mitochondria and enriched for genes expressed in post-meiotic stages of spermatogenesis.