Transposable elements (TEs) allow rewiring of regulatory networks, and the recent amplification of the ISX element dispersed 77 functional but suboptimal binding sites for the dosage compensation complex to a newly formed X chromosome in Drosophila. Here we identify two linked refining mutations within ISX that interact epistatically to increase binding affinity to the dosage compensation complex. Selection has increased the frequency of this derived haplotype in the population, which is fixed at 30% of ISX insertions and polymorphic among another 41%. Sharing of this haplotype indicates that high levels of gene conversion among ISX elements allow them to ‘crowd-source’ refining mutations, and a refining mutation that occurs at any single ISX element can spread in two dimensions: horizontally across insertion sites by non-allelic gene conversion, and vertically through the population by natural selection. These results describe a novel route by which fully functional regulatory elements can arise rapidly from TEs and implicate non-allelic gene conversion as having an important role in accelerating the evolutionary fine-tuning of regulatory networks. ELife digest. Mutations change genes and provide the raw material for evolution.

Genes are sections of DNA that contain the instructions for making proteins or other molecules, and so determine the physical characteristics of each organism. Genetic mutations that increase an organism's number of offspring and chances of survival are more likely to be passed on to future generations. Changes to when or where a gene is switched on (so-called regulatory mutations) can also provide fitness benefits and can therefore be selected for during evolution.

Transposable elements are sequences of DNA that are also called ‘jumping genes’ because they can make copies of themselves and these copies of the transposable element can move to other locations in the genome. Some transposable elements contain sequences that switch on nearby genes. If different copies of a transposable element that contains such a regulatory sequence insert themselves in more than one place, it can result in a network of genes that can all be controlled in the same way.

The regulatory sequences contained within transposable elements are not always optimal, but they can be fine-tuned through evolution. A fruit fly called Drosophila miranda has a transposable element called ISX that has, over time, placed up to 77 regulatory sequences around one of this species' sex chromosomes.

Just as in humans, female flies are XX and males are XY; but having only one copy of the X chromosome means that male flies need to increase the expression of certain genes to produce a full-dose of the molecules made by the genes. This process is called dosage compensation and in 2013 the 77 ISX regulatory sequences on the fruit fly's X chromosome were shown to help recruit the molecular machinery that carries out dosage compensation to nearby genes, albeit inefficiently. Now Ellison and Bachtrog—who also conducted the 2013 study—report how these transposable elements have been fine-tuned to make them more effective for dosage compensation. Ellison and Bachtrog uncovered two mutations that make the ISX transposable element better at recruiting the dosage compensation molecular machinery.

ISX spread around different locations along the fly's X chromosome before these mutations arose; this means that initially none of the 77 insertions carried the two mutations, but now 30% of the 77 elements have the mutations in all flies, and 41% have them in only some flies. The same mutations have spread between the different ISX elements because transposable elements with the mutations have been used to directly convert other ISX elements without them. These mutations have also become more common in the fruit fly population by being passed on to offspring and increasing their survival. These two routes have accelerated the fine-tuning of these transposable elements for use in gene regulation. This implies that regulatory sequences derived from transposable elements evolve in a way that is fundamentally different from those that arise by other means, as the direct conversion between these insertions allows fine-tuning mutations to spread more rapidly. A substantial portion of animal genomes is composed of repetitive sequences, including gene duplicates, satellite DNA, and transposable elements. Gene conversion is a major force shaping the evolution of repetitive regions, and interlocus or non-allelic gene conversion between sequence duplicates has been studied extensively for its role in concerted evolution (; ).

Non-allelic gene conversion also affects selection operating in gene families. Compared to single-copy genes, a family of gene duplicates presents a larger mutational target, and a mutation arising in any gene copy can be spread among copies by non-allelic gene conversion, thereby increasing the efficiency of both positive and purifying selection (). Non-allelic gene conversion homogenizes the arrays of ribosomal DNA gene copies present in the genomes of most organisms (), has generated allelic diversity within the human leukocyte antigen gene family (), and has allowed palindromic genes on the human Y chromosome to escape degeneration (). Transposable elements give rise to families of duplicate sequences. A propensity for some TEs to carry regulatory motifs and to insert adjacent to coding sequence gives them the potential for being potent modulators of gene regulatory networks (; ).

