Conservation of gene order, structure and sequence between three closely linked genes in Drosophila melanogaster and Drosophila virilis

In Drosophila melanogaster, the apparently unrelated genes anon-66Da, RpL14, and anon-66Db (from telomere to centromere) are located on a 5547 bp genomic fragment on chromosome arm 3L at cytological position 66D8. The three genes are tightly linked, and flanked by two relatively large genes with unknown function. We have taken a comparative genomic approach to investigate the evolutionary history of the three genes. To this end we isolated a Drosophila virilis 7.3 kb genomic fragment which is homologous to a 5.5 kb genomic region of D. melanogaster. Both fragments map to Muller's element D, namely to section 66D in D. melanogaster and to section 32E in D. virilis, and harbor the genes anon-66Da, RpL14, and anon-66Db. We demonstrate that the three genes exhibit a high conservation of gene topography in general and in detail. While most introns and intergenic regions reveal sequence divergences, there are, however, a number of interspersed conserved sequence motifs. In particular, two introns of the RpL14 gene contain a short, highly conserved 60 nt long sequence located at corresponding positions. This sequence represents a novel Drosophila small nucleolar RNA, which is homologous to human U49. Whereas DNA flanking the three genes shows no significant interspecies homologies, the 3'-flanking region in D. virilis contains sequences from the transposable element Penelope. The Penelope family of transposable elements has been shown to promote chromosomal rearrangements in the D. virilis species group. The presence of Penelope sequences in the D. virilis 7.3 kb genomic fragment may be indicative for a transposon-induced event of transposition which did not yet scramble the order of the three genes but led to the breakdown of sequence identity of the flanking DNA.     

Functional and comparative genomic analysis of the piebald deletion region of mouse chromosome 14

Several developmentally important genomic regions map within the piebald deletion complex on distal mouse chromosome 14. We have combined computational gene prediction and comparative sequence analysis to characterize an approximately 4.3-Mb segment of the piebald region to identify candidate genes for the phenotypes presented by homozygous deletion mice. As a result we have ordered 13 deletion breakpoints, integrated the sequence with markers from a bacterial artificial chromosome (BAC) physical map, and identified 16 known or predicted genes and >1500 conserved sequence elements (CSEs) across the region. The candidate genes identified include Phr1 (formerly Pam) and Spry2, which are mouse homologs of genes required for development in Drosophila melanogaster. Gene content, order, and position are highly conserved between mouse chromosome 14 and the orthologous region of human chromosome 13. Our studies combining computational gene prediction with genetic and comparative genomic analyses provide insight regarding the functional composition and organization of this defined chromosomal region.     

Using the P[wHy] hybrid transposable element to disrupt genes in region 54D-55B in Drosophila melanogaster

Understanding the function of each gene in the genome of a model organism such as Drosophila melanogaster is an important goal. The development of improved methods for uncovering the mutant phenotypes of specific genes can accelerate achievement of this goal. The P[wHy] hybrid transposable element can be used to generate nested sets of precisely mapped deletions in a given region of the Drosophila genome. Here we use the P[wHy] method to generate overlapping, molecularly defined deletions from a set of three P[wHy] insertions in the 54E-F region of chromosome 2. Deletions that span a total of 0.5 Mb were identified and molecularly mapped precisely. Using overlapping deletions, the mutant phenotypes of nine previously uncharacterized genes in a 101-kb region were determined, including identification of new loci required for viability and female fertility. In addition, the deletions were used to molecularly map previously isolated lethal mutations. Thus, the P[wHy] method provides an efficient method for systematically determining the phenotypes of genes in a given region of the fly genome.     

Colonization of heterochromatic genes by transposable elements in Drosophila

As a further step toward understanding transposable element-host genome interactions, we investigated the molecular anatomy of introns from five heterochromatic and 22 euchromatic protein-coding genes of Drosophila melanogaster. A total of 79 kb of intronic sequences from heterochromatic genes and 355 kb of intronic sequences from euchromatic genes have been used in Blast searches against Drosophila transposable elements (TEs). The results show that TE-homologous sequences belonging to 19 different families represent about 50% of intronic DNA from heterochromatic genes. In contrast, only 0.1% of the euchromatic intron DNA exhibits homology to known TEs. Intraspecific and interspecific size polymorphisms of introns were found, which are likely to be associated with changes in TE-related sequences. Together, the enrichment in TEs and the apparent dynamic state of heterochromatic introns suggest that TEs contribute significantly to the evolution of genes located in heterochromatin.     

