Molecular population genetics of Drosophila subtelomeric DNA

Telomerase-independent survival leads to a mosaic of complex subtelomere rearrangements in Chlamydomonas reinhardtii

Telomeres and subtelomeres, the genomic regions located at chromosome extremities, are essential for genome stability in eukaryotes. In the absence of the canonical maintenance mechanism provided by telomerase, telomere shortening induces genome instability. The landscape of the ensuing genome rearrangements is not accessible by short-read sequencing. Here, we leverage Oxford Nanopore Technologies long-read sequencing to survey the extensive repertoire of genome rearrangements in telomerase mutants of the model green microalga Chlamydomonas reinhardtii In telomerase-mutant strains grown for hundreds of generations, most chromosome extremities were capped by short telomere sequences that were either recruited de novo from other loci or maintained in a telomerase-independent manner. Other extremities did not end with telomeres but only with repeated subtelomeric sequences. The subtelomeric elements, including rDNA, were massively rearranged and involved in breakage-fusion-bridge cycles, translocations, recombinations, and chromosome circularization. These events were established progressively over time and displayed heterogeneity at the subpopulation level. New telomere-capped extremities composed of sequences originating from more internal genomic regions were associated with high DNA methylation, suggesting that de novo heterochromatin formation contributes to the restoration of chromosome end stability in C. reinhardtii The diversity of alternative strategies present in the same organism to maintain chromosome integrity and the variety of rearrangements found in telomerase mutants are remarkable, and illustrate genome plasticity at short timescales.

Architecture and evolution of subtelomeres in the unicellular green alga Chlamydomonas reinhardtii

In most eukaryotes, subtelomeres are dynamic genomic regions populated by multi-copy sequences of different origins, which can promote segmental duplications and chromosomal rearrangements. However, their repetitive nature has complicated the efforts to sequence them, analyse their structure and infer how they evolved. Here, we use recent genome assemblies of Chlamydomonas reinhardtii based on long-read sequencing to comprehensively describe the subtelomere architecture of the 17 chromosomes of this model unicellular green alga. We identify three main repeated elements present at subtelomeres, which we call Sultan, Subtile and Suber, alongside three chromosome extremities with ribosomal DNA as the only identified component of their subtelomeres. The most common architecture, present in 27 out of 34 subtelomeres, is a heterochromatic array of Sultan elements adjacent to the telomere, followed by a transcribed Spacer sequence, a G-rich microsatellite and transposable elements. Sequence similarity analyses suggest that Sultan elements underwent segmental duplications within each subtelomere and rearranged between subtelomeres at a much lower frequency. Analysis of other green algae reveals species-specific repeated elements that are shared across subtelomeres, with an overall organization similar to C. reinhardtii. This work uncovers the complexity and evolution of subtelomere architecture in green algae.

Architecture and evolution of subtelomeres in the unicellular green alga Chlamydomonas reinhardtii

In most eukaryotes, subtelomeres are dynamic genomic regions populated by multi-copy sequences of different origins, which can promote segmental duplications and chromosomal rearrangements. However, their repetitive nature has complicated the efforts to sequence them, analyse their structure and infer how they evolved. Here, we use recent genome assemblies of Chlamydomonas reinhardtii based on long-read sequencing to comprehensively describe the subtelomere architecture of the 17 chromosomes of this model unicellular green alga. We identify three main repeated elements present at subtelomeres, which we call Sultan, Subtile and Suber, alongside three chromosome extremities with ribosomal DNA as the only identified component of their subtelomeres. The most common architecture, present in 27 out of 34 subtelomeres, is a heterochromatic array of Sultan elements adjacent to the telomere, followed by a transcribed Spacer sequence, a G-rich microsatellite and transposable elements. Sequence similarity analyses suggest that Sultan elements underwent segmental duplications within each subtelomere and rearranged between subtelomeres at a much lower frequency. Analysis of other green algae reveals species-specific repeated elements that are shared across subtelomeres, with an overall organization similar to C. reinhardtii. This work uncovers the complexity and evolution of subtelomere architecture in green algae.

Keeping Balance Between Genetic Stability and Plasticity at the Telomere and Subtelomere of Trypanosoma brucei

Telomeres, the nucleoprotein complexes at chromosome ends, are well-known for their essential roles in genome integrity and chromosome stability. Yet, telomeres and subtelomeres are frequently less stable than chromosome internal regions. Many subtelomeric genes are important for responding to environmental cues, and subtelomeric instability can facilitate organismal adaptation to extracellular changes, which is a common theme in a number of microbial pathogens. In this review, I will focus on the delicate and important balance between stability and plasticity at telomeres and subtelomeres of a kinetoplastid parasite, Trypanosoma brucei, which causes human African trypanosomiasis and undergoes antigenic variation to evade the host immune response. I will summarize the current understanding about T. brucei telomere protein complex, the telomeric transcript, and telomeric R-loops, focusing on their roles in maintaining telomere and subtelomere stability and integrity. The similarities and differences in functions and underlying mechanisms of T. brucei telomere factors will be compared with those in human and yeast cells.

Identification and Characterization of Breakpoints and Mutations on Drosophila melanogaster Balancer Chromosomes

Balancers are rearranged chromosomes used in Drosophila melanogaster to maintain deleterious mutations in stable populations, preserve sets of linked genetic elements and construct complex experimental stocks. Here, we assess the phenotypes associated with breakpoint-induced mutations on commonly used third chromosome balancers and show remarkably few deleterious effects. We demonstrate that a breakpoint in p53 causes loss of radiation-induced apoptosis and a breakpoint in Fucosyltransferase A causes loss of fucosylation in nervous and intestinal tissue—the latter study providing new markers for intestinal cell identity and challenging previous conclusions about the regulation of fucosylation. We also describe thousands of potentially harmful mutations shared among X or third chromosome balancers, or unique to specific balancers, including an Ankyrin 2 mutation present on most TM3 balancers, and reiterate the risks of using balancers as experimental controls. We used long-read sequencing to confirm or refine the positions of two inversions with breakpoints lying in repetitive sequences and provide evidence that one of the inversions, In(2L)Cy, arose by ectopic recombination between foldback transposon insertions and the other, In(3R)C, cleanly separates subtelomeric and telomeric sequences and moves the subtelomeric sequences to an internal chromosome position. In addition, our characterization of In(3R)C shows that balancers may be polymorphic for terminal deletions. Finally, we present evidence that extremely distal mutations on balancers can add to the stability of stocks whose purpose is to maintain homologous chromosomes carrying mutations in distal genes. Overall, these studies add to our understanding of the structure, diversity and effectiveness of balancer chromosomes.