The regulatory elements provided by these TEs, however, may be suboptimal in function, and subject to subsequent fine-tuning (). Unlike regulatory elements where short binding motifs (10 basepairs on average for transcription factors; ) evolve de novo via point mutation or microsatellite expansion, binding sites that evolve from TEs are initially almost identical in sequence and are nested within a larger repeat unit (hundreds or thousands of basepairs in size), and may thus be subject to non-allelic gene conversion. Re-wiring of the dosage compensation network in Drosophila miranda was driven by TE-mediated amplification of a functional but suboptimal binding motif ().

Here we show that non-allelic gene conversion is catalyzing the rapid fine-tuning of these suboptimal motifs by allowing sequence variants that optimize binding affinity to spread among elements. Dosage compensation in Drosophila is mediated by a male-specific ribonucleoprotein complex (the male-specific lethal or MSL complex) that binds to a GA-rich sequence motif (the MSL recognition motif) at a number of chromatin entry sites on the X chromosome (; ).

We previously studied the acquisition of novel chromatin entry sites on newly formed X chromosomes in D. Miranda, a species where two independent sex chromosome/autosome fusions resulted in a karyotype composed of three X chromosome arms, each of a different age (; ). XL is homologous to the X chromosome of Drosophila melanogaster and has been a sex chromosome for at least 60 million years (); chromosome XR formed roughly 15 million years ago when an autosome (Muller element D) fused to XL (), and the neo-X/neo-Y chromosome pair originated around 1.5 million years ago when the Y fused to another autosome (Muller element C) (). Dosage compensation evolved on both XR and the neo-X shortly after their emergence, through acquisition of novel chromatin entry sites and co-option of the MSL regulatory network (; ). Interestingly, we discovered that the acquisition of dosage compensation on both XR and the neo-X chromosome was in part mediated by the independent domestication of helitron transposable elements that contained MSL recognition motifs, which we have termed ISXR and ISX, respectively (). ISX is highly enriched on the neo-X chromosome of D.

Phone, Apple TV are trademarks of Apple Inc. The Carter 3 Rapidshare Download Accelerator. I bought rapidshare premium account, how can i. STEP 3: High lite some rapidshare links. How to use Download Accelerator Plus with Rapidshare. Behavioral Retargeting Advertising.

Miranda and is derived from the abundant ISY element. Compared to ISY, ISX contains a 10 basepair deletion that creates a MSL recognition motif, thereby allowing it to act as a chromatin entry site (). Our previous study showed that while amplification of ISX about 1 million years ago provided dozens of functional chromatin entry sites on the neo-X chromosome of D. Miranda, the motif dispersed by ISX is distinct from the canonical motif that is enriched within chromatin entry sites on XL and XR, and shows significantly lower affinity to the MSL complex compared to motifs on XL and XR (). For these reasons, we postulated that the ISX binding motif is suboptimal, and predicted that refining mutations should accumulate within each MSL recognition motif until the neo-X chromosome becomes fully dosage compensated ().

To identify potential refining mutations that optimize MSL-binding at chromatin entry sites derived from the ISX element, we characterized sequence variation within the MSL recognition motifs and flanking sequence regions for all 77 insertions of the ISX element on the neo-X chromosome in the sequenced reference strain MSH22 (). Because we have previously demonstrated that ISX contains a functional MSL recognition motif but the closely related ISY element does not (), we sought to identify sequence variants that were present in multiple ISX elements but rare or absent in ISY elements from the same chromosome. ( A) The MSL recognition motif (MRE) plus 20 basepairs of flanking sequence were extracted from all 77 ISX transposable elements located on the neo-X chromosome in the MSH22 reference genome assembly.

