Comparison of dot chromosome sequences from D. melanogaster and D. virilis reveals an enrichment of DNA transposon sequences in heterochromatic domains

Background

Chromosome four of Drosophila melanogaster, known as the dot chromosome, is largely heterochromatic, as shown by immunofluorescent staining with antibodies to heterochromatin protein 1 (HP1) and histone H3K9me. In contrast, the absence of HP1 and H3K9me from the dot chromosome in D. virilis suggests that this region is euchromatic. D. virilis diverged from D. melanogaster 40 to 60 million years ago.

Results

Here we describe finished sequencing and analysis of 11 fosmids hybridizing to the dot chromosome of D. virilis (372,650 base-pairs) and seven fosmids from major euchromatic chromosome arms (273,110 base-pairs). Most genes from the dot chromosome of D. melanogaster remain on the dot chromosome in D. virilis, but many inversions have occurred. The dot chromosomes of both species are similar to the major chromosome arms in gene density and coding density, but the dot chromosome genes of both species have larger introns. The D. virilis dot chromosome fosmids have a high repeat density (22.8%), similar to homologous regions of D. melanogaster (26.5%). There are, however, major differences in the representation of repetitive elements. Remnants of DNA transposons make up only 6.3% of the D. virilis dot chromosome fosmids, but 18.4% of the homologous regions from D. melanogaster; DINE-1 and 1360 elements are particularly enriched in D. melanogaster. Euchromatic domains on the major chromosomes in both species have very few DNA transposons (less than 0.4 %).

Conclusion

Combining these results with recent findings about RNAi, we suggest that specific repetitive elements, as well as density, play a role in determining higher-order chromatin packaging.     

Satellite Ic: a possible link between the satellite DNAs of D. virilis and D. melanogaster.

In this study, we isolated and characterized a previously undetected cryptic satellite DNA comprising 0.1% of the total nuclear genome of D. virilis. This satellite is hidden from detection in neutral CsCl by satellite I and is therefore designated cryptic satellite I or Ic. Sequence analysis reveals that Ic is the repeating heptanucleotide [poly d(AATATAG): d(CTATATT)]. It is more closely related to the three simple sequence satellite DNAs of D. melanogaster, a distantly related species, than it is to any of the major D. virilis satellite DNA sequences. Ic may therefore be a link between the simple sequence satellites of D. virilis and D. melanogaster. As an extension of this theory, we have constructed a "family tree" linking the satellites of D. virilis and D. melanogaster by a series of "simple" operations. Only one intermediate required by this evolutionary scheme has not yet been identified.     

A comparative study of male meiosis in Drosophila melanogaster and D. virilis

Male meiosis in D. melanogaster cytologically follows the usual pattern, whereas in D. melanogaster and in D. virilis oocytes the chromosomes clump into a karyosphere at early meiotic prophase and remain so up to metaphase I.Male meiosis in D. virilis spermatocytes has an intermediate character: a part of the chromatin clumps together in a karyosphere at early prophase, whereas the other part of the chromatin remains diffuse all through prophase. At the end of prophase, the diffuse chromatin becomes integrated into the karyosphere before metaphase I. During the meiotic divisions the chromosomes have the same clumped aspect as those in Drosophila oocytes and thus differ strikingly from the dividing chromosomes in D. melanogaster spermatocytes.In D. virilis spermatocytes the nucleolus exhibits changes during the meiotic prophase that may be related to synthetical activities. The DNA specific staining with the fluorochrome DAPI reveals the existence of extrachromosomal DNA in the later prophase. Other striking differences in meiotic events between the two Drosophila species concern the centrioles and spermiogenesis.     

Comparative Analysis of Pdf-Mediated Circadian Behaviors Between Drosophila melanogaster and D. virilis

