Genomic heterogeneity of background substitutional patterns in Drosophila melanogaster

Mutation is the underlying force that provides the variation upon which evolutionary forces can act. It is important to understand how mutation rates vary within genomes and how the probabilities of fixation of new mutations vary as well. If substitutional processes across the genome are heterogeneous, then examining patterns of coding sequence evolution without taking these underlying variations into account may be misleading. Here we present the first rigorous test of substitution rate heterogeneity in the Drosophila melanogaster genome using almost 1500 nonfunctional fragments of the transposable element DNAREP1_DM. Not only do our analyses suggest that substitutional patterns in heterochromatic and euchromatic sequences are different, but also they provide support in favor of a recombination-associated substitutional bias toward G and C in this species. The magnitude of this bias is entirely sufficient to explain recombination-associated patterns of codon usage on the autosomes of the D. melanogaster genome. We also document a bias toward lower GC content in the pattern of small insertions and deletions (indels). In addition, the GC content of noncoding DNA in Drosophila is higher than would be predicted on the basis of the pattern of nucleotide substitutions and small indels. However, we argue that the fast turnover of noncoding sequences in Drosophila makes it difficult to assess the importance of the GC biases in nucleotide substitutions and small indels in shaping the base composition of noncoding sequences.     

Molecular evolution of P transposable elements in the genus Drosophila. I. The saltans and willistoni species groups

A phylogenetic survey using the polymerase chain reaction (PCR) has identified four major P element subfamilies in the saltans and willistoni species groups of Drosophila. One subfamily, containing about half of the sequences studied, consists of elements that are very similar to the canonical (and active) P element from D. melanogaster. Within this subfamily, nucleotide sequence differentiation among different copies from the same species and among elements from different species is relatively low. This observation suggests that the canonical elements are relatively recent additions to the genome or, less likely, are evolving slowly relative to the other subfamilies. Elements belonging to the three noncanonical lineages are distinct from the canonical elements and from one another. Furthermore, there is considerably more sequence variation, on the average, within the noncanonical subfamilies compared to the canonical elements. Horizontal transfer and the coexistence of multiple, independently evolving element subfamilies in the same genome may explain the distribution of P elements in the saltans and willistoni species groups. Such explanations are not mutually exclusive, and each may be involved to varying degrees in the maintenance of P elements in natural populations of Drosophila.      

Diversity in isochore structure among cold-blooded vertebrates based on GC content of coding and non-coding sequences

To review the general consideration about the different compositional structure of warm and cold-blooded vertebrates genomes, we used of the increasing number of genetic sequences, including coding (exons) and non-coding (introns) regions, that have been deposited on the databases throughout last years. The nucleotide distributions of the third codon positions (GC3) have been analyzed in 1510 coding sequences (CDS) of fish, 1414 CDS of amphibians and 320 CDS of reptiles. Also, the relationship between GC content of 74, 56 and 25 CDS of fish, amphibians and reptiles, respectively and that of their corresponding introns (GCI) have been considerated. In accordance with recent data, sequence analysis showed the presence of very GC3-rich CDS in these poikilotherm vertebrates. However, very high diversity in compositional patterns among different orders of fish, amphibians and reptiles was found. Significant positive correlations between GC3 and GCI was also confirmed for the genes analyzed. Nevertheless, introns resulted to be poorer in GC than their corresponding CDS, this difference being larger than in human genome. Because the limited number of available sequences including exons and introns we must be cautious about the results derived from them. However, the indicious of higher GC richness of coding sequences than of their corresponding introns could aid to understand the discrepancy of sequence analysis with the ultracentrifugation studies in cold-blooded vertebrates that did not predict the existence of GC-rich isochores.     

Phylogeny of the Drosophila saltans species group based on combined analysis of nuclear and mitochondrial DNA sequences