Sequence analysis of wheat subtelomeres reveals a high polymorphism among homoeologous chromosomes

Bread wheat, Triticum aestivum L., is one of the most important crops in the world. Understanding its genome organization (allohexaploid; AABBDD; 2n = 6x = 42) is essential for geneticists and plant breeders. Particularly, the knowledge of how homologous chromosomes (equivalent chromosomes from the same genome) specifically recognize each other to pair at the beginning of meiosis, the cellular process to generate gametes in sexually reproducing organisms, is fundamental for plant breeding and has a big influence on the fertility of wheat plants. Initial homologous chromosome interactions contribute to specific recognition and pairing between homologues at the onset of meiosis. Understanding the molecular basis of these critical processes can help to develop genetic tools in a breeding context to promote interspecific chromosome associations in hybrids or interspecific genetic crosses to facilitate the transfer of desirable agronomic traits from related species into a crop like wheat. The terminal regions of chromosomes, which include telomeres and subtelomeres, participate in chromosome recognition and pairing. We present a detailed molecular analysis of subtelomeres of wheat chromosome arms 1AS, 4AS, 7AS, 7BS and 7DS. Results showed a high polymorphism in the subtelomeric region among homoeologues (equivalent chromosomes from related genomes) for all the features analyzed, including genes, transposable elements, repeats, GC content, predicted CpG islands, recombination hotspots and targeted sequence motifs for relevant DNA-binding proteins. These polymorphisms might be the molecular basis for the specificity of homologous recognition and pairing in initial chromosome interactions at the beginning of meiosis in wheat.

The Telomere Paradox: Stable Genome Preservation with Rapidly Evolving Proteins

Telomeres ensure chromosome length homeostasis and protection from catastrophic end-to-end chromosome fusions. All eukaryotes require this essential, strictly conserved telomere-dependent genome preservation. However, recent evolutionary analyses of mammals, plants, and flies report pervasive rapid evolution of telomere proteins. The causes of this paradoxical observation - that unconserved machinery underlies an essential, conserved function - remain enigmatic. Indeed, these fast-evolving telomere proteins bind, extend, and protect telomeric DNA, which itself evolves slowly in most systems. We hypothesize that the universally fast-evolving subtelomere - the telomere-adjacent, repetitive sequence - is a primary driver of the 'telomere paradox'. Under this model, radical sequence changes in the subtelomere perturb subtelomere-dependent, telomere functions. Compromised telomere function then spurs adaptation of telomere proteins to maintain telomere length homeostasis and protection. We propose an experimental framework that leverages both protein divergence and subtelomeric sequence divergence to test the hypothesis that subtelomere sequence evolution shapes recurrent innovation of telomere machinery.

Altered Localization of Hybrid Incompatibility Proteins in Drosophila

Understanding the molecular basis of hybrid incompatibilities is a fundamental pursuit in evolutionary genetics. In crosses between Drosophila melanogaster females and Drosophila simulans males, an interaction between at least three genes is necessary for hybrid male lethality: Hmr ^mel, Lhr ^sim, and gfzf ^sim. Although HMR and LHR physically bind each other and function together in a single complex, the connection between gfzf and either of these proteins remains mysterious. Here, we show that GFZF localizes to many regions of the genome in both D. melanogaster and D. simulans, including at telomeric retrotransposon repeats. We find that GFZF localization at telomeres is significantly different between these two species, reflecting the rapid evolution of telomeric retrotransposon copy number composition between the two species. Next, we show that GFZF and HMR normally do not colocalize in D. melanogaster. In interspecies hybrids, however, HMR shows extensive mis-localization to GFZF sites, thus uncovering a new molecular interaction between these hybrid incompatibility factors. We find that spreading of HMR to GFZF sites requires gfzf ^sim but not Lhr ^sim, suggesting distinct roles for these factors in the hybrid incompatibility. Finally, we find that overexpression of HMR and LHR within species is sufficient to mis-localize HMR to GFZF binding sites, indicating that HMR has a natural low affinity for GFZF sites. Together, these studies provide the first insights into the different properties of gfzf between D. melanogaster and D. simulans, and uncover a molecular interaction between gfzf and Hmr in the form of altered protein localization.

Telomere and Subtelomere R-loops and Antigenic Variation in Trypanosomes

Trypanosoma brucei is a kinetoplastid parasite that causes African trypanosomiasis, which is fatal if left untreated. T. brucei regularly switches its major surface antigen, VSG, to evade the host immune responses. VSGs are exclusively expressed from subtelomeric expression sites (ESs) where VSG genes are flanked by upstream 70 bp repeats and downstream telomeric repeats. The telomere downstream of the active VSG is transcribed into a long-noncoding RNA (TERRA), which forms RNA:DNA hybrids (R-loops) with the telomeric DNA. At an elevated level, telomere R-loops cause more telomeric and subtelomeric double-strand breaks (DSBs) and increase VSG switching rate. In addition, stabilized R-loops are observed at the 70 bp repeats and immediately downstream of ES-linked VSGs in RNase H defective cells, which also have an increased amount of subtelomeric DSBs and more frequent VSG switching. Although subtelomere plasticity is expected to be beneficial to antigenic variation, severe defects in subtelomere integrity and stability increase cell lethality. Therefore, regulation of the telomere and 70 bp repeat R-loop levels is important for the balance between antigenic variation and cell fitness in T. brucei. In addition, the high level of the active ES transcription favors accumulation of R-loops at the telomere and 70 bp repeats, providing an intrinsic mechanism for local DSB formation, which is a strong inducer of VSG switching.

Functional diversification accompanies gene family expansion of MED2 homologs in Candida albicans

