Whole genome resequencing reveals natural target site preferences of transposable elements in Drosophila melanogaster

Transposable elements are mobile DNA sequences that integrate into host genomes using diverse mechanisms with varying degrees of target site specificity. While the target site preferences of some engineered transposable elements are well studied, the natural target preferences of most transposable elements are poorly characterized. Using population genomic resequencing data from 166 strains of Drosophila melanogaster, we identified over 8,000 new insertion sites not present in the reference genome sequence that we used to decode the natural target preferences of 22 families of transposable element in this species. We found that terminal inverted repeat transposon and long terminal repeat retrotransposon families present clade-specific target site duplications and target site sequence motifs. Additionally, we found that the sequence motifs at transposable element target sites are always palindromes that extend beyond the target site duplication. Our results demonstrate the utility of population genomics data for high-throughput inference of transposable element targeting preferences in the wild and establish general rules for terminal inverted repeat transposon and long terminal repeat retrotransposon target site selection in eukaryotic genomes.     

Two distinct domains in Drosophila melanogaster telomeres

Telomeres are generally considered heterochromatic. On the basis of DNA composition, the telomeric region of Drosophila melanogaster contains two distinct subdomains: a subtelomeric region of repetitive DNA, termed TAS, and a terminal array of retrotransposons, which perform the elongation function instead of telomerase. We have identified several P-element insertions into this retrotransposon array and compared expression levels of transgenes with similar integrations into TAS and euchromatic regions. In contrast to insertions in TAS, which are silenced, reporter genes in the terminal HeT-A, TAHRE, or TART retroelements did not exhibit repressed expression in comparison with the same transgene construct in euchromatin. These data, in combination with cytological studies, provide evidence that the subtelomeric TAS region exhibits features resembling heterochromatin, while the terminal retrotransposon array exhibits euchromatic characteristics.     

Fixation of transposable elements in the Drosophila melanogaster genome

We have investigated at the molecular level four cases in which D. melanogaster middle repetitive DNA probes consistently hybridized to a particular band on chromosomes sampled from a D. melanogaster natural population. Two corresponded to true fixations of a roo and a Stalker element, and the others were artefacts of the in situ hybridization technique caused by the presence of genomic DNA flanking the transposable elements (TEs) in the probes. The two fixed elements are located in the beta-heterochromatin (20A and 80B, respectively) and are embedded in large clusters of other elements, many of which may also be fixed. We also found evidence that this accumulation is an ongoing process. These results support the hypothesis that TEs accumulate in the non-recombining part of the genome. Their implications for the effects of TEs on determining the chromatin structure of the host genomes are discussed in the light of recent evidence for the role of TE-derived small interfering-RNAs as cis -acting determinants of heterochromatin formation.     

Principles of genome evolution in the Drosophila melanogaster species group

That closely related species often differ by chromosomal inversions was discovered by Sturtevant and Plunkett in 1926. Our knowledge of how these inversions originate is still very limited, although a prevailing view is that they are facilitated by ectopic recombination events between inverted repetitive sequences. The availability of genome sequences of related species now allows us to study in detail the mechanisms that generate interspecific inversions. We have analyzed the breakpoint regions of the 29 inversions that differentiate the chromosomes of Drosophila melanogaster and two closely related species, D. simulans and D. yakuba, and reconstructed the molecular events that underlie their origin. Experimental and computational analysis revealed that the breakpoint regions of 59% of the inversions (17/29) are associated with inverted duplications of genes or other nonrepetitive sequences. In only two cases do we find evidence for inverted repetitive sequences in inversion breakpoints. We propose that the presence of inverted duplications associated with inversion breakpoint regions is the result of staggered breaks, either isochromatid or chromatid, and that this, rather than ectopic exchange between inverted repetitive sequences, is the prevalent mechanism for the generation of inversions in the melanogaster species group. Outgroup analysis also revealed evidence for widespread breakpoint recycling. Lastly, we have found that expression domains in D. melanogaster may be disrupted in D. yakuba, bringing into question their potential adaptive significance.     

Distribution of dinucleotide microsatellites in the Drosophila melanogaster genome.