A Physical Map of the X Chromosome of Drosophila Melanogaster: Cosmid Contigs and Sequence Tagged Sites

A physical map of the euchromatic X chromosome of Drosophila melanogaster has been constructed by assembling contiguous arrays of cosmids that were selected by screening a library with DNA isolated from microamplified chromosomal divisions. This map, consisting of 893 cosmids, covers ~64% of the euchromatic part of the chromosome. In addition, 568 sequence tagged sites (STS), in aggregate representing 120 kb of sequenced DNA, were derived from selected cosmids. Most of these STSs, spaced at an average distance of ~35 kb along the euchromatic region of the chromosome, represent DNA tags that can be used as entry points to the fruitfly genome. Furthermore, 42 genes have been placed on the physical map, either through the hybridization of specific probes to the cosmids or through the fact that they were represented among the STSs. These provide a link between the physical and the genetic maps of D. melanogaster. Nine novel genes have been tentatively identified in Drosophila on the basis of matches between STS sequences and sequences from other species.     

Insertion-deletion variation at the yellow-achaete-scute region in two natural populations of Drosophila melanogaster

We have surveyed the region of the X chromosome of Drosophila melanogaster which encodes the yellow, achaete and scute genes for restriction map variation. Two natural populations, one from North Carolina, U.S.A. and the other from southern Spain were screened for variation at about 70 restriction sites and for variation due to DNA insertion or deletion events in 120 kilobases of DNA. Mean heterozygosity per nucleotide was estimated to be 0.0024 and 15 large insertions were found in the 49 chromosomes screened. Extensive disequilibrium between polymorphic sites were found across much of the region in the North Carolina population. The frequency of large insertions, which usually correspond to transposable genetic elements, is significantly lower than has been observed in autosomal regions of the genome. This is predicted for X-linked loci by certain models of transposable element evolution, where copy number is restricted by virtue of the recessive deleterious effects of the insertions. Our results appear to support such models. The deficiency of insertions may in this case be enhanced by hitch-hiking effects arising from the high level of disequilibrium.     

Genetic and bioinformatic analysis of 41C and the 2R heterochromatin of Drosophila melanogaster: a window on the heterochromatin-euchromatin junction

Genomic sequences provide powerful new tools in genetic analysis, making it possible to combine classical genetics with genomics to characterize the genes in a particular chromosome region. These approaches have been applied successfully to the euchromatin, but analysis of the heterochromatin has lagged somewhat behind. We describe a combined genetic and bioinformatics approach to the base of the right arm of the Drosophila melanogaster second chromosome, at the boundary between pericentric heterochromatin and euchromatin. We used resources provided by the genome project to derive a physical map of the region, examine gene density, and estimate the number of potential genes. We also carried out a large-scale genetic screen for lethal mutations in the region. We identified new alleles of the known essential genes and also identified mutations in 21 novel loci. Fourteen complementation groups map proximal to the assembled sequence. We used PCR to map the endpoints of several deficiencies and used the same set of deficiencies to order the essential genes, correlating the genetic and physical map. This allowed us to assign two of the complementation groups to particular "computed/curated genes" (CGs), one of which is Nipped-A, which our evidence suggests encodes Drosophila Tra1/TRRAP.     

Molecular and bioinformatic analysis of the FB-NOF transposable element

The Drosophila melanogaster transposable element FB-NOF is known to play a role in genome plasticity through the generation of all sort of genomic rearrangements. Moreover, several insertional mutants due to FB mobilizations have been reported. Its structure and sequence, however, have been poorly studied mainly as a consequence of the long, complex and repetitive sequence of FB inverted repeats. This repetitive region is composed of several 154 bp blocks, each with five almost identical repeats. In this paper, we report the sequencing process of 2 kb long FB inverted repeats of a complete FB-NOF element, with high precision and reliability. This achievement has been possible using a new map of the FB repetitive region, which identifies unambiguously each repeat with new features that can be used as landmarks. With this new vision of the element, a list of FB-NOF in the D. melanogaster genomic clones has been done, improving previous works that used only bioinformatic algorithms. The availability of many FB and FB-NOF sequences allowed an analysis of the FB insertion sequences that showed no sequence specificity, but a preference for A/T rich sequences. The position of NOF into FB is also studied, revealing that it is always located after a second repeat in a random block. With the results of this analysis, we propose a model of transposition in which NOF jumps from FB to FB, using an unidentified transposase enzyme that should specifically recognize the second repeat end of the FB blocks.     