The multiple sequence alignment of these 77 sequence regions (arranged from top-to-bottom in the order in which they are found on the chromosome) shows that there is sequence variation among elements both within and adjacent to the 21 basepair MRE motif. Each variant has been classified as ancestral or derived based on its frequency in the ISX progenitor element, ISY. The derived allele frequency for each variant in this region is shown for ISX as well as 139 ISY elements from the neo-X chromosome (see for ISY alignment). Red arrows point to the derived TT haplotype that is common among ISX elements but rare in ISY. ( B) Barplot showing the frequencies of all haplotypes at the GA/TT sites, for ISY and ISX elements separately. Two haplotypes are present within ISX elements (GA and TT) and the two alleles within each haplotype are in perfect linkage disequilibrium.

In contrast, the majority of ISY elements harbor the GA haplotype, but these two alleles are not in perfect linkage disequilibrium among ISY elements. Suzuki Serial Numbers Year. Rather, five additional allelic combinations are present at low frequencies in this location among ISY, but not ISX elements.

Using these criteria, we identified a sequence haplotype adjacent to the MSL recognition motif that is common among MSH22 ISX insertions and rare among ISY elements: 57% of ISX elements carry this haplotype vs 0.7% of neo-X ISY insertions, an asymmetry significantly different from that expected by chance (Fisher's Exact Test; p. ( A) MSL3 ChIP-seq data from D. Miranda strain MSH22 shows that the ISX insertions carrying the TT haplotype recruit significantly higher levels of MSL complex compared to those with the GA haplotype (Wilcoxon test p = 0.01).

( B) Engineered ISX elements that differ only with respect to the TT haplotype bind different levels of MSL complex. There is an epistatic interaction between the two ‘T’ alleles such that separately, they decrease MSL complex binding relative to the ancestral allele, but together in the TT haplotype, they increase MSL complex binding (Wilcoxon Test p = 0.028 for both comparisons [GT vs TT and TA vs TT]). The rectangles and error bars show the average and standard deviation of values from four biological replicates for each condition. We previously demonstrated that insertion of an ISX element in the D. Melanogaster genome results in recruitment of the MSL complex to an ectopic autosomal location (). We used this same system to dissect the relationship between the TT alleles and MSL complex binding affinity. Starting with a cloned ISX element (), we used site-directed mutagenesis to create variants of this element that differ only with respect to the TT haplotype.

Each of the four possible haplotypes (GA, GT, TA, and TT) was engineered and inserted onto D. Melanogaster chromosome 2L at cytosite 38F1 using recombinase mediated cassette exchange (RMCE) ().

We then measured the effect of each of the derived variants by quantifying allele-specific binding levels of the MSL complex in F1 hybrids between the ancestral haplotype (GA) and each of the derived haplotypes (GT, TA, and TT). Interestingly, each T allele, when assayed separately, has a negative effect on MSL binding levels compared to the ancestral G or A allele (). However, when combined, the TT haplotype results in significantly increased levels of MSL complex binding, relative to the ancestral GA haplotype (Wilcoxon Test p = 0.0289; ). These results suggest that there is sign epistasis between the two alleles and that the high frequency TT haplotype represents a refining/fine-tuning adaptation, since recruitment of MSL complex to the adjacent MSL recognition motif is increased. Non-allelic gene conversion is spreading the TT haplotype among ISX insertions.

It is unlikely that the TT haplotype arose multiple times by parallel mutation, and there are two possibilities that could explain its prevalence among MSH22 ISX insertions. First, this double mutation may have occurred early during the process of ISX amplification, thus giving rise to two lineages of ISX: one that carries the ancestral GA haplotype, and the other with the TT haplotype. The TT-harboring elements in MSH22 would then all be descendants from the latter ISX lineage. Alternatively, this mutation may have occurred only after the GA-containing ISX element was fixed in the population at all 77 neo-X insertion sites, at which point it was spread among independent ISX elements via non-allelic gene conversion. We can distinguish between these possibilities by examining patterns of sequence polymorphism for each ISX insertion across multiple strains of D.

A canonical signature of non-allelic gene conversion is the presence of shared polymorphisms across sequence duplicates (; ). If gene conversion is spreading the TT haplotype among ISX insertions, we expect it to be polymorphic among individuals at several ISX insertion sites, whereas we do not expect the TT haplotype to be polymorphic at individual ISX insertions under the alternative scenario. To genotype multiple wild-derived individuals at each of the MSH22 ISX insertions, we used paired-end Illumina genomic resequencing data from 23 inbred lines of D.

Miranda, including MSH22. We aligned all reads to the MSH22 reference genome and identified mate-pairs where one mate was anchored in unique sequence flanking an ISX insertion. We then assembled these reads to generate a contig spanning the 5′ flank of the ISX element insertion, which contains the MSL recognition motif, for each inbred line.

Using this approach we generated population data for 66 insertions out of the 77 total ISX insertions present in the MSH22 reference genome assembly. Uneven sequence coverage between insertions and individuals meant that not all insertions could be assembled for each individual. However, the majority of individuals are represented in the majority of datasets: each insertion dataset contained ∼20 lines on average (see Dataset S1 in Dryad: ). Almost all ISX insertions are fixed among strains (65 of 66) and insertion sites are identical between lines, suggesting that independent parallel insertions are unlikely to be present within our dataset. We performed PCR and Sanger sequencing on a subset of these regions and estimate the base-calling error rate of our Illumina contigs to be ∼0.1%.

Consistent with non-allelic gene conversion spreading the TT haplotype, we observe a strong signal of allele-sharing within the sequence region flanking the MSL recognition motif among ISX insertions (). On average, 68.9% of polymorphisms observed within a given insertion are shared among other insertions (though most polymorphisms are shared only between a few elements). The TT haplotype is especially striking in this regard as it is polymorphic in 41% of insertions (, and ). If population subdivision contributes to this excess of allele sharing, we would expect individuals to cluster by allele state at the TT locus, across all polymorphic ISX insertions. Instead, we find that different individuals contribute to the TT polymorphism at each of these ISX insertions (), suggesting that abundant non-allelic gene conversion is the most likely explanation for this observation.

Interestingly, the population frequency of the TT haplotype is similar among insertions that are near each other on the chromosome (permutation test p = 0.018; ). This is consistent with higher gene conversion rates between more closely linked ISX elements generating correlated population frequencies among adjacent elements ().

For each ISX insertion identified within the D. Miranda MSH22 reference genome assembly (alignment shown at left, see also ), we characterized sequence variation across D. Miranda individuals. Richard Schneider Jr Dreamlike Land Rar Download. The TT haplotype (magenta lines) was fixed across individuals at 30% of insertions (see example alignment, top right), polymorphic at 41% of insertions (example shown middle right), and absent at 29% of insertions (bottom right).

Allele sharing between insertions occurs at sites other than the TT haplotype, but these sites tend to be shared across fewer insertions (see heatmap, bottom right). Shows the population alignment across all ISX insertions on the neo-X.

To test if selection has acted to increase the frequency of the TT haplotype in the population, we examined patterns of polymorphisms at GA- and TT-containing ISX elements. The TT haplotype harbors significantly less linked variation than the ancestral GA haplotype, across insertion sites and individuals (haplotype diversity = 0.53 vs 0.81; resampling p. ( A) Haplotype diversity across all ISX sequences. Assembled ISX contigs were combined for all insertions and individuals.

The 25 basepairs flanking each side of the TT region were extracted from a total of 1291 sequences and split into two groups based on whether they contained the TT or GA haplotype. Haplotype diversity was then calculated for each group. The difference between groups is significantly larger than expected by chance (resampling p.

Recent work in a variety of eukaryotes suggests that transposable elements may be major drivers of regulatory evolution (; ). Their high transposition rate and ability to supply ready-to use regulatory elements across the genome implies that they may rapidly wire new genes into regulatory networks (). We recently showed that domesticated TEs contribute to rewiring of the dosage compensation network in D. Miranda, but appear to supply only suboptimal binding sites for the MSL complex ().

Here, we identify a derived haplotype with two mutations that interact epistatically to increase binding affinity for the MSL complex. We show that these fine-tuning mutations spread among independent ISX insertions by non-allelic gene conversion, and through the population by natural selection (). Relative to regulatory elements that evolve in isolation, a family of regulatory motifs dispersed by TEs presents a larger mutational target, and a mutation arising in any element contained within a larger repeat unit (the TE) can spread among copies by non-allelic gene conversion.

Consequently, the rate of evolutionary fine-tuning at such regulatory elements can be greatly accelerated by increasing their effective population size (). Thus, transposable elements can ‘crowd-source’ beneficial mutations to rapidly fine-tune regulatory networks. Shared polymorphism of the TT haplotype among ISX insertions suggests a model where a mutation that refines regulatory activity arose once at a single TE-derived regulatory element, and spread across elements via non-allelic gene conversion. Over evolutionary time, such a mutation spreads in two dimensions: horizontally among TE-derived regulatory elements and vertically through the population, until it is fixed across elements and across individuals.

The TT haplotype is at the midpoint of this process. Across ISX insertions, it is fixed, absent, and polymorphic, in approximately equal proportions. Our transgenic experiments show that each individual T allele actually decreases the binding affinity for the MSL complex relative to the ancestral GA haplotype. Thus, TA or GT haplotypes should be selected against in the population if present on a functional ISX element. Consistent with the deleterious effect of individual T alleles, the TA and GT haplotypes are present on some ISY elements but completely absent from ISX, that is, we find the two T mutations to be in perfect linkage disequilibrium among ISX elements but not ISY ().

While most ISY elements carry the ancestral GA haplotype, a small fraction (0.7% of neo-X ISY insertions) instead carry the derived TT haplotype. It is therefore possible that the TT haplotype was introduced onto the ISX background by non-allelic gene conversion from ISY. Under this scenario, the large family of ISY elements in the D. Miranda genome could be acting as a reservoir of natural variation, where complex mutations can accumulate in the absence of epistasis. Non-allelic gene conversion could then transfer these haplotypes to related repetitive elements (such as ISX).

While many of these haplotypes are likely to be neutral or deleterious, some may be beneficial, as in the case of the TT haplotype. Such a scenario avoids the waiting time for a double mutation, as well as the fitness valley that would have to be traversed if the two mutations were to occur sequentially on the ISX background. To conclude, our findings suggest that TE-dispersed binding motifs follow an evolutionary trajectory that is fundamentally different from those that arise by other means.

The complementary roles of TEs in dispersing regulatory motifs, and gene conversion in spreading subsequent refining mutations, combine to allow for the rapid rewiring and fine-tuning of gene regulatory networks. This process adds a new layer of complexity onto how TEs influence regulatory innovation, as well as a new context in which gene conversion affects genome evolution. Materials and methods. Resequencing data were mapped to version 2.2 of the D. Miranda MSH22 reference assembly (GenBank: AJMI00000000.2) using bowtie2 ().

ISX locations were identified in. Paired-end read alignments were evaluated within 2 kilobase windows flanking each ISX insertion and reads with mapping quality of 20 or greater were extracted along with their mate. The extracted mate pairs were then assembled using IDBA-UD, for each line separately ().

Contigs were aligned using FSA () and visualized with Jalview (). A custom Perl script (available at ) was used to identify sequence variants within the alignments. We also PCR amplified eight of the ISX insertions where the TT haplotype was polymorphic. We confirmed that this polymorphism was present at each of these insertions and estimated the base-calling accuracy of the assemblies by sequencing the PCR products using Sanger technology. We used the QuikChange Lightning site-directed mutagenesis kit from Agilent Technologies (Santa Clara, CA) and the ISX element cloned in to engineer four ISX variants that differed only with respect to the TT haplotype: ISX-GA, ISX-GT, ISX-TA, and ISX-TT.

Each construct was injected by BestGene Inc. (Chino Hills, CA) into D. Melanogaster embryos carrying a RMCE landing site at cytosite 38F1 on chromosome 2L (Bloomington Drosophila Stock Center strain #27388). Transformants were verified by PCR and Sanger sequencing. Quantification of allele-specific binding levels of the MSL complex. Male third instar larvae (∼250 mg) were collected from F1 hybrids between ISX-GA and each of the other three engineered lines: ISX-GT, ISX-TA, and ISX-TT.

Chromatin immunoprecipitation was performed for four biological replicates of each of these lines using the MSL2 d-300 primary antibody from Santa Cruz Biotechnology Inc. (Santa Cruz, CA) and the protocol described in. Primers flanking the ISX MRE region were used to generate heterozygous amplicons from the MSL2 IP and input control.

Sanger chromatograms were used in conjunction with polySNP software () to calculate relative abundance of ISX alleles within the IP and input control amplicons. Abundance of the ‘T’ alleles in the IP amplicons relative to the ancestral G/A alleles was calculated and normalized by the same values from the input control. Permutation and resampling tests.

To determine if the TT frequency among neighboring ISX elements was correlated, we clustered elements within 100 kb of each other and calculated the standard deviation in TT allele frequency within clusters. We then compared these values to 1000 permutations where TT allele frequency was randomly shuffled between ISX locations.

The haplotype diversity and allele frequency spectrum resampling tests were performed by drawing, without replacement, two groups of size 617 and 674, respectively, from the pool of 1291 ISX sequences. The intergroup difference in haplotype diversity, as well as the number of derived variants with frequency of 0.75 or greater, was calculated for each of 1000 replicates and compared to the difference between the TT and GA groups.

ELife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see ). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.

Thank you for sending your work entitled “Transposable elements crowd-source mutations to rapidly fine-tune regulatory networks” for consideration at eLife. Your article has been favorably evaluated by Chris Ponting (Senior editor) and 3 reviewers, one of whom is a member of our Board of Reviewing Editors. One reviewer, Hideki Innan, has agreed to reveal his identity. The Reviewing editor and the other reviewers discussed their comments before we reached this decision, and the Reviewing editor has assembled the following comments to help you prepare a revised submission.

There was general agreement that this is an exciting paper, and a substantial contribution to how our understanding of how regulatory evolution might occur. We also thought you could relatively easily do more to strengthen the case for selection having acted to increase the frequency of the TT haplotype. Simple tests based on the extent of haplotype sharing, the allele frequency distribution, etc., appear not to have been carried out. There is an obvious a priori hypothesis to test here, and multiple loci to test it with. There was also general agreement that your title, albeit catchy, is not professional enough. Perhaps use this phrase in the discussion instead, and use a more descriptive title here.

How about something like “Non-allelic gene conversion between repetitive sequences enable rapid evolutionary change at multiple regulatory sites encoded by transposable elements”? Reviewer #2 minor comments: 1) It is a bit unclear on the origin of the TT holpotype. According to the description in the beginning of the subsection “MRE variation among ISX insertions in D.

Miranda strain MSH22”, it seems that the TT haplotype already existed in the ISY family, and so did the GT and TA haplotypes? But because they were rare, the authors assumed that the TT haplotype is derived? The very origin of ISX was GA, and then TT was transferred from ISX to ISY? Do I understand correctly?

More detailed description would be nice. 2) I found an interesting (but not very emphasized) observation that the two sites are in strong LD in ISX, but not in ISY. A small sub-table or something in might help to convince readers. 3) I like that the authors looked at polymorphism.

What does the frequency spectrum of TT (using all sites shown in ) look like? Is it what expected under a directional selection model? Are there any sites that are not fixed in the population yet, but there is already GA/TT polymorphism? 4) I wonder how the spatial distribution of shared sites across the TE region looks. Does your hypothesis predict an excess of shared sites around the two sites? Maybe so, but I'm not sure.

