Polymorphism of canonical and noncanonical <Emphasis Type="Italic">gypsy</Emphasis> sequences in different species of <Emphasis Type="Italic">Drosophila </Emphasis><Emphasis Type="Italic">melanogaster</Emphasis> subgroup: possible evolutionary relations

Mobile genetic elements constitute a substantial part of eukaryotic genome and play an important role in its organization and functioning. Co-evolution of retrotransposons and their hosts resulted in the establishment of control systems employing mechanisms of RNA interference that seem to be impossible to evade. However, “active” copies of endogenous retrovirus gypsy escape cellular control in some cases, while its evolutionary elder “inactive” variants do not. To clarify the evolutionary relationship between “active” and “inactive” gypsy we combined two approaches: the analysis of gypsy sequences, isolated from G32 Drosophila melanogaster strain and from different Drosophila species of the melanogaster subgroup, as well as the study of databases, available on the Internet. No signs of “intermediate” (between “active” and “inactive”) gypsy form were found in GenBank, and four full-size G32 gypsy copies demonstrated a convergence that presumably involves gene conversion. No “active” gypsy were revealed among PCR generated gypsy ORF3 sequences from the various Drosophila species indicating that “active” gypsy appeared in some population of D. melanogaster and then started to spread out. Analysis of sequences flanking gypsy variants in G32 revealed their predominantly heterochromatic location. Discrepancy between the structure of actual gypsy sites in G32 and corresponding sequences in database might indicate significant inter-strain heterochromatin diversity. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Retrotransposon gtwin: structural analysis and distribution in Drosophila melanogaster strains

A search for noncanonical variants of the gypsy retrotransposon (MDG4) in the genome of the Drosophila melanogaster strain G32 led to the cloning of four copies of the poorly studied 7411-bp gtwin element. Sequence analysis showed that gtwin belongs to a family of endogeneous retroviruses, which are widespread in the Drosophila genome and have recently been termed insect erantiviruses. The gtwin retrotransposon is evolutionarily closest to MDG4, as evident from a good alignment of their nucleotide sequences including ORF1 (the pol gene) and ORF3 (the env gene), as well as the amino acid sequences of their protein products. These regions showed more than 75% homology. The distribution of gtwin was studied in several strains of the genus Drosophila. While strain G32 contained more than 20 copies of the element, ten other D. melanogaster strains carried gtwin in two to six copies per genome. The gtwin element was not detected in D. hydei or D. virilis. Comparison of the cloned gtwin sequences with the gtwin sequence available from the D. melanogaster genome database showed that the two variants of the mobile element differ by the presence or absence of a stop codon in the central region of ORF3. Its absence from the gtwin copies cloned from the strain G32 may indicate an association between the functional state of ORF3 and amplification of the element.     

<Emphasis Type="Italic">HB</Emphasis> mobile element in the <Emphasis Type="Italic">Drosophila melanogaster</Emphasis> genome: Structural and functional analyses

The structure was analyzed for 60 annotated copies of the mobile genetic element (MGE) HB from the Drosophila melanogaster genome. The genomic distribution of HB copies was studied, and preferential insertion sites (hot spots) were identified, which presumably amount to several kilobases. Structural analysis of the open reading frame (ORF) and terminal repeats of HB was performed. All 26 HB copies retaining the ORF sequence have a stop codon in the same position. Consequently, the HB ORF proved indeed to code for an enzyme of 148 amino acid residues, relatively small for Tc1-family transposases. The ORF consensus sequence was established. HB{}1185 was identified as the only HB copy potentially coding for a functional protein. All 37 repeat-containing HB copies were analyzed. Of these, only four had functional terminal sequences, lacking, however, a functional transposase gene. A new 7762-bp copy of MGE roo was found in the D. melanogaster genome; the copy was earlier unavailable from databases and represents an insert in the HB{}1605 sequence.     

Evolutionary pattern of the gtwin retrotransposon in the Drosophila melanogaster subgroup

The gtwin retrotransposon was recently discovered in the Drosophila melanogaster genome and it is evolutionarily closer to gypsy endogenous retrovirus. This study has identified gtwin homologous sequences in the genome of D. simulans, D. sechellia, D. erecta and D. yakuba by performing homology searches against the public genome database of Drosophila species. The phylogenetic analyses of the gtwin env gene sequences of these species have shown some incongruities with the host species phylogeny, suggesting some horizontal transfer events for this retroelement. Moreover, we reported the existence of DNA sequences putatively encoding full-length Env proteins in the genomes of Drosophila species other than D. melanogaster. The results suggest that the gtwin element may be an infectious retrovirus able to invade the genome of new species, supporting the gtwin evolutionary picture shown in this work.     