PAR proteins (partitioning defective) are major regulators of cell polarity and asymmetric cell division. One of the par genes, par-1, encodes a Ser/Thr kinase that is conserved from yeast to mammals. In Caenorhabditis elegans, par-1 governs asymmetric cell division by ensuring the polar distribution of cell fate determinants. However the precise mechanisms by which PAR-1 regulates asymmetric cell division in C. elegans remain to be elucidated. We performed a genomewide RNAi screen and identified six genes that specifically suppress the embryonic lethal phenotype associated with mutations in par-1. One of these suppressors is mpk-1, the C. elegans homolog of the conserved mitogen activated protein (MAP) kinase ERK. Loss of function of mpk-1 restored embryonic viability, asynchronous cell divisions, the asymmetric distribution of cell fate specification markers, and the distribution of PAR-1 protein in par-1 mutant embryos, indicating that this genetic interaction is functionally relevant for embryonic development. Furthermore, disrupting the function of other components of the MAPK signaling pathway resulted in suppression of par-1 embryonic lethality. Our data therefore indicates that MAP kinase signaling antagonizes PAR-1 signaling during early C. elegans embryonic polarization.ASYMMETRIC cell division, a process in which a mother cell divides in two different daughter cells, is a fundamental mechanism to achieve cell diversity during development. We use the early embryo of Caenorhabditis elegans as a model system to study asymmetric cell division. The C. elegans one-cell embryo divides asymmetrically along its anteroposterior axis, generating two cells of different sizes and fates: the larger anterior daughter cell will generate somatic tissues while the smaller posterior daughter cell will generate the germline (Sulston et al. 1983).A group of proteins called PAR proteins (partitioning defective) is required for asymmetric cell division in C. elegans (Kemphues et al. 1988). Depletion of any of the seven par genes (par-1 to -6 and pkc-3) leads to defects in asymmetric cell division and embryonic lethality (Kemphues et al. 1988; Kirby et al. 1990; Tabuse et al. 1998; Hung and Kemphues 1999; Hao et al. 2006). PAR-3 and PAR-6 are conserved proteins that contain PDZ-domains and form a complex with PKC-3 (Etemad-Moghadam et al. 1995; Izumi et al. 1998; Tabuse et al. 1998; Hung and Kemphues 1999). This complex becomes restricted to the anterior cortex of the embryo in response to spatially defined actomyosin contractions occurring in the embryo upon fertilization (Goldstein and Hird 1996; Munro et al. 2004). The posterior cortex of the embryo that becomes devoid of the anterior PAR proteins is occupied by the RING protein PAR-2 and the Ser/Thr kinase PAR-1 (Guo and Kemphues 1995; Boyd et al. 1996; Cuenca et al. 2003). Once polarized, the anterior and posterior PAR proteins mutually exclude each other from their respective cortices (Etemad-Moghadam et al. 1995; Boyd et al. 1996; Cuenca et al. 2003; Hao et al. 2006). Loss of function of the gene par-1, as opposed to loss of most other par genes, results in embryos that exhibit only subtle effects on the polarized cortical domains occupied by the other PAR proteins (Cuenca et al. 2003). However defects in this gene are associated with a more symmetric division in size, an aberrant distribution of cell fate specification markers, altered cell fates of the daughter cells of the embryo, and ultimately embryonic lethality (Kemphues et al. 1988; Guo and Kemphues 1995).PAR-1 controls asymmetric cell division and cell fate specification by regulating the localization of the two highly similar CCCH-type zinc-finger proteins MEX-5 and MEX-6 (referred to as MEX-5/6). MEX-5 and MEX-6 are 70% identical in their amino acid sequence and fulfill partially redundant functions in the embryo (Schubert et al. 2000). In wild-type animals, endogenous MEX-5 and GFP fusions of MEX-6 localize primarily to the anterior of the embryo while both proteins are evenly distributed in par-1 mutant embryos (Schubert et al. 2000; Cuenca et al. 2003). This suggests that in wild-type animals, PAR-1 acts in part by restricting MEX-5 and MEX-6 to the anterior of the embryo. The precise mechanism of this regulation is not known, but an elegant study performed for MEX-5 indicates that differential protein mobility in the anterior and posterior cytoplasm of the one-cell embryo contributes to this asymmetry (Tenlen et al. 2008). While increased mobility in the posterior of the one-cell embryo correlates with a par-1- and par-4-dependent phosphorylation on MEX-5, the kinase directly phosphorylating MEX-5 remains to be identified (Tenlen et al. 2008).Some of the phenotypes associated with loss of par-1 function are dependent on the function of mex-5 and mex-6. First, loss of function of par-1 leads to a decreased stability and aberrant localization of the posterior cell fate specification marker PIE-1, a protein that is usually inherited by the posterior daughter cell in wild-type animals and ensures the correct specification of the germline (Mello et al. 1996; Seydoux et al. 1996). This decreased stability is dependent on mex-5/6 function as PIE-1 levels are restored, albeit with symmetrical distribution, in mex-6(RNAi); mex-5(RNAi); par-1(b274) embryos (Schubert et al. 2000; Cuenca et al. 2003; Derenzo et al. 2003). Second, embryos lacking par-1 function exhibit decreased amounts of P granules in the one-cell embryo, while these markers are present in mex-6(pk440); mex-5(zu199); par-1(RNAi) embryos of comparable age (Cheeks et al. 2004). Third, in par-1(RNAi) one-cell embryos the posterior cortical domain occupied by the polarity protein PAR-2 is extended anteriorly, when compared to wild-type embryos (Cuenca et al. 2003). This anterior extension is rescued in embryos deficient for both par-1 and mex-5/6 (Cuenca et al. 2003). Taken together, these results indicate that par-1 acts in the embryo—at least in part—by regulating the localization and/or activity of the proteins MEX-5 and MEX-6. However, it remains unclear whether other proteins can modulate PAR-1 function to affect MEX-5/6 activity.To gain insight into the mechanisms of par-1 function in the embryo, we sought to identify genes that act together with par-1 during embryonic development. We performed an RNAi-based screen for genetic interactors of the temperature-sensitive allele par-1(zu310), using the embryonic lethal phenotype of this mutant as a readout. This method has proven successful in previous screens to identify genes involved in early embryonic processes (Labbé et al. 2006; O''Rourke et al. 2007). We were able to identify six genes that, upon disruption of their function, suppress the embryonic lethal phenotype of par-1 mutant embryos. One of these genes is mpk-1, the C. elegans homolog of the highly conserved MAP kinase ERK. Closer analysis subsequently showed that reduction of function of mpk-1 not only increases viability of par-1 mutant embryos, but also reverts several polarity phenotypes associated with loss of function of par-1. Our data indicate that mpk-1 antagonizes par-1 activity to regulate polarization and asymmetric cell divisions in the early embryo.     