Once the TT haplotypes was introduced in the ISX family, TT haplotypes may be preferentially transferred within the neo-X. In such a case, gene conversion does not matter I don't know any theory that handles such a complicated case. Reviewer #3 minor comments: In the beginning of the subsection “Non-allelic gene conversion is spreading the TT haplotype among ISX insertions”: Actually the decreased linked variation in the TT haplotype is more likely to be the result of non-allelic gene conversion, rather than selection. Gene conversion may happen across hundreds nucleotides, and this in turn can be responsible of the observed decreased variation, on the scale of course of a hundreds bp around the TT sites.: The top-to-bottom orientation of the sequence in this figure to me is kind of counter intuitive, and is different from the more usual left-to-right orientation in. I think might work better if rotated 90° counterclockwise; the authors can try and see if they also like it better. Author response.

There was general agreement that this is an exciting paper, and a substantial contribution to how our understanding of how regulatory evolution might occur. We also thought you could relatively easily do more to strengthen the case for selection having acted to increase the frequency of the TT haplotype. Simple tests based on the extent of haplotype sharing, the allele frequency distribution, etc., appear not to have been carried out. There is an obvious a priori hypothesis to test here, and multiple loci to test it with. There was also general agreement that your title, albeit catchy, is not professional enough. Perhaps use this phrase in the discussion instead (to encourage popular science writers to pick it up), and use a more descriptive title here.