Duplicate Gene Evolution Toward Multiple Fates at the <Emphasis Type="Italic">Drosophila melanogaster HIP</Emphasis>/<Emphasis Type="Italic">HIP-Replacement</Emphasis> Locus

Hsc/Hsp70-interacting protein (HIP) is a rapidly evolving Hsp70 cofactor. Analyses of multiple Drosophila species indicate that the HIP gene is duplicated only in D. melanogaster. The HIP region, in fact, contains seven distinctly evolving duplicated genes. The regional duplication occurred in two steps, fixed rapidly, and illustrates multiple modes of duplicate gene evolution. HIP and its duplicate HIP-R are adaptively evolving in a manner unique to the region: they exhibit elevated divergence from other drosophilids and low polymorphism within D. melanogaster. HIP and HIP-R are virtually identical, share polymorphisms, and are subject to gene conversion. In contrast, two other duplicate genes in the region, CG33221 and GP-CG32779, are pseudogenes, and the chimeric gene Crg1 is subject to balancing selection. HIP and HIP-R are evolving rapidly and adaptively; however, positive selection is not sufficient to explain the molecular evolution of the region as a whole. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Evolutionary History and Mode of the <Emphasis Type="Italic">amylase</Emphasis> Multigene Family in <Emphasis Type="Italic">Drosophila</Emphasis>

Previous studies indicate that the tandemly repeated members of the amylase (Amy) gene family evolved in a concerted manner in the melanogaster subgroup and in some other species. In this paper, we analyzed all of the 49 active and complete Amy gene sequences in Drosophila, mostly from subgenus Sophophora. Phylogenetic analysis indicated that the two types of diverged Amy genes in the Drosophila montium subgroup and Drosophila ananassae, which are located in distant chromosomal regions from each other, originated independently in different evolutionary lineages of the melanogaster group after the split of the obscura and melanogaster groups. One of the two clusters was lost after duplication in the melanogaster subgroup. Given the time, 24.9 mya, of divergence between the obscura and the melanogaster groups (Russo et al. 1995), the two duplication events were estimated to occur at about 13.96 ± 1.93 and 12.38 ± 1.76 mya in the montium subgroup and D. ananassae, respectively. An accelerated rate of amino acid changes was not observed in either lineage after these gene duplications. However, the G+C contents at the third codon positions (GC3) decreased significantly along one of the two Amy clusters both in the montium subgroup and in D. ananassae right after gene duplication. Furthermore, one of the two types of the Amy genes with a lower GC3 content has lost a specific regulatory element within the montium subgroup species and D. ananassae. While the tandemly repeated members evolved in a concerted manner, the two types of diverged Amy genes in Drosophila experienced frequent gene duplication, gene loss, and divergent evolution following the model of a birth-and-death process.     

The <Emphasis Type="Italic">Caenorhabditis</Emphasis><Emphasis Type="Italic">briggsae</Emphasis> genome contains active <Emphasis Type="Italic">CbmaT1</Emphasis> and <Emphasis Type="Italic">Tcb1</Emphasis> transposons

The maT clade of transposons is a group of transposable elements intermediate in sequence and predicted protein structure to mariner and Tc transposons, with a distribution thus far limited to a few invertebrate species. We present evidence, based on searches of publicly available databases, that the nematode Caenorhabditis briggsae has several maT-like transposons, which we have designated as CbmaT elements, dispersed throughout its genome. We also describe two additional transposon sequences that probably share their evolutionary history with the CbmaT transposons. One resembles a fold back variant of a CbmaT element, with long (380-bp) inverted terminal repeats (ITRs) that show a high degree (71%) of identity to CbmaT1. The other, which shares only the 26-bp ITR sequences with one of the CbmaT variants, is present in eight nearly identical copies, but does not have a transposase gene and may therefore be cross mobilised by a CbmaT transposase. Using PCR-based mobility assays, we show that CbmaT1 transposons are capable of excising from the C. briggsae genome. CbmaT1 excised approximately 500 times less frequently than Tcb1 in the reference strain AF16, but both CbmaT1 and Tcb1 excised at extremely high frequencies in the HK105 strain. The HK105 strain also exhibited a high frequency of spontaneous induction of unc-22 mutants, suggesting that it may be a mutator strain of C. briggsae.     

