Genome-wide analysis of retrogene polymorphisms in Drosophila melanogaster

Gene duplication via retrotransposition has been shown to be an important mechanism in evolution, affecting gene dosage and allowing for the acquisition of new gene functions. Although fixed retrotransposed genes have been found in a variety of species, very little effort has been made to identify retrogene polymorphisms. Here, we examine 37 Illumina-sequenced North American Drosophila melanogaster inbred lines and present the first ever data set and analysis of polymorphic retrogenes in Drosophila. We show that this type of polymorphism is quite common, with any two gametes in the North American population differing in the presence or absence of six retrogenes, accounting for ~13% of gene copy-number heterozygosity. These retrogenes were identified by a straightforward method that can be applied using any type of DNA sequencing data. We also use a variant of this method to conduct a genome-wide scan for intron presence/absence polymorphisms, and show that any two chromosomes in the population likely differ in the presence of multiple introns. We show that these polymorphisms are all in fact deletions rather than intron gain events present in the reference genome. Finally, by leveraging the known location of the parental genes that give rise to the retrogene polymorphisms, we provide direct evidence that natural selection is responsible for the excess of fixations of retrogenes moving off of the X chromosome in Drosophila. Further efforts to identify retrogene and intron presence/absence polymorphisms will undoubtedly improve our understanding of the evolution of gene copy number and gene structure.     

Global regulation of X chromosomal genes by the MSL complex in Drosophila melanogaster

A long-standing model postulates that X-chromosome dosage compensation in Drosophila occurs by twofold up-regulation of the single male X, but previous data cannot exclude an alternative model, in which male autosomes are down-regulated to balance gene expression. To distinguish between the two models, we used RNA interference to deplete Male-Specific Lethal (MSL) complexes from male-like tissue culture cells. We found that expression of many genes from the X chromosome decreased, while expression from the autosomes was largely unchanged. We conclude that the primary role of the MSL complex is to up-regulate the male X chromosome.     

Single nucleotide polymorphism markers for genetic mapping in Drosophila melanogaster

For nearly a century, genetic analysis in Drosophila melanogaster has been a powerful tool for analyzing gene function, yet Drosophila lacks the molecular genetic mapping tools that recently have revolutionized human, mouse, and plant genetics. Here, we describe the systematic characterization of a dense set of molecular markers in Drosophila by using a sequence tagged site-based physical map of the genome. We identify 474 biallelic markers in standard laboratory strains of Drosophila that span the genome. Most of these markers are single nucleotide polymorphisms and sequences for these variants are provided in an accessible format. The average density of the new markers is one per 225 kb on the autosomes and one per megabase on the X chromosome. We include in this survey a set of P-element strains that provide additional use for high-resolution mapping. We show one application of the new markers in a simple set of crosses to map a mutation in the hedgehog gene to an interval of <1 Mb. This new map resource significantly increases the efficiency and resolution of recombination mapping and will be of immediate value to the Drosophila research community.     

Biased clustered substitutions in the human genome: the footprints of male-driven biased gene conversion

We examined fixed substitutions in the human lineage since divergence from the common ancestor with the chimpanzee, and determined what fraction are AT to GC (weak-to-strong). Substitutions that are densely clustered on the chromosomes show a remarkable excess of weak-to-strong "biased" substitutions. These unexpected biased clustered substitutions (UBCS) are common near the telomeres of all autosomes but not the sex chromosomes. Regions of extreme bias are enriched for genes. Human and chimp orthologous regions show a striking similarity in the shape and magnitude of their respective UBCS maps, suggesting a relatively stable force leads to clustered bias. The strong and stable signal near telomeres may have participated in the evolution of isochores. One exception to the UBCS pattern found in all autosomes is chromosome 2, which shows a UBCS peak midchromosome, mapping to the fusion site of two ancestral chromosomes. This provides evidence that the fusion occurred as recently as 740,000 years ago and no more than approximately 3 million years ago. No biased clustering was found in SNPs, suggesting that clusters of biased substitutions are selected from mutations. UBCS is strongly correlated with male (and not female) recombination rates, which explains the lack of UBCS signal on chromosome X. These observations support the hypothesis that biased gene conversion (BGC), specifically in the male germline, played a significant role in the evolution of the human genome.     