How about something like “Non-allelic gene conversion between repetitive sequences enable rapid evolutionary change at multiple regulatory sites encoded by transposable elements”? We have strengthened the case for selection having acted to increase the frequency of the TT haplotype, and in addition to haplotype diversity include two additional tests of directional selection using the allele frequency spectrum and nucleotide diversity.

Consistent with models of directional selection, the frequency spectrum of TT-carrying haplotypes shows an excess of high frequency derived alleles, compared to the GA-carrying ISX elements, and nucleotide diversity is reduced at ISX elements where the TT haplotype is fixed, relative to GA-containing insertion sites. The results of the three statistical tests of selection are now presented in the new. We have also modified our title as suggested. Reviewer #2 minor comments: 1) It is a bit unclear on the origin of the TT holpotype. According to the description in the beginning of the subsection “MRE variation among ISX insertions in D.

Miranda strain MSH22”, it seems that the TT haplotype already existed in the ISY family, and so did the GT and TA haplotypes? But because they were rare, the authors assumed that the TT haplotype is derived? The very origin of ISX was GA, and then TT was transferred from ISX to ISY? Do I understand correctly? More detailed description would be nice.

There are multiple possibilities of how the TT haplotype could have originated—either by a double mutation in a GA-containing ISX element, or the two mutations might have occurred subsequently (and the population went through an adaptive valley), or it might have originated in an ISY element, and gene converted onto ISX. We cannot distinguish between these possibilities, but add some discussion on its origin in the Discussion section of the manuscript. 2) I found an interesting (but not very emphasized) observation that the two sites are in strong LD in ISX, but not in ISY. A small sub-table or something in might help to convince readers. We added a second panel to () showing the frequencies of all observed haplotypes for ISY and ISX elements and modified the figure legend to emphasize the strong LD of the TT haplotype on the ISX background, and also in the Discussion section of the paper.