Interspecies comparison of a gene pair with partially redundant function: the rst and kirre genes in D. virilis and D. melanogaster

Abstract The D. melanogaster rst and kirre genes encode two highly related immunoglobulin-like cell adhesion molecules that function redundantly during embryonic muscle development. The two genes appear to be derived from a common ancestor by gene duplication. Gene duplications have been proposed to be of major evolutionary significance since duplicated redundant sequences can accumulate mutations without detrimental effects for the organism and leave the duplicated genes free to assume novel functions. To address the issue of conservation of the duplicated sequences and their putative redundancy, as well as to identify putative functional divergence of the paralogs during drosophilid evolution, we performed an interspecies comparison of the rst and kirre genes from D. virilis and D. melanogaster. The D. virilis genome contains orthologues of both rst and kirre and hence the duplication took place before the split of the two lineages and has subsequently been conserved. However, whilst the Rst orthologues show a high degree of sequence similarity, this similarity is lower in Kirre orthologues. Especially the intracellular domains of D. virilis and D. melanogaster Kirre sequences are highly divergent: the D. virilis kirre gene lacks the 3′-most exon present in D. melanogaster, which contains motifs conserved between kirre and rst in D. melanogaster. Hence, while each of the two genes is highly conserved at the level of its exon-intron organization, the selection forces acting on the rst and kirre coding sequences are different. These findings are discussed in the light of general evolutionary mechanisms.     

Delimiting regulatory sequences of the Drosophila melanogaster Ddc gene.

We delimited sequences necessary for in vivo expression of the Drosophila melanogaster dopa decarboxylase gene Ddc. The expression of in vitro-altered genes was assayed following germ line integration via P-element vectors. Sequences between -209 and -24 were necessary for normally regulated expression, although genes lacking these sequences could be expressed at 10 to 50% of wild-type levels at specific developmental times. These genes showed components of normal developmental expression, which suggests that they retain some regulatory elements. All Ddc genes lacking the normal immediate 5'-flanking sequences were grossly deficient in larval central nervous system expression. Thus, this upstream region must contain at least one element necessary for this expression. A mutated Ddc gene without a normal TATA boxlike sequence used the normal RNA start points, indicating that this sequences is not required for start point specificity.     

Coding and potential regulatory sequences of a cluster of chorion genes in Drosophila melanogaster