Gene duplication facilitates functional diversification and provides greater phenotypic flexibility to an organism. Expanded gene families arise through repeated gene duplication but the extent of functional divergence that accompanies each paralogous gene is generally unexplored because of the difficulty in isolating the effects of single family members. The telomere-associated (TLO) gene family is a remarkable example of gene family expansion, with 14 members in the more pathogenic Candida albicans relative to two TLO genes in the closely-related species C. dubliniensis. TLO genes encode interchangeable Med2 subunits of the major transcriptional regulatory complex Mediator. To identify biological functions associated with each C. albicans TLO, expression of individual family members was regulated using a Tet-ON system and the strains were assessed across a range of phenotypes involved in growth and virulence traits. All TLOs affected multiple phenotypes and a single phenotype was often affected by multiple TLOs, including simple phenotypes such as cell aggregation and complex phenotypes such as virulence in a Galleria mellonella model of infection. No phenotype was regulated by all TLOs, suggesting neofunctionalization or subfunctionalization of ancestral properties among different family members. Importantly, regulation of three phenotypes could be mapped to individual polymorphic sites among the TLO genes, including an indel correlated with two phenotypes, growth in sucrose and macrophage killing. Different selective pressures have operated on the TLO sequence, with the 5’ conserved Med2 domain experiencing purifying selection and the gene/clade-specific 3’ end undergoing extensive positive selection that may contribute to the impact of individual TLOs on phenotypic variability. Therefore, expansion of the TLO gene family has conferred unique regulatory properties to each paralog such that it influences a range of phenotypes. We posit that the genetic diversity associated with this expansion contributed to C. albicans success as a commensal and opportunistic pathogen.

Gene duplication is a rapid mechanism to generate additional sequences for natural selection to act upon and confer greater organismal fitness. If additional copies of the gene are beneficial, this process may be repeated to produce an expanded gene family containing many copies of related sequences. Following duplication, individual gene family members may retain functions of the ancestral gene or acquire new functions through mutation. How functional diversification accompanies expansion into large gene families remains largely unexplored due to the difficulty in assessing individual genes in the presence of the remaining family members. Here, we addressed this question using an inducible promoter to regulate expression of individual genes of the TLO gene family in the commensal yeast and opportunistic pathogen Candida albicans, which encode components of a major transcriptional regulator. Induced expression of individual TLOs affected a wide range of phenotypes such that significant functional overlap occurred among TLO genes and most phenotypes were affected by more than one TLO. Induced expression of individual TLOs did not produce massive phenotypic effects in most cases, suggesting that functional overlap among TLO genes may buffer new mutations that arise. Specific sequence variants among the TLO genes correlated with certain phenotypes and these sequence variants did not necessarily correlate with sequence similarity across the entire gene. Therefore, individual TLO family members evolved specific functional roles following duplication that likely reflect a combination of inherited function and new mutation.

The Molecular and Genetic Characterization of Second Chromosome Balancers in Drosophila melanogaster

Balancer chromosomes are multiply inverted and rearranged chromosomes used in Drosophila melanogaster for many tasks, such as maintaining mutant alleles in stock and complex stock construction. Balancers were created before molecular characterization of their breakpoints was possible, so the precise locations of many of these breakpoints are unknown. Here, we report or confirm the positions of the 14 euchromatic breakpoints on the 2^nd chromosome balancers SM1, SM5, CyO, and SM6a This total includes three breakpoints involved in a complex rearrangement on SM5 that is associated with the duplication of two genomic regions. Unbiased sequencing of several balancers allowed us to identify stocks with incorrectly identified balancers as well as single and double crossover events that had occurred between 2^nd chromosome balancers and their homologs. The confirmed crossover events that we recovered were at least 2 Mb from the closest inversion breakpoint, consistent with observations from other balancer chromosomes. Balancer chromosomes differ from one another both by large tracts of sequence diversity generated by recombination and by small differences, such as single nucleotide polymorphisms (SNPs). Therefore, we also report loss-of-function mutations carried by these chromosomes and unique SNP and InDel polymorphisms present on only single balancers. These findings provide valuable information about the structure of commonly used 2^nd chromosome balancers and extend recent work examining the structure of X and 3^rd chromosome balancers. Finally, these observations provide new insights into how the sequences of individual balancers have diverged over time.

Short and long-term evolutionary dynamics of subtelomeric piRNA clusters in Drosophila

Two Telomeric Associated Sequences, TAS-R and TAS-L, form the principal subtelomeric repeat families identified in Drosophila melanogaster. They are PIWI-interacting RNA (piRNA) clusters involved in repression of Transposable Elements. In this study, we revisited TAS structural and functional dynamics in D. melanogaster and in related species. In silico analysis revealed that TAS-R family members are composed of previously uncharacterized domains. This analysis also showed that TAS-L repeats are composed of arrays of a region we have named "TAS-L like" (TLL) identified specifically in one TAS-R family member, X-TAS. TLL were also present in other species of the melanogaster subgroup. Therefore, it is possible that TLL represents an ancestral subtelomeric piRNA core-cluster. Furthermore, all D. melanogaster genomes tested possessed at least one TAS-R locus, whereas TAS-L can be absent. A screen of 110 D. melanogaster lines showed that X-TAS is always present in flies living in the wild, but often absent in long-term laboratory stocks and that natural populations frequently lost their X-TAS within 2 years upon lab conditioning. Therefore, the unexpected structural and temporal dynamics of subtelomeric piRNA clusters demonstrated here suggests that genome organization is subjected to distinct selective pressures in the wild and upon domestication in the laboratory.

Targeted identification of TE insertions in a Drosophila genome through hemi-specific PCR

BACKGROUND

Transposable elements (TEs) are major components of eukaryotic genomes and drivers of genome evolution, producing intraspecific polymorphism and interspecific differences through mobilization and non-homologous recombination. TE insertion sites are often highly variable within species, creating a need for targeted genome re-sequencing (TGS) methods to identify TE insertion sites.

METHODS

We present a hemi-specific PCR approach for TGS of P-elements in Drosophila genomes on the Illumina platform. We also present a computational framework for identifying new insertions from TGS reads. Finally, we describe a new method for estimating the frequency of TE insertions from WGS data, which is based precise insertion sites provided by TGS annotations.

RESULTS

By comparing our results to TE annotations based on whole genome re-sequencing (WGS) data for the same Drosophilamelanogaster strain, we demonstrate that TGS is powerful for identifying true insertions, even in repeat-rich heterochromatic regions. We also demonstrate that TGS offers enhanced annotation of precise insertion sites, which facilitates estimation of TE insertion frequency.

CONCLUSIONS