3) I like that the authors looked at polymorphism. What does the frequency spectrum of TT (using all sites shown in ) look like? Is it what expected under a directional selection model? Are there any sites that are not fixed in the population yet, but there is already GA/TT polymorphism?

We have added a supplementary figure (Figure 4–figure supplement 3) that compares the frequency spectrum across TT-carrying ISX sequences versus GA-carrying ISX sequences. Consistent with directional selection models of incomplete hitchhiking due to recombination and/or gene conversion, the frequency spectrum of TT-carrying haplotypes shows an excess of high frequency derived alleles, compared to the GA-carrying ISX elements. 4) I wonder how the spatial distribution of shared sites across the TE region looks. Does your hypothesis predict an excess of shared sites around the two sites? Maybe so, but I'm not sure. Once the TT haplotypes was introduced in the ISX family, TT haplotypes may be preferentially transferred within the neo-X.

In such a case, gene conversion does not matter I don't know any theory that handles such a complicated case. The spatial distribution of shared polymorphism across the ISX element is shown on the bottom right of. We are also not aware of any theory that would handle gene conversion and selection in multigene families, and make predictions about patterns of polymorphism at linked sites.

Reviewer #3 minor comments: In the beginning of the subsection “Non-allelic gene conversion is spreading the TT haplotype among ISX insertions”: Actually the decreased linked variation in the TT haplotype is more likely to be the result of non-allelic gene conversion, rather than selection. Gene conversion may happen across hundreds nucleotides, and this in turn can be responsible of the observed decreased variation, on the scale of course of a hundreds bp around the TT sites. We don't think this is the case. Non-allelic gene conversion can certainly reduce variation among insertions. However, while non-allelic gene conversion may decrease variation within the TT haplotype, it should equally decrease variation within the GA haplotype (assuming that they are subject to similar levels of gene conversion).