We have characterized at the nucleotide level a 4.8-kilobase pair segment of the third chromosome of Droophila melanogaster, which contains a cluster of three chorion genes, s 18-1, s 15-1 and s 19-1. These genes are tandemly oriented and share the same basic organization: a small and a large exon separated by a short intron in the signal peptide region. In the coding region, limited similarities at the DNA and protein level suggest a common but distant evolutionary origin. The flanking sequences were searched for elements that might be involved in controlling the tissue-specific and temporally regulated expression and the selective amplification of the chorion genes. A good candidate for a cis-regulatory element is the hexamer, TCACGT, which is found in all three genes in a highly significant position, 23 to 27 nucleotides upstream of the TATA-box, accompanied by additional, less exact similarities. Palindromes and short inverted repeats that are found in the vicinity of their complement are non-uniformly distributed: they are most concentrated in the 3 flanking part of all three genes, in and near regions of unusually high A and T content. The highest number of dyad symmetries, remiiscent of sequences that function as viral replication origins, is found associated with the T- and A-rich regions between genes s18-1 and s15-1.     

Extra sequences found at P element excision sites in Drosophila melanogaster.

Summary. We have previously established a transgenic Drosophila line with a highly transposable P element insertion. Using this strain we analyzed transposition and excision of the P element at the molecular level. We examined sequences flanking the new insertion sites and those of the remnants after excision. Our results on mobilization of the P element demonstrate that target-site duplication at the original insertion site does not play a role in forward excision and transposition. After P element excision an 8 by target-site duplication and part of the 31 by terminal inverted repeat (5–18 bp) remained in all the strains examined. Moreover, in 11 out of 28 strains, extra sequences were found between the two remaining inverted repeats. The double-strand gap repair model does not explain the origin of these extra sequences. The mechanism creating them may be similar to the hairpin model proposed for the transposon Tam in Antirrhinum majus.     

Nucleotide variation at the runt locus in Drosophila melanogaster and Drosophila simulans.

Intra- and interspecific nucleotide variation for the major developmental gene runt in Drosophila was studied in D. melanogaster and D. simulans. The 1.5-kb protein-coding region and the 0.4-kb intron of the runt gene were sequenced for 11 alleles in each species. The D. melanogaster alleles originated from east Africa. Estimated parameters of intraspecific variation in D. melanogaster (exons: theta = 0.018, pi = 0.018; intron: theta = 0.014, pi = 0.014) and D. simulans (exons: theta = 0.007, pi = 0.005; intron: theta = 0.008, pi = 0.005) were below average for other X-linked genes, while divergence between species (exons: D = 0.094; intron: D = 0.069) fell within the normal range for both silent and replacement changes. This estimate for runt, along with published values for three other genes in regions of normal recombination, show east African D. melanogaster to be roughly twice as polymorphic as D. simulans. The majority of nucleotide variation, silent and replacement, in both species was found to be selectively neutral using various statistical tests (HKA, McDonald-Kreitman, Tajima, and Fu and Li tests). Monte Carlo simulations of the coalescent process significantly rejected a Wright-Fisher model with respect to an amino acid polymorphism and the distribution of polymorphic sites among the D. simulans lines. This indicated an old lineage and may reflect ancestral population substructuring in D. simulans.     

Tracing the evolution of H-2 D region genes using sequences associated with a repetitive element

The class I genes in the murine MHC are genetically divided into the K, D, Qa, and T1a region subfamilies. These genes presumably arose by duplication from a common class I ancestor. Oligonucleotide probes specific for sequences associated with a moderately repetitive B2 SINE element, which is inserted into the 3' untranslated region of the H-2D and H-2L genes, were used to examine the evolutionary relationship between these classically defined D region genes (H-2D and H-2L) and the other members of the class I gene family. Hybridization analyses of recombinant cosmid and genomic DNA indicated that the D region genes separated genetically from the other members of the class I gene family 12 to 14 million years ago. The evidence suggests that during this time frame the chromosomal segment harboring the characteristic insertion became fixed in the ancestral population which gave rise to Mus domesticus. Previous studies have shown that the number of genes present in the Qa and T1a regions varies among inbred strains and among laboratory stocks of wild mice derived from more distant species on the genus Mus. No evidence was found in this study to support the hypothesis that variation in class I gene number is the result of recent duplications of the functionally defined class I genes of the D region, H-2D and H-2L.     

Divergence of U2 snRNA sequences in the genome of D. melanogaster.

Four different U2-snRNA genes/related sequences of D. melanogaster were cloned and characterized. The sequences of all four genes suggest that they were generated by a DNA-mediated mechanism. These genes/related sequences were found to be located in two loci, each locus containing two U2 snRNA sequences. Using coding sequences as well as flanking sequences as hybridization probes against polytene chromosomes of D. melanogaster Oregon R we were able to map these loci separately at positions 34BC and 84C. By Northern analysis we observed that the quantities of U2- and U1-snRNA are coordinated and change during the embryonic development of the fly.