Analysis of the chromocenter DNA composition in polytene chromosomes of <Emphasis Type="Italic">Drosophila orena</Emphasis> ovarian nurse cells

Chromocenter DNA fragments of polytene chromosomes of Drosophila orena ovarian nurse cells were cloned from a region-specific library (Dore1) in a plasmid vector to yield 133 clones. A total of 76 clones were selected and sequenced. The total length of the sequenced fragments was 23940 bp. Analysis with several software packages revealed various repetitive sequences among the fragments of the Dore1 library, including mobile genetic elements (25 fragments homologous to various LTR retrotransposons, five fragments homologous to LINEs, three fragments homologous to Helitrons, one fragment homologous to Polinton, and one fragment homologous to the mini-me non-LTR retrotransposon), four minisatellites, a satellite (SAR_DM), the (TATATG)n simple sequence repeat, and a low-complexity T-rich repeat. Sequences homologous to protein-coding genes were also found in the Dore1 library. Various repetitive DNA sequences and gene homologs were identified as conserved sequences of pericentric heterochromatin of polytene chromosomes of ovarian nurse cells in nine species of the melanogaster species subgroup.     

First evidence of methylation in the genome of <Emphasis Type="Italic">Drosophila willistoni</Emphasis>

DNA methylation has been studied abundantly in vertebrates and recent evidence confirms that this phenomenon could be disseminated among some invertebrates groups, including Drosophila species. In this paper, we used the Methylation-Sensitive Restriction Endonuclease (MSRE) technique and Southern blot with specific probes, to detect methylation in the Drosophila willistoni species. We found differential cleavage patterns between males and females that cannot be explained by Mendelian inheritance, pointing to a DNA methylation phenomenon different from the Drosophila melanogaster one. The sequencing of some of these bands showed that these fragments were formed by different DNA elements, among which rDNA. We also characterized the D. willitoni dDnmt2 sequence, through a Mega Blast search against the D. willistoni Trace Archive Database using the D. melanogaster dDnmt2 nucleotide sequence as query. The complete analysis of D. willistoni dDnmt2 sequence showed that its promoter region is larger, its dDnmt2 nucleotide sequence is 33% divergent from the D. melanogaster one, Inverted Terminal Repeats (ITRs) are absent and only the B isoform of the enzyme is produced. In contrast, ORF2 is more conserved. Comparing the D. willistoni and D. melanogaster dDnmt2 protein sequences, we found higher conservation in motifs from the large domain, responsible for the catalysis of methyl transfer, and great variability in the region that carries out the recognition of specific DNA sequences (TRD). Globally, our results reveal that methylation of the D. willistoni genome could be involved in a singular process of species-specific dosage compensation and that the DNA methylation in the Drosophila genus can have diverse functions. This could be related to the evolutionary history of each species and also to the acquisition time of the dDnmt2 gene.     

Contrasting Patterns of Sequence Evolution at the Functionally Redundant <Emphasis Type="Italic">bric à brac</Emphasis> Paralogs in <Emphasis Type="Italic">Drosophila melanogaster</Emphasis>

Genes with overlapping expression and function may gradually diverge despite retaining some common functions. To test whether such genes show distinct patterns of molecular evolution within species, we examined sequence variation at the bric à brac (bab) locus of Drosophila melanogaster. This locus is composed of two anciently duplicated paralogs, bab1 and bab2, which are involved in patterning the adult abdomen, legs, and ovaries. We have sequenced the 148 kb genomic region spanning the bab1 and bab2 genes from 94 inbred lines of D. melanogaster sampled from a single location. Two non-coding regions, one in each paralog, appear to be under selection. The strongest evidence of directional selection is found in a region of bab2 that has no known functional role. The other region is located in the bab1 paralog and is known to contain a cis-regulatory element that controls sex-specific abdominal pigmentation. The coding region of bab1 appears to be under stronger functional constraint than the bab2 coding sequences. Thus, the two paralogs are evolving under different selective regimes in the same natural population, illuminating the different evolutionary trajectories of partially redundant duplicate genes.     

Identification of <Emphasis Type="Italic">DUF538</Emphasis> cDNA clone from <Emphasis Type="Italic">Celosia cristata</Emphasis> expressed sequences of nonstressed and stressed leaves

