Dorcas J. Orengo

Detecting the Footprint of Positive Selection in a European Population of Drosophila melanogaster: Multilocus Pattern of Variation and Distance to Coding Regions

The effects on nucleotide variation of adaptations to temperate habitats and of the possible bottleneck associated with the origin of European populations of Drosophila melanogaster should be detectable in DNA sequences given the short time elapsed relative to the species population size. We surveyed nucleotide variation in 109 fragments distributed across the X chromosome in a European population of D. melanogaster to detect the footprint of positive selection. Fragments were located primarily in large noncoding regions. Multilocus tests based on Tajimas D statistic revealed a significant departure from neutral expectations in a stationary panmictic population, with an important contribution from both positive and negative D values. A positive relationship between Tajimas D values and distance to coding region was detected, with a comparative excess of significantly negative D values in the subset of fragments closer to coding regions. Also, there was a significant heterogeneity in the polymorphism to divergence ratio, with 12 fragments contributing 42% to the test statistic. Moreover, these fragments were comparatively closer to coding regions. These findings would imply positive selection events, and thus selective sweeps, during the species expansion to Europe.

Molecular Population Genetics of the Insulin/TOR Signal Transduction Pathway: A Network-Level Analysis in Drosophila melanogaster

The IT-insulin/target of rapamycin (TOR)-signal transduction pathway is a relatively well-characterized pathway that plays a central role in fundamental biological processes. Network-level analyses of DNA divergence in Drosophila and vertebrates have revealed a clear gradient in the levels of purifying selection along this pathway, with the downstream genes being the most constrained. Remarkably, this feature does not result from factors known to affect selective constraint such as gene expression, codon bias, protein length, and connectivity. The present work aims to establish whether the selective constraint gradient detected along the IT pathway at the between-species level can also be observed at a shorter time scale. With this purpose, we have surveyed DNA polymorphism in Drosophila melanogaster and divergence from D. simulans along the IT pathway. Our network-level analysis shows that DNA polymorphism exhibits the same polarity in the strength of purifying selection as previously detected at the divergence level. This equivalent feature detected both within species and between closely and distantly related species points to the action of a general mechanism, whose action is neither organism specific nor evolutionary time dependent. The detected polarity would be, therefore, intrinsic to the IT pathway architecture and function.

Relationship Between Chromosomal Polymorphism and Wing Size in a Natural Population of Drosophila Subobscura

Chromosomal polymorphism and wing size (as a measure of body size) were analysed simultaneously in two samples of Drosophila subobscura from Barcelona, Spain. The very rich chromosomal polymorphism of this species makes it difficult to detect the relationship of this polymorphism with any phenotypical character. However, a positive significant regression of wing size on the percentage of the autosome length with standard arrangement was found. Furthermore, for each polymorphic chromosome, except for the J chromosome, an association between the most frequent arrangements and wing size was observed. This trend, which was the same in the two samples, was that expected according to the latitudinal clines of both characters.

WING-SIZE HERITABILITY IN A NATURAL POPULATION OF DROSOPHILA SUBOBSCURA

Heritability of wing size was determined in a natural population of Drosophila subobscura for two consecutive year samples. In the 1988 sample, heritability in the laboratory environment was around 15%, whereas the lower bound in nature was around 0.1%. On the other hand, in the 1989 sample, heritability in the laboratory was around 90% and in nature around 15%. Differences between the two years could be caused by the more variable climate in which the males used as fathers of the 1988 sample developed. This indicates the importance of determining the environment in which parents developed in nature before discussing the heritability values obtained.

Characterization of the breakpoints of a polymorphic inversion complex detects strict and broad breakpoint reuse at the molecular level

Inversions are an integral part of structural variation within species, and they play a leading role in genome reorganization across species. Work at both the cytological and genome sequence levels has revealed heterogeneity in the distribution of inversion breakpoints, with some regions being recurrently used. Breakpoint reuse at the molecular level has mostly been assessed for fixed inversions through genome sequence comparison, and therefore rather broadly. Here, we have identified and sequenced the breakpoints of two polymorphic inversions-E1 and E2 that share a breakpoint-in the extant Est and E1 + 2 chromosomal arrangements of Drosophila subobscura. The breakpoints are two medium-sized repeated motifs that mediated the inversions by two different mechanisms: E1 via staggered breaks and subsequent repair and E2 via repeat-mediated ectopic recombination. The fine delimitation of the shared breakpoint revealed its strict reuse at the molecular level regardless of which was the intermediate arrangement. The occurrence of other rearrangements in the most proximal and distal extended breakpoint regions reveals the broad reuse of these regions. This differential degree of fragility might be related to their sharing the presence outside the inverted region of snoRNA-encoding genes.

A molecular perspective on a complex polymorphic inversion system with cytological evidence of multiply reused breakpoints

Genome sequence comparison across the Drosophila genus revealed that some fixed inversion breakpoints had been multiply reused at this long timescale. Cytological studies of Drosophila inversion polymorphism had previously shown that, also at this shorter timescale, some breakpoints had been multiply reused. The paucity of molecularly characterized polymorphic inversion breakpoints has so far precluded contrasting whether cytologically shared breakpoints of these relatively young inversions are actually reused at the molecular level. The E chromosome of Drosophila subobscura stands out because it presents several inversion complexes. This is the case of the E1+2+9+3 arrangement that originated from the ancestral Est arrangement through the sequential accumulation of four inversions (E1, E2, E9 and E3) sharing some breakpoints. We recently identified the breakpoints of inversions E1 and E2, which allowed establishing reuse at the molecular level of the cytologically shared breakpoint of these inversions. Here, we identified and sequenced the breakpoints of inversions E9 and E3, because they share breakpoints at sections 58D and 64C with those of inversions E1 and E2. This has allowed establishing that E9 and E3 originated through the staggered-break mechanism. Most importantly, sequence comparison has revealed the multiple reuse at the molecular level of the proximal breakpoint (section 58D), which would have been used at least by inversions E2, E9 and E3. In contrast, the distal breakpoint (section 64C) might have been only reused once by inversions E1 and E2, because the distal E3 breakpoint is displaced >70 kb from the other breakpoint limits.