Emergence of male-biased genes on the chicken Z-chromosome: sex-chromosome contrasts between male and female heterogametic systems

There has been extensive traffic of male-biased genes out of the mammalian and Drosophila X-chromosomes, and there are also reports of an under-representation of male-biased genes on the X. This may reflect an adaptive process driven by natural selection where an autosomal location of male-biased genes is favored since male genes are only exposed to selection one-third of the time when X-linked. However, there are several alternative explanations to "out-of-the-X" gene movement, including mutational bias and a means for X-linked genes to escape meiotic sex chromosome inactivation (MSCI) during spermatogenesis. As a critical test of the hypothesis that genomic relocation of sex-biased genes is an adaptive process, I examined the emergence, and loss, of genes on the chicken Z-chromosome, i.e., a female heterogametic system (males ZZ, females ZW). Here, the analogous prediction would be an emergence of male-biased genes onto, not a loss from, the Z-chromosome because Z is found more often in males than autosomes are. I found that genes expressed in testis but not in ovary are highly over-represented among genes that have emerged on the Z-chromosome during avian evolution. Moreover, genes with male-biased expression are similarly over-represented among new Z-chromosomal genes. Interestingly, genes with female-biased expression have more often moved from than to the Z-chromosome. These observations show that male and female heterogametic organisms display opposing directionalities in the emergence and loss of sex-biased genes on sex chromosomes. This is consistent with theoretical models on the evolution of sexually antagonistic genes in which new mutations are at least partly dominant.     

FISH mapping of microsatellite loci from <Emphasis Type="Italic">Drosophila subobscura</Emphasis> and its comparison to related species

Microsatellites are highly polymorphic markers that are distributed through all the genome being more abundant in non-coding regions. Whether they are neutral or under selection, these markers if localized can be used as co-dominant molecular markers to explore the dynamics of the evolutionary processes. Their cytological localization can allow identifying genes under selection, inferring recombination from a genomic point of view, or screening for the genomic reorganizations occurring during the evolution of a lineage, among others. In this paper, we report for the first time the localization of microsatellite loci by fluorescent in situ hybridization on Drosophila polytene chromosomes. In Drosophila subobscura, 72 dinucleotide microsatellite loci were localized by fluorescent in situ hybridization yielding unique hybridization signals. In the sex chromosome, microsatellite distribution was not uniform and its density was higher than in autosomes. We identified homologous segments to the sequence flanking the microsatellite loci by browsing the genome sequence of Drosophila pseudoobscura and Drosophila melanogaster. Their localization supports the conservation of Muller’s chromosomal elements among Drosophila species and the existence of multiple intrachromosomal rearrangements within each evolutionary lineage. Finally, the lack of microsatellite repeats in the homologous D. melanogaster sequences suggests convergent evolution for high microsatellite density in the distal part of the X chromosome.     

Retroelements (LINEs and SINEs) in vole genomes: Differential distribution in the constitutive heterochromatin

The chromosomal distribution of mobile genetic elements is scarcely known in Arvicolinae species, but could be of relevance to understand the origin and complex evolution of the sex chromosome heterochromatin. In this work we cloned two retrotransposon sequences, L1 and SINE-B1, from the genome of Chionomys nivalis and investigated their chromosomal distribution on several arvicoline species. Our results demonstrate first that both retroelements are the most abundant repeated DNA sequences in the genome of these species. L1 elements, in most species, are highly accumulated in the sex chromosomes compared to the autosomes. This favoured L1 insertion could have played an important role in the origin of the enlarged heterochromatic blocks existing in the sex chromosomes of some Microtus species. Also, we propose that L1 accumulation on the X heterochromatin could have been the consequence of different, independent and rapid amplification processes acting in each species. SINE elements, however, were completely lacking from the constitutive heterochromatin, either in autosomes or in the heterochromatic blocks of sex chromosomes. These data could indicate that some SINE elements are incompatible with the formation of heterochromatic complexes and hence are necessarily missing from the constitutive heterochromatin.     