The test of selection we employed shows that given its frequency, the TT haplotype harbors significantly less variation than the GA haplotype (a test that is similar in spirit to the haplotype test of selection introduced by Hudson et al, Genetics 1994). We now also show that the frequency spectrum of TT-carrying haplotypes shows an excess of high frequency derived alleles, compared to the GA-carrying ISX elements (Figure 4–figure supplement 3), which further strengthens the case that selection is operating on the TT haplotype.: The top-to-bottom orientation of the sequence in this figure to me is kind of counter intuitive, and is different from the more usual left-to-right orientation in. I think might work better if rotated 90° counterclockwise; the authors can try and see if they also like it better.

Is oriented the same way as, i.e. The nucleotides of individual ISX elements are aligned from left-to-right, and the 77 different ISX insertions within the MSH22 individual are aligned top-to-bottom. We rephrase our figure legend to state this more clearly, and also modify and include the position of the individual ISX elements in the reference MSH22 strain, to make this figure more intuitive (as done in ). Article and author information. ELife is a non-profit organisation inspired by research funders and led by scientists.

Our mission is to help scientists accelerate discovery by operating a platform for research communication that encourages and recognises the most responsible behaviours in science. ELife Sciences Publications, Ltd is a limited liability non-profit non-stock corporation incorporated in the State of Delaware, USA, with company number 5030732, and is registered in the UK with company number FC030576 and branch number BR015634 at the address: eLife Sciences Publications, Ltd 1st Floor, 24 Hills Road Cambridge CB2 1JP UK.

Cover: ARTIST: LIl Wayne ALBUM: Tha Carter III SOURCE: CDDA RIP DATE: 5/30/2008 TOTAL TRACKS: 18 PLAYTIME: 01:26:56 STREET DATE: 5/30/2008 RLS TYPE: Album BITRATE: 166kbps GENRE: Rap CBR/VBR?: VBR LABEL: Universal SIZE: 95.92 MB Tracklist: 1. Carter (Featuring Jay-Z) 3. Got Money (Featuring T-Pain) 5. Comfortable (Featuring Babyface) 6. Phone Home 7. Tie My Hands (Featuring Robin Thicke) 9.

Shoot Me Down 10. Playin With Fire 11. Lollipop (Featuring Static Major) 12.

La La (Featuring Brisco & Busta Rhymes) 13. Nothin On Me (Featuring Juelz Santana & Fabolous) 14.

Let The Beat Build 15. Officer (Featuring Bobby Valentino) 16. Misunderstood 17. Lollipop (Remix) 18. Got Money Feat.

T-Pain & Pitbull (Remix).