Sequences Upstream of the Homologous cis-elements of the Adh Adult Enhancer of Drosophila Are Required for Maximal Levels of Adh Gene Transcription in Adults of Scaptodrosophila lebanonensis

The evolution of cis-regulatory elements is of particular interest for our understanding of the evolution of gene regulation. The Adh gene of Drosophilidae shows interspecific differences in tissue-specific expression and transcript levels during development. In Scaptodrosophila lebanonensis adults, the level of distal transcripts is maximal between the fourth and eighth day after eclosion and is around five times higher than that in D. melanogaster AdhS. To examine whether these quantitative differences are regulated by sequences lying upstream of the distal promoter, we performed in vitro deletion mutagenesis of the Adh gene of S. lebanonensis, followed by P-element-mediated germ-line transformation. All constructs included, as a cotransgene, a modified Adh gene of D. melanogaster (dAdh) in a fixed position and orientation that acted as a chromosomal position control. Using this approach, we have identified a fragment of 1.5 kb in the 5′ region, 830 bp upstream of the distal start site, which is required to achieve maximal levels of distal transcript in S. lebanonensis. The presence of this fragment produces a 3.5-fold higher level of distal mRNA (as determined by real time quantitative PCR) compared with the D. melanogaster dAdh cotransgene. This region contains the degenerated end of a minisatellite sequence expanding farther upstream and does not correspond to the Adh adult enhancer (AAE) of D. melanogaster. Indeed, the cis-regulatory elements of the AAE have been identified by phylogenetic footprinting within the region 830 bp upstream of the distal start site of S. lebanonensis. Furthermore, the deletions Δ-830 and Δ-2358 yield the same pattern of tissue-specific expression, indicating that all tissue-specific elements are contained within the region 830 bp upstream of the distal start site.

The origin of chromosomal inversions as a source of segmental duplications in the Sophophora subgenus of Drosophila

Chromosomal inversions can contribute to the adaptation of organisms to their environment by capturing particular advantageous allelic combinations of a set of genes included in the inverted fragment and also by advantageous functional changes due to the inversion process itself that might affect not only the expression of flanking genes but also their dose and structure. Of the two mechanisms originating inversions —ectopic recombination, and staggered double-strand breaks and subsequent repair— only the latter confers the inversion the potential to have dosage effects and/or to generate advantageous chimeric genes. In Drosophila subobscura, there is ample evidence for the adaptive character of its chromosomal polymorphism, with an important contribution of some warm-climate arrangements such as E1+2+9+12. Here, we have characterized the breakpoints of inversion E12 and established that it originated through the staggered-break mechanism like four of the five inversions of D. subobscura previously studied. This mechanism that also predominates in the D. melanogaster lineage might be prevalent in the Sophophora subgenus and contribute to the adaptive character of the polymorphic and fixed inversions of its species. Finally, we have shown that the D. subobscura inversion breakpoint regions have generally been disrupted by additional structural changes occurred at different time scales.

Uncovering the Footprint of Positive Selection on the X Chromosome of Drosophila melanogaster

A usual approach to detect the spatial footprint left by recent adaptive events has been to follow up putative candidates emerging from multilocus scans of variation by sequencing additional fragments. We have used a similar experimental and analytical approach to study variation at 15 independently evolving and randomly chosen regions of the X chromosome of Drosophila melanogaster. These incompletely sequenced regions, each extending over approximately 40 kb, were subjected to two tests of positive selection that take into account the spatial distribution of nucleotide variation. Our analysis of variation at these genomic regions in a European population of D. melanogaster has allowed us to uncover a candidate region for positive selection and to empirically evaluate the comparative performance of the two tests of selection under a bottleneck scenario. Moreover, the boundaries here estimated for both the rate of adaptive substitution (delta) and the average selection coefficient (s) would support previous estimates obtained by maximum likelihood that suggest rather strong but uncommon positive selection.

Multiple and diverse structural changes affect the breakpoint regions of polymorphic inversions across the Drosophila genus

Chromosomal polymorphism is widespread in the Drosophila genus, with extensive evidence supporting its adaptive character in diverse species. Moreover, inversions are the major contributors to the genus chromosomal evolution. The molecular characterization of a reduced number of polymorphic inversion breakpoints in Drosophila melanogaster and Drosophila subobscura supports that their inversions would have mostly originated through a mechanism that generates duplications —staggered double-strand breaks— and has thus the potential to contribute to their adaptive character. There is also evidence for inversion breakpoint reuse at different time scales. Here, we have characterized the breakpoints of two inversions of D. subobscura —O4 and O8— involved in complex arrangements that are frequent in the warm parts of the species distribution area. The duplications detected at their breakpoints are consistent with their origin through the staggered-break mechanism, which further supports it as the prevalent mechanism in D. subobscura. The comparative analysis of inversions breakpoint regions across the Drosophila genus has revealed several genes affected by multiple disruptions due not only to inversions but also to single-gene transpositions and duplications.