Disentangling the relationship between sex-biased gene expression and X-linkage

X chromosomes are preferentially transmitted through females, which may favor the accumulation of X-linked alleles/genes with female-beneficial effects. Numerous studies have shown that genes with sex-biased expression are under- or over-represented on the X chromosomes of a wide variety of organisms. The patterns, however, vary between different animal species, and the causes of these differences are unresolved. Additionally, genes with sex-biased expression tend to be narrowly expressed in a limited number of tissues, and narrowly expressed genes are also non-randomly X-linked in a taxon-specific manner. It is therefore unclear whether the unique gene content of the X chromosome is the result of selection on genes with sex-biased expression, narrowly expressed genes, or some combination of the two. To address this problem, we measured sex-biased expression in multiple Drosophila species and at different developmental time points. These data were combined with available expression measurements from Drosophila melanogaster and mouse to reconcile the inconsistencies in X-chromosome content among taxa. Our results suggest that most of the differences between Drosophila and mammals are confounded by disparate data collection/analysis approaches as well as the correlation between sex bias and expression breadth. Both the Drosophila and mouse X chromosomes harbor an excess of genes with female-biased expression after controlling for the confounding factors, suggesting that the asymmetrical transmission of the X chromosome favors the accumulation of female-beneficial mutations in X-linked genes. However, some taxon-specific patterns remain, and we provide evidence that these are in part a consequence of constraints imposed by the dosage compensation mechanism in Drosophila.     

An analysis of the gene complement of a marsupial, Monodelphis domestica: evolution of lineage-specific genes and giant chromosomes

The newly sequenced genome of Monodelphis domestica not only provides the out-group necessary to better understand our own eutherian lineage, but it enables insights into the innovative biology of metatherians. Here, we compare Monodelphis with Homo sequences from alignments of single nucleotides, genes, and whole chromosomes. Using PhyOP, we have established orthologs in Homo for 82% (15,250) of Monodelphis gene predictions. Those with single orthologs in each species exhibited a high median synonymous substitution rate (d(S) = 1.02), thereby explaining the relative paucity of aligned regions outside of coding sequences. Orthology assignments were used to construct a synteny map that illustrates the considerable fragmentation of Monodelphis and Homo karyotypes since their therian last common ancestor. Fifteen percent of Monodelphis genes are predicted, from their low divergence at synonymous sites, to have been duplicated in the metatherian lineage. The majority of Monodelphis-specific genes possess predicted roles in chemosensation, reproduction, adaptation to specific diets, and immunity. Using alignments of Monodelphis genes to sequences from either Homo or Trichosurus vulpecula (an Australian marsupial), we show that metatherian X chromosomes have elevated silent substitution rates and high G+C contents in comparison with both metatherian autosomes and eutherian chromosomes. Each of these elevations is also a feature of subtelomeric chromosomal regions. We attribute these observations to high rates of female-specific recombination near the chromosomal ends and within the X chromosome, which act to sustain or increase G+C levels by biased gene conversion. In particular, we propose that the higher G+C content of the Monodelphis X chromosome is a direct consequence of its small size relative to the giant autosomes.     

Autosomal location of a new subtype of 1.688 satellite DNA ofDrosophila melanogaster