Evolutionary rearrangement of the amylase genomic regions between Drosophila melanogaster and Drosophila pseudoobscura

Two Drosophila pseudoobscura genomic clones have sequence similarity to the Drosophila melanogaster amylase region that maps to the 53CD region on the D. melanogaster cytogenetic map. The two clones with similarity to amylase map to sections 73A and 78C of the D. pseudoobscura third chromosome cytogenetic map. The complete sequences of both the 73A and 78C regions were compared to the D. melanogaster genome to determine if the coding region for amylase is present in both regions and to determine the evolutionary mechanism responsible for the observed distribution of the amylase gene or genes. The D. pseudoobscura 73A and 78C linkage groups are conserved with the D. melanogaster 41E and 53CD regions, respectively. The amylase gene, however, has not maintained its conserved linkage between the two species. These data indicate that amylase has moved via a transposition event in the D. melanogaster or D. pseudoobscura lineage. The predicted genes within the 73A and 78C regions show patterns of molecular evolution in synonymous and nonsynonymous sites that are consistent with previous studies of these two species.     

A complete and a truncated U1 snRNA gene of Drosophila melanogaster are found as inverted repeats at region 82E of the polytene chromosomes.

A phage containing two sequences homologous to U1 snRNA was isolated from a Drosophila melanogaster genomic library, and identified with a previously cloned D. melanogaster U1 snRNA gene. DNA sequence analysis showed that complete and truncated U1 snRNA genes are present, both of which have base substitutions relative to U1 snRNA. These genes show conservation of 5' and 3' flanking regions relative to other U1 and U2 snRNA genes of Drosophila. Intramolecular renaturation experiments and electron microscope mapping demonstrates that the two U1 snRNA sequences are present as inverted repeats about 2.7kb apart, separated by a smaller pair of inverted repeats of an unrelated sequence. These U1 snRNA sequences were located by in situ hybridization at 82E, and related sequences were found at 21D and 95C on the polytene chromosome map. The results are discussed with reference to the origin and function of snRNAs.     

Recurrent adaptation in RNA interference genes across the Drosophila phylogeny

RNA interference (RNAi) is quickly emerging as a vital component of genome organization, gene regulation, and immunity in Drosophila and other species. Previous studies have suggested that, as a whole, genes involved in RNAi are under intense positive selection in Drosophila melanogaster. Here, we characterize the extent and patterns of adaptive evolution in 23 known Drosophila RNAi genes, both within D. melanogaster and across the Drosophila phylogeny. We find strong evidence for recurrent protein-coding adaptation at a large number of RNAi genes, particularly those involved in antiviral immunity and defense against transposable elements. We identify specific functional domains involved in direct protein-RNA interactions as particular hotspots of recurrent adaptation in multiple RNAi genes, suggesting that targeted coadaptive arms races may be a general feature of RNAi evolution. Our observations suggest a predictive model of how selective pressures generated by evolutionary arms race scenarios may affect multiple genes across protein interaction networks and other biochemical pathways.     

Coding region deletions associated with the major form of rDNA interruption in Drosophila

The nucleotide sequences at and around the termini of 5 kb type 1 interruptions in three separate clones of D. melanogaster rDNA repeats have been determined, and have been compared with the sequence of the corresponding region of an insertion-free rDNA repeat. All three interrupted rDNA repeats contain a small deletion of 28S rRNA coding material at the left coding/insertion sequence junction. A second deletion was found in one of the three clones, ad other aberrations were suggested by the results of restriction enzyme digestions of unfractionated rDNA. The termini of 5 kb type 1 rDNA insertions in D. melanogaster were also compared with the corresponding regions of 28S rDNA interruptions in D. virilis: the insertion site is identical in the two species, but the termini of the two species' interruptions show no homology. I sequenced a 1.1 kb region of the 5 kb type 1 D. melanogaster rDNA interruption that covers the sequences of the 1 kb and 0.5 kb insertions. There is 98% homology between the rightmost 1 kb of the 5 kb interruption and the sequences of the shorter insertions. Data suggest that Drosophila rDNA interruptions arose as a transposable element, and that divergence had included length alterations generated by unequal crossing over.