Nucleotide sequences from two nuclear loci, alcohol dehydrogenase and internal transcribed spacer-1 of the nuclear ribosomal DNA repeats, and two mitochondrial genes, cytochrome oxidase I and cytochrome oxidase II, were determined from nine species in the Drosophila saltans species group. The partition homogeneity test and partitioned Bremer support were used to measure incongruence between phylogenetic hypotheses generated from individual partitions. Individual loci were generally congruent with each other and consistent with the previously proposed morphological hypothesis, although they differed in level of resolution. Since extreme conflict between partitions did not exist, the data were combined and analyzed simultaneously. The total evidence method gave a more resolved and highly supported phylogeny, as indicated by bootstrap proportions and decay indices, than did any of the individual analyses. The cordata and elliptica subgroups, considered to have diverged early in the history of the D. saltans group, were sister taxa to the remainder of the saltans group. The sturtevanti subgroup, represented by D. milleri and D. sturtevanti, occupies an intermediate position in this phylogeny. The saltans and parasaltans subgroups are sister clades and occupy the most recently derived portion of the phylogeny. As with previous morphological studies, phylogenetic relationships within the saltans subgroup were not satisfactorily resolved by the molecular data.      

The Drosophila mobile element jockey is a LINE element and contains coding sequences homologous to retroviral proteins

A detailed investigation of Drosophila melanogaster mobile dispersed repetitive element jockey is performed. Its structural features resemble those of LINE elements. Sequencing of the complete jockey 5020 bp in length revealed two long open reading frames ORF1 and ORF2 overlapping with a frameshift-1. Judging by amino acid homologies, ORF1 encodes a nucleic acid binding protein, characteristic of replication competent retroviruses; the 3' part of ORF2 encodes an RNA-dependent DNA polymerase which has an amino acid sequence, similar to recently published sequences of LINE elements of Drosophila, Trypanosoma and mammals. This fact demonstrates their evolutionary relationship. Sequencing of several deleted copies of jockey revealed the absence of the major part of ORF2, though the rest of the element, including the ends, is highly conservative.     

Kin grouping is insufficient to explain the inclusive fitness gains of conspecific brood parasitism in the common eider

Conspecific brood parasitism allows females to exploit other females' nests and enhance their reproductive output. Here, we test a recent theoretical model of how host females gain inclusive fitness from brood parasitism. High levels of relatedness between host and parasitizer can be maintained either by: (a) kin recognizing and parasitizing each other as a form of cooperative breeding or (b) natal philopatry and nest site fidelity facilitating the formation of kin groups, thereby increasing the probability of parasitism between relatives nesting in close proximity. To address these two hypotheses we genotyped feathers and hatch membranes of common eiders (Somateria mollissima) from western Hudson Bay, Canada, using a noninvasive sampling methodology. We found that most instances of brood parasitism do result in inclusive fitness gains. Furthermore, females with failed nests moved an average of 492 m from their previous year's nest site, while successful females only moved an average of 13 m. Therefore, we observed host–parasite relatedness can occur at levels higher than would be expected by chance even in the absence of kin grouping, suggesting that closely related females nesting near one another is not essential to maintain high host–parasitizer relatedness. In addition, kin grouping is only a transient phenomenon that cannot occur every year due to the propensity for females of failed nests to nest farther away from their nest site in subsequent years than females with successful nests, which provides support for kin recognition as a more likely mechanism to maintain high host–parasitizer relatedness over time.     

Divergence of the Regulation of alpha-Amylase Activity in Drosophila melanogaster, Drosophila funebris, and Drosophila saltans

The regulation of amylase activity in threeDrosophila species, D. melanogaster,D. funebris and D. saltans, wasanalyzed by measuring the specific activity levels infour dietary environments, cornmeal, glucose, 5% starch, and 10% starch, at threedevelopmental stages, i.e., the third-instar larval,pupal, and 2-day-old adult stages. The developmentalprofiles of amylase activity for the threeDrosophila species showed that the level of activity washigh at the larval and adult stages but substantiallylow at the pupal stage, suggesting thatDrosophila does not utilize starch at the pupalstage. Divergence in the regulation of amylase was observed amongthe three Drosophila species on the followingpoints. (1) The order of amylase specific activity wasD. melanogaster > D. funebris >D. saltans. (2) The response pattern to the dietary environment varied amongthe species and changed during development. (3) Thetiming of the switch in the response pattern to thedietary environment during development was before pupation in D. funebris and D.saltans but after pupation in D.melanogaster. The significance of the divergence inthe regulation of amylase activity for adaptation to astarch environment in Drosophila is discussed.     