TGS by hemi-specific PCR is a powerful approach for identifying TE insertions of particular TE families in species with a high-quality reference genome, at greatly reduced cost as compared to WGS. It may therefore be ideal for population genomic studies of particular TE families. Additionally, TGS and WGS can be used as complementary approaches, with TGS annotations identifying more annotated insertions with greater precision for a target TE family, and WGS data allowing for estimates of TE insertion frequencies, and a broader picture of the location of non-target TEs across the genome.

Contrasting evolutionary genome dynamics between domesticated and wild yeasts

Structural rearrangements have long been recognized as an important source of genetic variation, with implications in phenotypic diversity and disease, yet their detailed evolutionary dynamics remain elusive. Here we use long-read sequencing to generate end-to-end genome assemblies for 12 strains representing major subpopulations of the partially domesticated yeast Saccharomyces cerevisiae and its wild relative Saccharomyces paradoxus. These population-level high-quality genomes with comprehensive annotation enable precise definition of chromosomal boundaries between cores and subtelomeres and a high-resolution view of evolutionary genome dynamics. In chromosomal cores, S. paradoxus shows faster accumulation of balanced rearrangements (inversions, reciprocal translocations and transpositions), whereas S. cerevisiae accumulates unbalanced rearrangements (novel insertions, deletions and duplications) more rapidly. In subtelomeres, both species show extensive interchromosomal reshuffling, with a higher tempo in S. cerevisiae. Such striking contrasts between wild and domesticated yeasts are likely to reflect the influence of human activities on structural genome evolution.

A Pooled Sequencing Approach Identifies a Candidate Meiotic Driver in Drosophila

Meiotic drive occurs when a selfish element increases its transmission frequency above the Mendelian ratio by hijacking the asymmetric divisions of female meiosis. Meiotic drive causes genomic conflict and potentially has a major impact on genome evolution, but only a few drive loci of large effect have been described. New methods to reliably detect meiotic drive are therefore needed, particularly for discovering moderate-strength drivers that are likely to be more prevalent in natural populations than strong drivers. Here, we report an efficient method that uses sequencing of large pools of backcross (BC1) progeny to test for deviations from Mendelian segregation genome-wide with single-nucleotide polymorphisms (SNPs) that distinguish the parental strains. We show that meiotic drive can be detected by a characteristic pattern of decay in distortion of SNP frequencies, caused by recombination unlinking the driver from distal loci. We further show that control crosses allow allele-frequency distortion caused by meiotic drive to be distinguished from distortion resulting from developmental effects. We used this approach to test whether chromosomes with extreme telomere-length differences segregate at Mendelian ratios, as telomeric regions are a potential hotspot for meiotic drive due to their roles in meiotic segregation and multiple observations of high rates of telomere sequence evolution. Using four different pairings of long and short telomere strains, we find no evidence that extreme telomere-length variation causes meiotic drive in Drosophila However, we identify one candidate meiotic driver in a centromere-linked region that shows an ∼8% increase in transmission frequency, corresponding to a ∼54:46 segregation ratio. Our results show that candidate meiotic drivers of moderate strength can be readily detected and localized in pools of BC1 progeny.

Recurrent Innovation at Genes Required for Telomere Integrity in Drosophila

Telomeres are nucleoprotein complexes at the ends of linear chromosomes. These specialized structures ensure genome integrity and faithful chromosome inheritance. Recurrent addition of repetitive, telomere-specific DNA elements to chromosome ends combats end-attrition, while specialized telomere-associated proteins protect naked, double-stranded chromosome ends from promiscuous repair into end-to-end fusions. Although telomere length homeostasis and end-protection are ubiquitous across eukaryotes, there is sporadic but building evidence that the molecular machinery supporting these essential processes evolves rapidly. Nevertheless, no global analysis of the evolutionary forces that shape these fast-evolving proteins has been performed on any eukaryote. The abundant population and comparative genomic resources of Drosophila melanogaster and its close relatives offer us a unique opportunity to fill this gap. Here we leverage population genetics, molecular evolution, and phylogenomics to define the scope and evolutionary mechanisms driving fast evolution of genes required for telomere integrity. We uncover evidence of pervasive positive selection across multiple evolutionary timescales. We also document prolific expansion, turnover, and expression evolution in gene families founded by telomeric proteins. Motivated by the mutant phenotypes and molecular roles of these fast-evolving genes, we put forward four alternative, but not mutually exclusive, models of intra-genomic conflict that may play out at very termini of eukaryotic chromosomes. Our findings set the stage for investigating both the genetic causes and functional consequences of telomere protein evolution in Drosophila and beyond.

Real-Time Evolution of a Subtelomeric Gene Family in Candida albicans

Subtelomeric regions of the genome are notable for high rates of sequence evolution and rapid gene turnover. Evidence of subtelomeric evolution has relied heavily on comparisons of historical evolutionary patterns to infer trends and frequencies of these events. Here, we describe evolution of the subtelomeric TLO gene family in Candida albicans during laboratory passaging for over 4000 generations. C. albicans is a commensal and opportunistic pathogen of humans and the TLO gene family encodes a subunit of the Mediator complex that regulates transcription and affects a range of virulence factors. We identified 16 distinct subtelomeric recombination events that altered the TLO repertoire. Ectopic recombination between subtelomeres on different chromosome ends occurred approximately once per 5000 generations and was often followed by loss of heterozygosity, resulting in the complete loss of one TLO gene sequence with expansion of another. In one case, recombination within TLO genes produced a novel TLO gene sequence. TLO copy number changes were biased, with some TLOs preferentially being copied to novel chromosome arms and other TLO genes being frequently lost. The majority of these nonreciprocal recombination events occurred either within the 3′ end of the TLO coding sequence or within a conserved 50-bp sequence element centromere-proximal to TLO coding sequence. Thus, subtelomeric recombination is a rapid mechanism of generating genotypic diversity through alterations in the number and sequence of related gene family members.

Sperm should evolve to make female meiosis fair

Genomic conflicts arise when an allele gains an evolutionary advantage at a cost to organismal fitness. Oögenesis is inherently susceptible to such conflicts because alleles compete for inclusion into the egg. Alleles that distort meiosis in their favor (i.e., meiotic drivers) often decrease organismal fitness, and therefore indirectly favor the evolution of mechanisms to suppress meiotic drive. In this light, many facets of oögenesis and gametogenesis have been interpreted as mechanisms of protection against genomic outlaws. That females of many animal species do not complete meiosis until after fertilization, appears to run counter to this interpretation, because this delay provides an opportunity for sperm-acting alleles to meddle with the outcome of female meiosis and help like alleles drive in heterozygous females. Contrary to this perceived danger, the population genetic theory presented herein suggests that, in fact, sperm nearly always evolve to increase the fairness of female meiosis in the face of genomic conflicts. These results are consistent with the apparent sperm dependence of the best characterized female meiotic driversin animals. Rather than providing an opportunity for sperm collaboration in female meiotic drive, the "fertilization requirement" indirectly protects females from meiotic drivers by providing sperm an opportunity to suppress drive.