Microsatellites, a special class of repetitive DNA, have become one of the most popular genetic markers. The progress of various genome projects has made it possible to study the genomic distribution of microsatellites and to evaluate the potential influence of several parameters on their genesis. We report the distribution of dinucleotide microsatellites in the genome of Drosophila melanogaster. When considering only microsatellites with five or more repeat units, the average length of dinucleotide repeats in D. melanogaster is 6.7 repeats. We tested a wide range of parameters which could potentially influence microsatellite density, and we did not detect a significant influence of recombination rate, number of exons, or total length of coding sequence. In concordance with the neutral expectation for the origin of microsatellites, a significant positive correlation between AT content and (AT/TA)n microsatellite density was detected. While this pattern may indicate that microsatellite genesis is a random process, we also found evidence for a nonrandom distribution of microsatellites. Average microsatellite density was higher on the X chromosome, but extreme heterogeneity was observed between different genomic regions. Such a clumping of microsatellites was also evident on a more local scale, as 38.9% of the contiguous sequences analyzed showed a deviation from a random distribution of microsatellites.     

Mobilization of retrotransposon copia in the Drosophila melanogaster genome]

Considerable heterogeneity of retrotransposon copia sites of location on polytene chromosomes was revealed in one of the substocks of the inbred Drosophila melanogaster stock. Heterogeneity of copia sites of location was found in no other substocks analyzed. The heterogeneity was shown to be caused by copia insertions in new sites. The frequency of insertions is about 12% per haploid genome per generation. The retrotransposon excisions and somatic transpositions were not observed. The location of retrotransposons mdg1, mdg2, mdg3, mdg4, 297 and H.M.S. Beagle appeared to be stable in all the stocks analyzed. Thus, a model system allowing to study mechanisms of retrotransposon copia transpositions in D. melanogaster tissues as well as phenotypic effects of copia mobilization is described.     

Cytochemistry of late ovarian chambers of Drosophila melanogaster

Summary Late ovarian chambers of Drosophila melanogaster have been examined by ultrastructural cytochemistry in an attempt to characterize some of the transformations which precede the completion of oogenesis. From stage 11 onward peroxidase activity is present in the endoplasmic reticulum of both nurse cells and oocyte, as well as in the egg-covering precursors of the columnar follicle cells. Catalase activity is restricted to the very last stages of oogenesis (stage 13–14) and appears to be located in membrane-bound organelles of the ooplasm which are continuous with the endoplasmic reticulum. Because of the presence of catalase as well as by their structural appearance, these organelles are to be identified as microperoxisomes. Catalase activity becomes cytochemically detectable in the ooplasm somehow in coincidence with the formation of glycogen. Furthermore, glycogen is first formed in intimate association with alpha-1 yolk platelets. On the basis of these findings it is suggested that glycogen synthesis occurs by a process of gluconeogenesis.     

Multiple Allelism of the net Gene in Drosophila melanogaster and Drosophila simulans

Thenetgene mutations are known to cause abnormal pattern of veining in all wing regions except for the first posterior cells. In natural populations of Drosophila melanogaster, the net alleles were identified, which differ in phenotypic expression from standard mutations. The mutants net-extra-analis from a population Belokurikha-2000 have only a single additional vein in the third posterior cell. A line from Chernobyl-1986 population have another nontypical allele net ^Ch86 and shows a lower degree of abnormalities than that usually observed. About 10% of these flies have an additional vein fragment in the first posterior cell. In both males and females ofD. simulans population Tashkent -2001, which exhibit net ^ST91 mutation, a net of additional veins is formed as a specific additional fragment in the first posterior cell. The pattern of veining conferred by alleles net-extra-analis and net ^Ch86 is altered to a lesser extent; these alleles are dominant with respect to alleles net ^2-45 and net ^ST91, which cause more abnormalities. The heterozygotes for alleles net ^ST9 and net ^Ch86 and for Df(2) net ⁶² deletion have an additional fragment in the first posterior cell and show similarly strong deviations from normal wing vein pattern. The naturalnet alleles correspond, presumably, to different molecular gene defects involved into uncertain local interactions with numerous modifying factors and other genes that specify the wing vein pattern.     