Evolution of the Adh locus in the Drosophila willistoni group: the loss of an intron, and shift in codon usage

We report here the DNA sequence of the alcohol dehydrogenase gene (Adh) cloned from Drosophila willistoni. The three major findings are as follows: (1) Relative to all other Adh genes known from Drosophila, D. willistoni Adh has the last intron precisely deleted; PCR directly from total genomic DNA indicates that the deletion exists in all members of the willistoni group but not in any other group, including the closely related saltans group. Otherwise the structure and predicted protein are very similar to those of other species. (2) There is a significant shift in codon usage, especially compared with that in D. melanogaster Adh. The most striking shift is from C to U in the wobble position (both third and first position). Unlike the codon-usage-bias pattern typical of highly biased genes in D. melanogaster, including Adh, D. willistoni has nearly 50% G + C in the third position. (3) The phylogenetic information provided by this new sequence is in agreement with almost all other molecular and morphological data, in placing the obscura group closer to the melanogaster group, with the willistoni group farther distant but still clearly within the subgenus Sophophora.      

Structural and sequence diversity of the transposon Galileo in the Drosophila willistoni genome

Background

Galileo is one of three members of the P superfamily of DNA transposons. It was originally discovered in Drosophila buzzatii, in which three segregating chromosomal inversions were shown to have been generated by ectopic recombination between Galileo copies. Subsequently, Galileo was identified in six of 12 sequenced Drosophila genomes, indicating its widespread distribution within this genus. Galileo is strikingly abundant in Drosophila willistoni, a neotropical species that is highly polymorphic for chromosomal inversions, suggesting a role for this transposon in the evolution of its genome.

Results

We carried out a detailed characterization of all Galileo copies present in the D. willistoni genome. A total of 191 copies, including 133 with two terminal inverted repeats (TIRs), were classified according to structure in six groups. The TIRs exhibited remarkable variation in their length and structure compared to the most complete copy. Three copies showed extended TIRs due to internal tandem repeats, the insertion of other transposable elements (TEs), or the incorporation of non-TIR sequences into the TIRs. Phylogenetic analyses of the transposase (TPase)-encoding and TIR segments yielded two divergent clades, which we termed Galileo subfamilies V and W. Target-site duplications (TSDs) in D. willistoni Galileo copies were 7- or 8-bp in length, with the consensus sequence GTATTAC. Analysis of the region around the TSDs revealed a target site motif (TSM) with a 15-bp palindrome that may give rise to a stem-loop secondary structure.

Conclusions

There is a remarkable abundance and diversity of Galileo copies in the D. willistoni genome, although no functional copies were found. The TIRs in particular have a dynamic structure and extend in different ways, but their ends (required for transposition) are more conserved than the rest of the element. The D. willistoni genome harbors two Galileo subfamilies (V and W) that diverged ~9 million years ago and may have descended from an ancestral element in the genome. Galileo shows a significant insertion preference for a 15-bp palindromic TSM.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-792) contains supplementary material, which is available to authorized users.     

Neurogenesis of the peripheral nervous system in Drosophila embryos: DNA replication patterns and cell lineages