During the screening of aDrosophila melanogaster YAC library with DNA from the minichromosomeDp(1;f)1187 we isolated a clone, yw20D5, which contains a new subtype of 1.688 satellite DNA. Although the sequences of several monomers subcloned from the YAC show a considerable variation in length, the derived consensus sequence is 356-bp long. This new subtype and the one constituted by the 353-bp repeats are both located on the left arm heterochromatin of chromosome 3, arranged in separate arrays. Despite their autosomal location, phylogenetic relationships among 1.688 satellite sequences suggest that they may have originated from the 359-bp repeats of the X chromosome heterochromatin. We have used the new 356-bp repeats to investigate whether sequences related to the 1.688 satellite are dispersed along the euchromatic arms of the autosomes in a similar way to that in which they are found along the X chromosome euchromatin.accepted for publication by D. Ward     

On the origin of new genes in Drosophila

Several mechanisms have been proposed to account for the origination of new genes. Despite extensive case studies, the general principles governing this fundamental process are still unclear at the whole-genome level. Here, we unveil genome-wide patterns for the mutational mechanisms leading to new genes and their subsequent lineage-specific evolution at different time nodes in the Drosophila melanogaster species subgroup. We find that (1) tandem gene duplication has generated approximately 80% of the nascent duplicates that are limited to single species (D. melanogaster or Drosophila yakuba); (2) the most abundant new genes shared by multiple species (44.1%) are dispersed duplicates, and are more likely to be retained and be functional; (3) de novo gene origination from noncoding sequences plays an unexpectedly important role during the origin of new genes, and is responsible for 11.9% of the new genes; (4) retroposition is also an important mechanism, and had generated approximately 10% of the new genes; (5) approximately 30% of the new genes in the D. melanogaster species complex recruited various genomic sequences and formed chimeric gene structures, suggesting structure innovation as an important way to help fixation of new genes; and (6) the rate of the origin of new functional genes is estimated to be five to 11 genes per million years in the D. melanogaster subgroup. Finally, we survey gene frequencies among 19 globally derived strains for D. melanogaster-specific new genes and reveal that 44.4% of them show copy number polymorphisms within a population. In conclusion, we provide a panoramic picture for the origin of new genes in Drosophila species.     

Does the human X contain a third evolutionary block? Origin of genes on human Xp11 and Xq28

Comparative gene mapping of human X-borne genes in marsupials defined an ancient conserved region and a recently added region of the eutherian X, and the separate evolutionary origins of these regions was confirmed by their locations on chicken chromosomes 4p and 1q, respectively. However, two groups of genes, from the pericentric region of the short arm of the human X (at Xp11) and a large group of genes from human Xq28, were thought to be part of a third evolutionary block, being located in a single region in fish, but mapping to chicken chromosomes other than 4p and 1q. We tested this hypothesis by comparative mapping of genes in these regions. Our gene mapping results show that human Xp11 genes are located on the marsupial X chromosome and platypus chromosome 6, indicating that the Xp11 region was part of original therian X chromosome. We investigated the evolutionary origin of genes from human Xp11 and Xq28, finding that chicken paralogs of human Xp11 and Xq28 genes had been misidentified as orthologs, and their true orthologs are represented in the chicken EST database, but not in the current chicken genome assembly. This completely undermines the evidence supporting a separate evolutionary origin for this region of the human X chromosome, and we conclude, instead, that it was part of the ancient autosome, which became the conserved region of the therian X chromosome 166 million years ago.     

Assembly of functionally active Drosophila origin recognition complex from recombinant proteins

In eukaryotes the sites for the initiation of chromosomal DNA replication are believed to be determined in part by the binding of a heteromeric origin recognition complex (ORC) to DNA. We have cloned the genes encoding the subunits of the Drosophila ORC. Each of the genes is unique and can be mapped to discrete chromosomal locations implying that the pattern and developmental regulation of origin usage in Drosophila is not regulated solely by a large family of different ORC proteins. The six-subunit ORC can be reconstituted with recombinant proteins into a complex that restores DNA replication in ORC-depleted Drosophila or Xenopus egg extracts.