The relations between recombination rate and patterns of molecular variation and evolution in Drosophila

Genetic recombination affects levels of variability and the efficacy of selection because natural selection acting at one site affects evolutionary processes at linked sites. The variation in local recombination rates across the Drosophila genome provides excellent material for testing hypotheses concerning the evolutionary consequences of recombination. The current state of knowledge from studies of Drosophila genomics and population genetics is reviewed here. Selection at linked sites has influenced the relations between recombination rates and patterns of molecular variation and evolution, such that higher rates of recombination are associated with both higher levels of variability and a greater efficacy of selection. It seems likely that background selection against deleterious mutations is a major factor contributing to these patterns in genome regions in which crossing over is rare or absent, whereas selective sweeps of positively selected mutations probably play an important role in regions with crossing over.

The relation between recombination rate and patterns of molecular evolution and variation in Drosophila melanogaster

Genetic recombination associated with sexual reproduction increases the efficiency of natural selection by reducing the strength of Hill–Robertson interference. Such interference can be caused either by selective sweeps of positively selected alleles or by background selection (BGS) against deleterious mutations. Its consequences can be studied by comparing patterns of molecular evolution and variation in genomic regions with different rates of crossing over. We carried out a comprehensive study of the benefits of recombination in Drosophila melanogaster, both by contrasting five independent genomic regions that lack crossing over with the rest of the genome and by comparing regions with different rates of crossing over, using data on DNA sequence polymorphisms from an African population that is geographically close to the putatively ancestral population for the species, and on sequence divergence from a related species. We observed reductions in sequence diversity in noncrossover (NC) regions that are inconsistent with the effects of hard selective sweeps in the absence of recombination. Overall, the observed patterns suggest that the recombination rate experienced by a gene is positively related to an increase in the efficiency of both positive and purifying selection. The results are consistent with a BGS model with interference among selected sites in NC regions, and joint effects of BGS, selective sweeps, and a past population expansion on variability in regions of the genome that experience crossing over. In such crossover regions, the X chromosome exhibits a higher rate of adaptive protein sequence evolution than the autosomes, implying a Faster-X effect.

Identification of two Penelope-like elements with different structures and chromosome localization in kuruma shrimp genome

Penelope, originally found as a key element responsible for the hybrid dysgenesis in Drosophila virilis, has been widely conserved throughout eukaryotic genomes. In other organisms, they are often referred to as Penelope-like elements or PLEs. In this study, we found two types of PLEs, designated MjPLE01 and MjPLE02, from kuruma shrimp, Marsupenaeus japonicus. There was no observed nucleotide similarity between MjPLE01 and 02, and both elements differed from each other in terms of their structure; MjPLE02 has a distinctive endonuclease (EN) domain at the C-terminus while MjPLE01 do not. A phylogenetic tree that includes publicly available PLEs and TERTs showed that MjPLE01 and 02 were closely related to Coprina elements, which have been reported as an EN-deficient PLE, and to Penelope-Poseidon group, which possess an EN domain, respectively. Genomic Southern blot analysis using MjPLE01 as a probe showed several multiple bands that differ among individual shrimps. On the other hand, two major identical bands were observed when MjPLE02 was used. Colony hybridization showed co-localization of MjPLE01 and GGTTA repeats, suggesting that MjPLE01 might be prevalent in subtelomeric regions of kuruma shrimp genome. These results suggest that the kuruma shrimp genome has at least two types of PLEs with different domain compositions, phylogenetic positions, and probably chromosomeal localization. Such distinctive types of PLEs in an organism have never been described and hence could be a potential source to understand how multiple PLE types evolved.

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.

Recombination is a process by which chromosomes exchange genetic material during meiosis. It is important in evolution because it provides offspring with new combinations of genes, and so estimating the rate of recombination is of fundamental importance in various population genomic inference problems. In this paper, we develop a new statistical method to enable robust estimation of fine-scale recombination maps of Drosophila, a genus of common fruit flies, in which the background recombination rate is high and natural selection has been prevalent. We apply our method to produce fine-scale recombination maps for a North American population and an African population of D. melanogaster. For both populations, we find extensive fine-scale variation in recombination rate throughout the genome. We provide a quantitative characterization of the similarities and differences between the recombination maps of the two populations; our study reveals high correlation at broad scales and low correlation at fine scales, as has been documented among human populations. We also examine the correlation between various genomic features. Furthermore, using a conservative approach, we find a handful of putative recombination “hotspot” regions with solid statistical support for a local elevation of at least 10 times the background recombination rate.

Telomere-targeted retrotransposons in the rice blast fungus Magnaporthe oryzae: agents of telomere instability

The fungus Magnaporthe oryzae is a serious pathogen of rice and other grasses. Telomeric restriction fragments in Magnaporthe isolates that infect perennial ryegrass (prg) are hotspots for genomic rearrangement and undergo frequent, spontaneous alterations during fungal culture. The telomeres of rice-infecting isolates are very stable by comparison. Sequencing of chromosome ends from a number of prg-infecting isolates revealed two related non-LTR retrotransposons (M. oryzae Telomeric Retrotransposons or MoTeRs) inserted in the telomere repeats. This contrasts with rice pathogen telomeres that are uninterrupted by other sequences. Genetic evidence indicates that the MoTeR elements are responsible for the observed instability. MoTeRs represent a new family of telomere-targeted transposons whose members are found exclusively in fungi.

Evaluating the evidence for transmission distortion in human pedigrees