Cell lineages that give rise to the PNS were studied using the thymidine analog 5-bromo-2'-deoxyuridine (BrdU) to visualized DNA replication immunocytochemically. The precursors of the PNS in the body segments of Drosophila embryos replicate their DNA in a spatially and temporally stereotyped pattern. The sequence of DNA replication within developing sensory organs suggests particular lineage relationships of the cells that constitute a sensory organ, i.e., neuron and associated support cells. In embryos that are mutant for the achaete-scute complex or daughterless, in which most or all of the PNS is missing, no BrdU-labeled cells were found in the appropriate regions, suggesting that these PNS precursors either do not form or fail to replicate. Thus, the BrdU technique allows determination as to whether a mutation affects the PNS precursors or terminal differentiation.     

Evolution of the GC content of the histone 3 gene in seven Drosophila species

The molecular evolution of the histone multigene family was studied by cloning and determining the nucleotide sequences of the histone 3 genes in seven Drosophila species, D. takahashii, D. lutescens, D. ficusphila, D. persimilis, D.pseudoobscura, D. americana and D. immigrans. CT repeats, a TATA box and an AGTG motif in the 5' region, and a hairpin loop and purine-rich motifs (CAA(T/G)GAGA) in the 3' region were conserved even in distantly related species. In D. hydei and D.americana, the GC content at the third codon position in the protein coding region was relatively low (49% and 45%), while in D. takahashii and D. lutescens it was relatively high (64% and 65%). The non- significant correlation between the GC contents in the 3' region and at the third codon position as well as the evidence of less constraint in the 3' region suggested that mutational bias may not be the major mechanism responsible for the biased nucleotide change at the third codon position or for codon usage bias.     

Interspecific and intraspecific comparisons of the period locus in the Drosophila willistoni sibling species

The period (per) locus has received much attention in molecular evolution studies because it is one of the best studied "behavioral genes" and because it offers insight into the evolution of repetitive sequences. We studied most of the coding region of per in Drosophila willistoni and confirmed previously observed patterns of conservation and divergence among distantly related species. Five regions are so highly diverged that they cannot be aligned, whereas a region encompassing the PAS domain is very conserved. Structural and nucleotide polymorphism patterns in the willistoni group are not the same as those observed in previously studied species. We sequenced the region homologous to the highly polymorphic threonine-glycine repeat of D. melanogaster in multiple strains of D. willistoni, as well as in other members of willistoni group, and found an unusual amount of conservation in this region. However, the next nonconserved region downstream in the sequence is quite variable and polymorphic for the number of repeated glycines. The glycine codon usage is significantly different in this glycine repeat as compared to other parts of the gene. We were able to plot the directionality of change in the glycine repeat region onto a phylogeny and find that the addition of glycines is the general trend with the diversification of the willistoni group.      

Localization of sequences coding for histone messenger RNA in the chromosomes of Drosophila melanogaster

In situ hybridization of sea urchin (Psammechinus miliaris, Lytechinus pictus and Strongylocentrotus purpuratus) histone messenger RNA has been used to map complementary sequences on polytene chromosomes from Drosophila melanogaster. The sea urchin RNA hybridizes to the polytene regions from 39D3 through 39E1-2, including both of these bands (39D2 may also be included). This region is identical to the one which hybridizes most heavily with non-polyadenylated cytoplasmic RNA from D. melanogaster tissues. Sea urchin mRNAs coding for several individual histones each hybridize across the entire region from 39D3 (or D2) through 39E1-2, as would be expected if the individual mRNA sequences are interspersed. In view of the apparently even distribution of sequences complementary to histone mRNA within the 39D3-39E1-2 region, the significance of the several polytene bands in this region remains an open question. Biochemical characterization of the hybrids between sea urchin histone mRNA and D. melanogaster DNA suggests that sea urchin mRNAs for several of the histone classes have some portions which retain enough sequence homology with the D. melanogaster sequences to form hybrids, although the hybrids have base pair mismatches. In situ hybridization of chromosomes in which region 39D-E is ectopically paired show no evidence of sequence homology in the chromosome region with which 39D-E is associated.     

Mitochondrial mutations and aging: random drift is insufficient to explain the accumulation of mitochondrial deletion mutants in short‐lived animals