DUF538 domain-containing protein family consists of several plant proteins of unknown functions. This protein family has already been discovered by genome annotation tools and cloned as an inducible gene product under various environmental stress conditions. For the first time, we presented a full length DUF538 cDNA (encoding 170 amino acid residues) clone, which was randomly isolated from Celosia cristata leaf cDNA library constructed under normal growth conditions and consistently amplified from leaf cDNA populations prepared from nonstressed and drought-stressed leaves. We predicted that a DUF538 gene product can be a putative candidate for common stress-related protein (regulatory factor) in the plant system. The nucleotide and deduced amino acid sequences of the isolated clone have been submitted to EMBL data bases under accession no. AJ535713.     

Analysis of a <Emphasis Type="Italic">c</Emphasis><Subscript>0</Subscript><Emphasis Type="Italic">t-1</Emphasis> library enables the targeted identification of minisatellite and satellite families in <Emphasis Type="Italic">Beta vulgaris</Emphasis>

Background  

Repetitive DNA is a major fraction of eukaryotic genomes and occurs particularly often in plants. Currently, the sequencing of the sugar beet (Beta vulgaris) genome is under way and knowledge of repetitive DNA sequences is critical for the genome annotation. We generated a c ₀ t-1 library, representing highly to moderately repetitive sequences, for the characterization of the major B. vulgaris repeat families. While highly abundant satellites are well-described, minisatellites are only poorly investigated in plants. Therefore, we focused on the identification and characterization of these tandemly repeated sequences.     

The first complete <Emphasis Type="Italic">Mag</Emphasis> family retrotransposons discovered in <Emphasis Type="Italic">Drosophila</Emphasis>

A retrotransposon of the Mag family was found in the Drosophila simulans genome for the first time. We also identified novel transposable elements representing the Mag family in seven Drosophila species. The high similarity between the 3’ and 5’ long terminal repeats in the found copies of transposable elements indicates that their retrotransposition has occurred relatively recently. Thus, the Mag family of retrotransposons is quite common for the genus Drosophila.     

Characterization of <Emphasis Type="Italic">EamaT1</Emphasis>, a member of <Emphasis Type="Italic">maT</Emphasis> family of transposable elements from the earthworm <Emphasis Type="Italic">Eisenia andrei</Emphasis> (Annelida,Oligochaeta)

The maT family is a unique clade within the Tc1-mariner superfamily, and their distribution is to date known as being limited to invertebrates. A novel transposon named EamaT1 is described from the genome of the earthworm Eisenia andrei. The full sized EamaT1 was obtained by degenerate and inverse PCR-based amplification. Sequence analysis of multiple copies of the EamaT1, which consisted of 0.9 and 1.4 kb elements, showed that the consensual EamaT1 with inverted terminal repeats (ITRs) of 69 bp was 1,422 bp long and flanked by a duplicated TA dinucleotide. The EamaT1 is present in approximately 120–250 copies per diploid genome but undergoes an inactivation process as a result of accumulating multiple mutations and is nonfunctional. The open reading frame (ORF) of the EamaT1 consensus encoding 356 amino acid sequences of transposase contained a DD37D signature and a conserved paired-like DNA binding motif for the transposition mechanism. The result of ITRs comparison confirmed their consensus terminal sequences (5′-CAGGGTG-3′) and AT-rich region on the internal bases for ITRs-transposase interaction.     

Prevalence and Acquisition of the Genes for Zoocin A and Zoocin A Resistance in <Emphasis Type="Italic">Streptococcus equi</Emphasis> subsp. <Emphasis Type="Italic">zooepidemicus</Emphasis>

Zoocin A is a streptococcolytic enzyme produced by Streptococcus equi subsp. zooepidemicus strain 4881. The zoocin A gene (zooA) and the gene specifying resistance to zoocin A (zif) are adjacent on the chromosome and are divergently transcribed. Twenty-four S. equi subsp. zooepidemicus strains were analyzed to determine the genetic difference among three previously characterized as zoocin A producers (strains 4881, 9g, and 9h) and the 21 nonproducers. LT-PCR and Southern hybridization studies revealed that none of the nonproducer strains possessed zooA or zif. RAPD and PFGE showed that the 24 strains were a genetically diverse population with eight RAPD profiles. S. equi subsp. zooepidemicus strains 9g and 9h appeared to be genetically identical to each other but quite different from strain 4881. Sequences derived from 4881 and 9g showed that zooA and zif were integrated into the chromosome adjacent to the gene flaR. A comparison of these sequences with the genome sequences of S. equi subsp. zooepidemicus strains H70 and MGCS10565 and S. equi subsp. equi strain 4047 suggests that flaR flanks a region of genome plasticity in this species. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.