Children of a heterozygous parent are expected to carry either allele with equal probability. Exceptions can occur, however, due to meiotic drive, competition among gametes, or viability selection, which we collectively term “transmission distortion” (TD). Although there are several well-characterized examples of these phenomena, their existence in humans remains unknown. We therefore performed a genome-wide scan for TD by applying the transmission disequilibrium test (TDT) genome-wide to three large sets of human pedigrees of European descent: the Framingham Heart Study (FHS), a founder population of European origin (HUTT), and a subset of the Autism Genetic Resource Exchange (AGRE). Genotyping error is an important confounder in this type of analysis. In FHS and HUTT, despite extensive quality control, we did not find sufficient evidence to exclude genotyping error in the strongest signals. In AGRE, however, many signals extended across multiple SNPs, a pattern highly unlikely to arise from genotyping error. We identified several candidate regions in this data set, notably a locus in 10q26.13 displaying a genome-wide significant TDT in combined female and male transmissions and a signature of recent positive selection, as well as a paternal TD signal in 6p21.1, the same region in which a significant TD signal was previously observed in 30 European males. Neither region replicated in FHS, however, and the paternal signal was not visible in sperm competition assays or as allelic imbalance in sperm. In maternal transmissions, we detected no strong signals near centromeres or telomeres, the regions predicted to be most susceptible to female-specific meiotic drive, but we found a significant enrichment of top signals among genes involved in cell junctions. These results illustrate both the potential benefits and the challenges of using the TDT to study transmission distortion and provide candidates for investigation in future studies.

Sequencing of pooled DNA samples (Pool-Seq) uncovers complex dynamics of transposable element insertions in Drosophila melanogaster

Transposable elements (TEs) are mobile genetic elements that parasitize genomes by semi-autonomously increasing their own copy number within the host genome. While TEs are important for genome evolution, appropriate methods for performing unbiased genome-wide surveys of TE variation in natural populations have been lacking. Here, we describe a novel and cost-effective approach for estimating population frequencies of TE insertions using paired-end Illumina reads from a pooled population sample. Importantly, the method treats insertions present in and absent from the reference genome identically, allowing unbiased TE population frequency estimates. We apply this method to data from a natural Drosophila melanogaster population from Portugal. Consistent with previous reports, we show that low recombining genomic regions harbor more TE insertions and maintain insertions at higher frequencies than do high recombining regions. We conservatively estimate that there are almost twice as many “novel” TE insertion sites as sites known from the reference sequence in our population sample (6,824 novel versus 3,639 reference sites, with on average a 31-fold coverage per insertion site). Different families of transposable elements show large differences in their insertion densities and population frequencies. Our analyses suggest that the history of TE activity significantly contributes to this pattern, with recently active families segregating at lower frequencies than those active in the more distant past. Finally, using our high-resolution TE abundance measurements, we identified 13 candidate positively selected TE insertions based on their high population frequencies and on low Tajima's D values in their neighborhoods.