Mitochondrial DNA deletions accumulate over the life course in post‐mitotic cells of many species and may contribute to aging. Often a single mutant expands clonally and finally replaces the wild‐type population of a whole cell. One proposal to explain the driving force behind this accumulation states that random drift alone, without any selection advantage, is sufficient to explain the clonal accumulation of a single mutant. Existing mathematical models show that such a process might indeed work for humans. However, to be a general explanation for the clonal accumulation of mtDNA mutants, it is important to know whether random drift could also explain the accumulation process in short‐lived species like rodents. To clarify this issue, we modelled this process mathematically and performed extensive computer simulations to study how different mutation rates affect accumulation time and the resulting degree of heteroplasmy. We show that random drift works for lifespans of around 100 years, but for short‐lived animals, the resulting degree of heteroplasmy is incompatible with experimental observations.     

Patterns of intron sequence evolution in Drosophila are dependent upon length and GC content

Background

Introns comprise a large fraction of eukaryotic genomes, yet little is known about their functional significance. Regulatory elements have been mapped to some introns, though these are believed to account for only a small fraction of genome wide intronic DNA. No consistent patterns have emerged from studies that have investigated general levels of evolutionary constraint in introns.

Results

We examine the relationship between intron length and levels of evolutionary constraint by analyzing inter-specific divergence at 225 intron fragments in Drosophila melanogaster and Drosophila simulans, sampled from a broad distribution of intron lengths. We document a strongly negative correlation between intron length and divergence. Interestingly, we also find that divergence in introns is negatively correlated with GC content. This relationship does not account for the correlation between intron length and divergence, however, and may simply reflect local variation in mutational rates or biases.

Conclusion

Short introns make up only a small fraction of total intronic DNA in the genome. Our finding that long introns evolve more slowly than average implies that, while the majority of introns in the Drosophila genome may experience little or no selective constraint, most intronic DNA in the genome is likely to be evolving under considerable constraint. Our results suggest that functional elements may be ubiquitous within longer introns and that these introns may have a more general role in regulating gene expression than previously appreciated. Our finding that GC content and divergence are negatively correlated in introns has important implications for the interpretation of the correlation between divergence and levels of codon bias observed in Drosophila.     

Further Study on the Esterase Patterns of Sibling Species in the Drosophila saltans Subgroup (saltans Group): Intraspecific and Interspecific Variations in the Development

Twenty of the 32 esterase bands previously detected in the adults of D. prosaltans, D. saltans and D.␣austrosaltans were found in larvae and pupae studied in this work. The results showed that, in addition to expressing the highest number of esterase bands, the adult stage of the three species exhibited the highest degree of expression (amount of synthesis) for most of the bands. Differences between larval and pupal stages were detected in the degree of expression (amount of synthesis) of the bands and in the frequency of samples expressing them. The frequencies of expression of the bands corresponding to genes in loci 1–3 were greater in pupae than in larvae while the frequencies of expression of the bands corresponding to genes in loci 4–9 were predominantly expressed in larvae or were equal in both developmental stages. Like the adults, larvae, pupae and empty pupal cases (which were also studied in this work) showed specific esterases. Taken together, the observations showed that, in the species studied, every developmental stage is characterized by specific bands and by specific frequency and degree of expression of the bands shared with other stages.     

Different patterns of evolution in alternative and constitutive coding regions of Drosophila alternatively spliced genes

Seven hundred and ninety Drosophila melanogaster genes, alternatively spliced in coding regions were considered together with their Drosophila pseudoobscura orthologs. It was found that nucleotide substitutions in alternative coding regions accumulate more intensively than in constitutive regions. Moreover, the evolutionary pattern of alternative regions depends on their inclusion mechanisms (use of alternative promoters, splicing sites or polyadenylation sites) significantly. The rate of synonymous substitutions varies is more dramatically than that of nonsynonymous substitutions. Nucleotide substitution patterns in different classes of alternative regions of mammalian and Drosophila genes have little in common.