Transposable elements as artisans of the heterochromatic genome in Drosophila melanogaster

Over 50 years ago Barbara McClintock discovered that maize contains mobile genetic elements, but her findings were at first considered nothing more than anomalies. Today it is widely recognized that transposable elements have colonized all eukaryotic genomes and represent a major force driving evolution of organisms. Our contribution to this special issue deals with the theme of transposable element-host genome interactions. We bring together published and unpublished work to provide a picture of the contribution of transposable elements to the evolution of the heterochromatic genome in Drosophila melanogaster. In particular, we discuss data on 1) colonization of constitutive heterochromatin by transposable elements, 2) instability of constitutive heterochromatin induced by the I factor, and 3) evolution of constitutive heterochromatin and heterochromatic genes driven by transposable elements. Drawing attention to these topics may have direct implications on important aspects of genome organization and gene expression.     

Computational identification of novel chitinase-like proteins in the Drosophila melanogaster genome

MOTIVATION: Multiple chitinases as well as lectins closely related to them have been characterized previously from many insect species and the corresponding genes/cDNAs have been cloned. However, the identification of the entire assortment of genes for chitinase family proteins and their differences in biochemical properties have not been carried out in any individual insect species. The completion of the entire DNA sequence of Drosophila melanogaster (fruit fly) genome and identification of open reading frames presents an opportunity to study the structures and functions of chitinase-like proteins, and also to identify new members of this family in DROSOPHILA: We are, therefore, interested in studying the functional genomics of chitinase-like gene families in insects. METHODS: We searched the Drosophila protein sequences database using fully characterized insect chitinase sequences and BLASTP software, identified all the putative chitinase-like proteins encoded in Drosophila genome, and predicted their structures using domain analysis tools. A phylogenetic analysis of the chitinase-like proteins from Drosophila and several other insect species was carried out. The structures of these chitinases were modeled using homology modeling software. RESULTS: Our analysis revealed the presence of 18 chitinase-like proteins in the Drosophila protein database. Among these are seven novel chitinase-like proteins that contain four signature amino acid sequences of chitinases belonging to family 18 glycosylhydrolases, including both acidic and hydrophobic amino acid residues critical for enzyme activity. All the proteins contain at least one catalytic domain with one having four catalytic domains. Phylogenetic analysis of chitinase-like proteins from Drosophila and other insects revealed an evolutionary relationship among all these proteins, which indicated gene duplication and domain shuffling to generate the observed diversity in the encoded proteins. Homology modeling showed that all the Drosophila chitinase-like proteins contain one or more catalytic domains with a (alpha/beta)8 barrel-like structure. Our results suggest that insects utilize multiple family 18 chitinolytic enzymes and also non-enzymatic chitinase-like proteins for degrading/remodeling/binding to chitin in different insect anatomical extracellular structures, such as the cuticle, peritrophic membrane, trachea and mouth parts during insect development, and possibly for other roles including chitin synthesis. AVAILABILITY: Perl program and supplementary material are available at http://www.ksu.edu/bioinformatics/supplementary.htm     

In silico identification of novel selenoproteins in the Drosophila melanogaster genome

In selenoproteins, incorporation of the amino acid selenocysteine is specified by the UGA codon, usually a stop signal. The alternative decoding of UGA is conferred by an mRNA structure, the SECIS element, located in the 3′-untranslated region of the selenoprotein mRNA. Because of the non-standard use of the UGA codon, current computational gene prediction methods are unable to identify selenoproteins in the sequence of the eukaryotic genomes. Here we describe a method to predict selenoproteins in genomic sequences, which relies on the prediction of SECIS elements in coordination with the prediction of genes in which the strong codon bias characteristic of protein coding regions extends beyond a TGA codon interrupting the open reading frame. We applied the method to the Drosophila melanogaster genome, and predicted four potential selenoprotein genes. One of them belongs to a known family of selenoproteins, and we have tested experimentally two other predictions with positive results. Finally, we have characterized the expression pattern of these two novel selenoprotein genes.