Transposable elements (TE's) are parasitic genetic elements that spread by replicating themselves within a host genome. Most organisms are burdened with transposable elements; in fact, up to 80% of some genomes can consist of TE–derived DNA. Here, we use new sequencing technology to examine variation in genomic TE composition within a population at a finer scale and in a more unbiased fashion than has been possible before. We study a Portuguese population of D. melanogaster and find a large number of TE insertions, most of which occur in few individuals. Our analysis confirms that TE insertions are subject to purifying selection that counteracts their spread, and it suggests that the genome records waves of past TE invasions, with recently active elements occurring at low population frequency. We also find indications that TE insertions may sometimes have beneficial effects.

Scrambling eggs: meiotic drive and the evolution of female recombination rates

Theories to explain the prevalence of sex and recombination have long been a central theme of evolutionary biology. Yet despite decades of attention dedicated to the evolution of sex and recombination, the widespread pattern of sex differences in the recombination rate is not well understood and has received relatively little theoretical attention. Here, we argue that female meiotic drivers--alleles that increase in frequency by exploiting the asymmetric cell division of oogenesis--present a potent selective pressure favoring the modification of the female recombination rate. Because recombination plays a central role in shaping patterns of variation within and among dyads, modifiers of the female recombination rate can function as potent suppressors or enhancers of female meiotic drive. We show that when female recombination modifiers are unlinked to female drivers, recombination modifiers that suppress harmful female drive can spread. By contrast, a recombination modifier tightly linked to a driver can increase in frequency by enhancing female drive. Our results predict that rapidly evolving female recombination rates, particularly around centromeres, should be a common outcome of meiotic drive. We discuss how selection to modify the efficacy of meiotic drive may contribute to commonly observed patterns of sex differences in recombination.

Inferring the history of interchromosomal gene transposition in Drosophila using n-dimensional parsimony

Gene transposition puts a new gene copy in a novel genomic environment. Moreover, genes moving between the autosomes and the X chromosome experience change in several evolutionary parameters. Previous studies of gene transposition have not utilized the phylogenetic framework that becomes possible with the availability of whole genomes from multiple species. Here we used parsimonious reconstruction on the genomic distribution of gene families to analyze interchromosomal gene transposition in Drosophila. We identified 782 genes that have moved chromosomes within the phylogeny of 10 Drosophila species, including 87 gene families with multiple independent movements on different branches of the phylogeny. Using this large catalog of transposed genes, we detected accelerated sequence evolution in duplicated genes that transposed when compared to the parental copy at the original locus. We also observed a more refined picture of the biased movement of genes from the X chromosome to the autosomes. The bias of X-to-autosome movement was significantly stronger for RNA-based movements than for DNA-based movements, and among DNA-based movements there was an excess of genes moving onto the X chromosome as well. Genes involved in female-specific functions moved onto the X chromosome while genes with male-specific functions moved off the X. There was a significant overrepresentation of proteins involving chromosomal function among transposed genes, suggesting that genetic conflict between sexes and among chromosomes may be a driving force behind gene transposition in Drosophila.

Genetic variation and DNA replication timing, or why is there late replicating DNA?

Mutation rates vary significantly within the genome and across species. Recent studies revealed a long suspected replication-timing effect on mutation rate, but the mechanisms that regulate the increase in mutation rate as the genome is replicated remain unclear. Evidence is emerging, however, that DNA repair systems, in general, are less efficient in late replicating heterochromatic regions compared to early replicating euchromatic regions of the genome. At the same time, mutation rates in both vertebrates and invertebrates have been shown to vary with generation time (GT). GT is correlated with genome size, which suggests a possible nucleotypic effect on species-specific mutation rates. These and other observations all converge on a role for DNA replication checkpoints in modulating generation times and mutation rates during the DNA synthetic phase (S phase) of the cell cycle. The following will examine the potential role of the intra-S checkpoint in regulating cell cycle times (GT) and mutation rates in eukaryotes. This article was published online on August 5, 2011. An error was subsequently identified. This notice is included in the online and print versions to indicate that both have been corrected October 4, 2011.

Differential maintenance of DNA sequences in telomeric and centromeric heterochromatin

Repeated DNA in heterochromatin presents enormous difficulties for whole-genome sequencing; hence, sequence organization in a significant portion of the genomes of multicellular organisms is relatively unknown. Two sequenced BACs now allow us to compare telomeric retrotransposon arrays from Drosophila melanogaster telomeres with an array of telomeric retrotransposons that transposed into the centromeric region of the Y chromosome >13 MYA, providing a unique opportunity to compare the structural evolution of this retrotransposon in two contexts. We find that these retrotransposon arrays, both heterochromatic, are maintained quite differently, resulting in sequence organizations that apparently reflect different roles in the two chromosomal environments. The telomere array has grown only by transposition of new elements to the chromosome end; the centromeric array instead has grown by repeated amplifications of segments of the original telomere array. Many elements in the telomere have been variably 5'-truncated apparently by gradual erosion and irregular deletions of the chromosome end; however, a significant fraction (4 and possibly 5 or 6 of 15 elements examined) remain complete and capable of further retrotransposition. In contrast, each element in the centromere region has lost ≥ 40% of its sequence by internal, rather than terminal, deletions, and no element retains a significant part of the original coding region. Thus the centromeric array has been restructured to resemble the highly repetitive satellite sequences typical of centromeres in multicellular organisms, whereas, over a similar or longer time period, the telomere array has maintained its ability to provide retrotransposons competent to extend telomere ends.

Three selections are better than one: clinal variation of thermal QTL from independent selection experiments in Drosophila

We report the results of two independent selection experiments that have exposed distinct populations of Drosophila melanogaster to different forms of thermal selection. A recombinant population derived from Arvin California and Zimbabwe isofemale lines was exposed to laboratory natural selection at two temperatures (T(AZ): 18°C and 28°C). Microsatellite mapping identified quantitative trait loci (QTL) on the X-chromosome between the replicate "Hot" and "Cold" populations. In a separate experiment, disruptive selection was imposed on an outbred California population for the "knockdown" temperature (T(KD)) in a thermal column. Microsatellite mapping of the "High" and "Low" populations also uncovered primarily X-linked QTL. Notably, a marker in the shaggy locus at band 3A was significantly differentiated in both experiments. Finer scale mapping of the 3A region has narrowed the QTL to the shaggy gene region, which contains several candidate genes that function in circadian rhythms. The same allele that was increased in frequency in the High T(KD) populations is significantly clinal in North America and is more common at the warm end of the cline (Florida vs. Maine; however, the cline was not apparent in Australia). Together, these studies show that independent selection experiments can uncover the same target of selection and that evolution in the laboratory can recapitulate putatively adaptive clinal variation in nature.

Rapid expansion and functional divergence of subtelomeric gene families in yeasts

BACKGROUND

Subtelomeres, regions proximal to telomeres, exhibit characteristics unique to eukaryotic genomes. Genes residing in these loci are subject to epigenetic regulation and elevated rates of both meiotic and mitotic recombination. However, most genome sequences do not contain assembled subtelomeric sequences, and, as a result, subtelomeres are often overlooked in comparative genomics.

RESULTS

We studied the evolution and functional divergence of subtelomeric gene families in the yeast lineage. Our computational results show that subtelomeric families are evolving and expanding much faster than families that do not contain subtelomeric genes. Focusing on three related subtelomeric MAL gene families involved in disaccharide metabolism that show typical patterns of rapid expansion and evolution, we show experimentally how frequent duplication events followed by functional divergence yield novel alleles that allow the metabolism of different carbohydrates.

CONCLUSIONS

Taken together, our computational and experimental analyses show that the extraordinary instability of eukaryotic subtelomeres supports rapid adaptation to novel niches by promoting gene recombination and duplication followed by functional divergence of the alleles.

A 'higher order' of telomere regulation: telomere heterochromatin and telomeric RNAs

Protection of chromosome ends from DNA repair and degradation activities is mediated by specialized protein complexes bound to telomere repeats. Recently, it has become apparent that epigenetic regulation of the telomric chromatin template critically impacts on telomere function and telomere-length homeostasis from yeast to man. Across all species, telomeric repeats as well as the adjacent subtelomeric regions carry features of repressive chromatin. Disruption of this silent chromatin environment results in loss of telomere-length control and increased telomere recombination. In turn, progressive telomere loss reduces chromatin compaction at telomeric and subtelomeric domains. The recent discoveries of telomere chromatin regulation during early mammalian development, as well as during nuclear reprogramming, further highlights a central role of telomere chromatin changes in ontogenesis. In addition, telomeres were recently shown to generate long, non-coding RNAs that remain associated to telomeric chromatin and will provide new insights into the regulation of telomere length and telomere chromatin. In this review, we will discuss the epigenetic regulation of telomeres across species, with special emphasis on mammalian telomeres. We will also discuss the links between epigenetic alterations at mammalian telomeres and telomere-associated diseases.

Unusual composition of a yeast chromosome arm is associated with its delayed replication

The 11.3-Mb genome of the yeast Lachancea (Saccharomyces) kluyveri displays an intriguing compositional heterogeneity: a region of approximately 1 Mb, covering almost the whole left arm of chromosome C (C-left), has an average GC content of 52.9%, which is significantly higher than the 40.4% global GC content of the rest of the genome. This region contains the MAT locus, which remains normal in composition. The excess of GC base pairs affects both coding and noncoding sequences, and thus is not due to selective pressure acting on protein sequences. It leads to a strong codon usage bias and alters the amino acid composition of the 457 proteins encoded on C-left that do not show obvious bias for functional categories, or the presence of paralogs or orthologs of essential genes of Saccharomyces cerevisiae. They share significant synteny conservation with other species of the Saccharomycetaceae, and phylogenetic analysis indicates that C-left originates from a Lachancea species. In contrast, there is a complete absence of transposable elements in C-left, whereas 18 elements per megabase are distributed across the rest of the genome. Comparative hybridization of synchronized cells using high-density genome arrays reveals that C-left is replicated later during S phase than the rest of the genome. Two possible primary causes of this major compositional heterogeneity are discussed: an ancient hybridization of two related species with very distinct GC composition, or an intrinsic mechanism, possibly associated with the loss of the silent cassettes from C-left that progressively increased the GC content and generated the delayed replication of this chromosomal arm.

The evolutionary significance of copy number variation in the human genome

Copy number variation provides the raw material for gene family expansion and diversification, which is an important evolutionary force. Moreover, copy number variants (CNVs) can influence gene transcriptional and translational levels and have been associated with complex disease susceptibility. Therefore, natural selection may have affected at least some of the greater than one thousand CNVs thus far discovered among the genomes of phenotypically normal humans. While identifying and understanding particular instances of natural selection may shed light on important aspects of human evolutionary history, our ability to analyze CNVs in traditional population genetic frameworks has been limited. However, progress has been made by adapting some of these frameworks for use with copy number data. Moving forward, these efforts will be aided by non-human organism studies of the population genetics of copy number variation, and by more direct comparisons of within-species copy number variation and between-species copy number fixation.

Molecular population genetics and evolution of Drosophila meiosis genes

While many functional elements of the meiotic process are well characterized in model organisms, the genetic basis of most of the natural phenotypic variation observed in meiotic pathways has not been determined. To begin to address this issue, we characterized patterns of polymorphism and divergence in the protein-coding regions of 33 genes across 31 lines of Drosophila melanogaster and 6 lines of Drosophila simulans. We sequenced genes known to be involved in chromosome segregation, recombination, DNA repair, and related heterochromatin binding. As expected, we found several of the genes to be highly conserved, consistent with purifying selection. However, a subset of genes showed patterns of polymorphism and divergence typical of other types of natural selection. Moreover, several intriguing differences between the two Drosophila lineages were evident: along the D. simulans lineage we consistently found evidence of adaptive protein evolution, whereas along the D. melanogaster lineage several loci exhibited patterns consistent with the maintenance of protein variation.

Silent but not static: accelerated base-pair substitution in silenced chromatin of budding yeasts

Subtelomeric DNA in budding yeasts, like metazoan heterochromatin, is gene poor, repetitive, transiently silenced, and highly dynamic. The rapid evolution of subtelomeric regions is commonly thought to arise from transposon activity and increased recombination between repetitive elements. However, we found evidence of an additional factor in this diversification. We observed a surprising level of nucleotide divergence in transcriptionally silenced regions in inter-species comparisons of Saccharomyces yeasts. Likewise, intra-species analysis of polymorphisms also revealed increased SNP frequencies in both intergenic and synonymous coding positions of silenced DNA. This analysis suggested that silenced DNA in Saccharomyces cerevisiae and closely related species had increased single base-pair substitution that was likely due to the effects of the silencing machinery on DNA replication or repair.

Many plants, fungi, pathogens, and animals have chromosome regions that are silenced. Special proteins change the chromosome structure in these domains, turning genes off or lowering their expression levels. We found an increased frequency of DNA mutations in these silenced regions of closely related yeasts. This increase is likely due to silencing proteins interfering with DNA repair or replication. Accurate replication of genetic information with minimal mutations is usually critical for the survival and fitness of an organism; however, there are examples where a high mutation rate is beneficial. The silenced regions of chromosomes are often associated with virus-like transposable elements, and with genes that are important in responding to environmental changes. Hence, it is possible that elevated DNA mutations in silenced regions contribute to genome defense against transposable elements or increased genetic diversity to cope with variation in surrounding conditions.

The subtelomere of Oryza sativa chromosome 3 short arm as a hot bed of new gene origination in rice

Despite general observations of non-random genomic distribution of new genes, it is unclear whether or not new genes preferentially occur in certain genomic regions driven by related molecular mechanisms. Using 1.5 Mb of genomic sequences from short arms of chromosome 3 of Oryza glaberrima and O. punctata, we conducted a comparative genomic analysis with the reference O. sativa ssp. japonica genome. We identified a 60-kb segment located in the middle of the subtelomeric region of chromosome 3, which is unique to the species O. sativa. The region contained gene duplicates that occurred in Asian cultivated rice species that diverged from the ancestor of Asian and African cultivated rice one million years ago (MYA). For the 12 genes and one complete retrotransposon identified in this segment in O. sativa ssp. japonica, we searched for their parental genes. The high similarity between duplicated paralogs further supports the recent origination of these genes. We found that this segment was recently generated through multiple independent gene recombination and transposon insertion events. Among the 12 genes, we found that five had chimeric gene structures derived from multiple parental genes. Nine out of the 12 new genes seem to be functional, as suggested by Ka/Ks analysis and the presence of cDNA and/or MPSS data. Furthermore, for the eight transcribed genes, at least two genes could be classified as defense or stress response-related genes. Given these findings, and the fact that subtelomeres are associated with high rates of recombination and transcription, it is likely that subtelomeres may facilitate gene recombination and transposon insertions and serve as hot spots for new gene origination in rice genomes.

Recurrent deletion and gene presence/absence polymorphism: telomere dynamics dominate evolution at the tip of 3L in Drosophila melanogaster and D. simulans

Although Drosophila melanogaster has been the subject of intensive analysis of polymorphism and divergence, little is known about the distribution of variation at the most distal regions of chromosomes arms. Here we report a survey of genetic variation on the tip of 3L in D. melanogaster and D. simulans. Levels of single nucleotide polymorphism in the most distal euchromatic sequence are approximately one order of magnitude less than that typically observed in genomic regions of normal crossing over, consistent with what might be expected under models of linked selection in regions of low crossing over. However, despite this reduced level of nucleotide variation, we found abundant deletion polymorphism. These deletions create at least three gene presence/absence polymorphisms within D. melanogaster: the putative G-protein coupled receptor mthl-8 (which is the most distal known or predicted gene on 3L) and the unannotated mRNAs AY060886 and BT006009. Strikingly, D. simulans is also segregating deletions that cause mthl8 presence/absence polymorphism. Breakpoint sequencing and tests of correlations with segregating SNPs in D. melanogaster suggest that each deletion is unique. Cloned breakpoint sequences revealed the presence of Het-A elements just distal to unique, canonical euchromatic sequences. This pattern suggests a model in which repeated telomeric deficiencies cause deletions of euchromatic sequence followed by subsequent "healing" by retrotranposition of Het-A elements. These data reveal the dominance of telomeric dynamics on the evolution of closely linked sequences in Drosophila.