Evolutionary dynamics of the LTR retrotransposons roo and rooA inferred from twelve complete Drosophila genomes

Background  

Roo is the most abundant retrotransposon in the fruit fly Drosophila melanogaster. Its evolutionary origins and dynamics are thus of special interest for understanding the evolutionary history of Drosophila genome organization. We here study the phylogenetic distribution and evolution of roo, and its highly diverged relative rooA in 12 completely sequenced genomes of the genus Drosophila.     

Sex-specific methylation in <Emphasis Type="Italic">Drosophila</Emphasis>: an investigation of the <Emphasis Type="Italic">Sophophora</Emphasis> subgenus

Epigenetic phenomena have been widely characterized in the genomes of vertebrates and DNA methylation is a key mechanism of epigenetic regulation. The DNA methylation systems of invertebrates and vertebrates show several notable differences. However, the evolutionary implications of those differences only recently began to be revealed. Our study investigated the recurrence of sex-specific methylation, as previously described for the species Drosophila willistoni, in other species of the Sophophora subgenus that present close evolutionary relationship. The MSRE and Southern blot techniques were used to analyze rDNA of some species of the willistoni, melanogaster, saltans and obscura groups of Drosophila and the results suggested that differential DNA methylation between sexes only occurs in Drosophila tropicalis and D. insularis, two sibling species of the willistoni subgroup. However, only using the MSRE technique we could detect sex-specific patterns of DNA methylation in all species of willistoni subgroup. These results indicate that DNA methylation may present important differences, even between closely related species, shedding new light on this Neotropical species complex.     

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.     

Characterization of Maltase Clusters in the Genus <Emphasis Type="Italic">Drosophila</Emphasis>

To reveal evolutionary history of maltase gene family in the genus Drosophila, we undertook a bioinformatics study of maltase genes from available genomes of 12 Drosophila species. Molecular evolution of a closely related glycoside hydrolase, the α-amylase, in Drosophila has been extensively studied for a long time. The α-amylases were even used as a model of evolution of multigene families. On the other hand, maltase, i.e., the α-glucosidase, got only scarce attention. In this study, we, therefore, investigated spatial organization of the maltase genes in Drosophila genomes, compared the amino acid sequences of the encoded enzymes and analyzed the intron/exon composition of orthologous genes. We found that the Drosophila maltases are more numerous than previously thought (ten instead of three genes) and are localized in two clusters on two chromosomes (2L and 2R). To elucidate the approximate time line of evolution of the clusters, we estimated the order and dated duplication of all the 10 genes. Both clusters are the result of ancient series of subsequent duplication events, which took place from 352 to 61 million years ago, i.e., well before speciation to extant Drosophila species. Also observed was a remarkable intron/exon composition diversity of particular maltase genes of these clusters, probably a result of independent intron loss after duplication of intron-rich gene ancestor, which emerged well before speciation in a common ancestor of all extant Drosophila species.     

Identification of One Intron Loss and Phylogenetic Evolution of Dfak Gene in the Drosophila melanogaster Species Group

Intron loss and its evolutionary significance have been noted in Drosophila. The current study provides another example of intron loss within a single-copy Dfak gene in Drosophila. By using polymerase chain reaction (PCR), we amplified about 1.3 kb fragment spanning intron 5–10, located in the position of Tyr kinase (TyK) domain of Dfak gene from Drosophila melanogaster species group, and observed size difference among the amplified DNA fragments from different species. Further sequencing analysis revealed that D. melanogaster and D. simulans deleted an about 60 bp of DNA fragment relative to other 7 Drosophila species, such as D. elegans, D. ficusphila, D. biarmipes, D. takahashii, D. jambulina, D. prostipennis and D. pseudoobscura, and the deleted fragment located precisely in the position of one intron. The data suggested that intron loss might have occurred in the Dfak gene evolutionary process of D. melanogaster and D. simulans of Drosophila melanogaster species group. In addition, the constructed phylogenetic tree based on the Dfak TyK domains clearly revealed the evolutionary relationships between subgroups of Drosophila melanogaster species group, and the intron loss identified from D. melanogaster and D. simulans provides a unique diagnostic tool for taxonomic classification of the melanogaster subgroup from other group of genus Drosophila.     

Reconstruction of ancestral karyotype illuminates chromosome evolution in the genus Cucumis

Karyotype dynamics driven by complex chromosome rearrangements constitute a fundamental issue in evolutionary genetics. The evolutionary events underlying karyotype diversity within plant genera, however, have rarely been reconstructed from a computed ancestral progenitor. Here, we developed a method to rapidly and accurately represent extant karyotypes with the genus, Cucumis, using highly customizable comparative oligo-painting (COP) allowing visualization of fine-scale genome structures of eight Cucumis species from both African-origin and Asian-origin clades. Based on COP data, an evolutionary framework containing a genus-level ancestral karyotype was reconstructed, allowing elucidation of the evolutionary events that account for the origin of these diverse genomes within Cucumis. Our results characterize the cryptic rearrangement hotspots on ancestral chromosomes, and demonstrate that the ancestral Cucumis karyotype (n = 12) evolved to extant Cucumis genomes by hybridizations and frequent lineage- and species-specific genome reshuffling. Relative to the African species, the Asian species, including melon (Cucumis melo, n = 12), Cucumis hystrix (n = 12) and cucumber (Cucumis sativus, n = 7), had highly shuffled genomes caused by large-scale inversions, centromere repositioning and chromothripsis-like rearrangement. The deduced reconstructed ancestral karyotype for the genus allowed us to propose evolutionary trajectories and specific events underlying the origin of these Cucumis species. Our findings highlight that the partitioned evolutionary plasticity of Cucumis karyotype is primarily located in the centromere-proximal regions marked by rearrangement hotspots, which can potentially serve as a reservoir for chromosome evolution due to their fragility.     

Genetics on the Fly: A Primer on the Drosophila Model System

Fruit flies of the genus Drosophila have been an attractive and effective genetic model organism since Thomas Hunt Morgan and colleagues made seminal discoveries with them a century ago. Work with Drosophila has enabled dramatic advances in cell and developmental biology, neurobiology and behavior, molecular biology, evolutionary and population genetics, and other fields. With more tissue types and observable behaviors than in other short-generation model organisms, and with vast genome data available for many species within the genus, the fly’s tractable complexity will continue to enable exciting opportunities to explore mechanisms of complex developmental programs, behaviors, and broader evolutionary questions. This primer describes the organism’s natural history, the features of sequenced genomes within the genus, the wide range of available genetic tools and online resources, the types of biological questions Drosophila can help address, and historical milestones.     

Characterization of a Tc1-like transposable element in zebrafish (Danio rerio)

We have characterized Tdr1, a family of Tc1-like transposable elements found in the genome of zebrafish (Danio rerio). The copy number and distribution of the sequence in the zebrafish genome have been determined, and by these criteria Tdr1 can be classified as a moderately repetitive, interspersed element. Examination of the sequences and structures of several copies of Tdr1 revealed that a particular deletion derivative, 1250 by long, of the transposon has been amplified to become the dominant form of Tdr1. The deletion in these elements encompasses sequences encoding the N-terminal portion of the putative Tdr1 transposase. Sequences corresponding to the deleted region were also detected, and thus allowed prediction of the nucleotide sequence of a hypothetical full-length element. Well conserved segments of Tc1-like transposons were found in the flanking regions of known fish genes, suggesting that these elements have a long evolutionary history in piscine genomes. Tdr1 elements have long, 208 by inverted repeats, with a short DNA motif repeated four times at the termini of the inverted repeats. Although different from that of the prototype C. elegans transposon Tc1, this inverted repeat structure is shared by transposable elements from salmonid fish species and two Drosophila species. We propose that these transposons form a subgroup within the Tc1-like family. Comparison of Tc1-like transposons supports the hypothesis that the transposase genes and their flanking sequences have been shaped by independent evolutionary constraints. Although Tc1-like sequences are present in the genomes of several strains of zebrafish and in salmonid fishes, these sequences are not conserved in the genus Danio, thus raising the possibility that these elements can be exploited for gene tagging and genome mapping.     

Chloroplast Genome Evolution in the Euglenaceae

Over the last few years multiple studies have been published outlining chloroplast genomes that represent many of the photosynthetic euglenid genera. However, these genomes were scattered throughout the euglenophyceaean phylogenetic tree, and focused on comparisons with Euglena gracilis. Here, we present a study exclusively on taxa within the Euglenaceae. Six new chloroplast genomes were characterized, those of Cryptoglena skujai, E. gracilis var. bacillaris, Euglena viridis, Euglenaria anabaena, Monomorphina parapyrum, and Trachelomonas volvocina, and added to six previously published chloroplast genomes to determine if trends existed within the family. With this study: at least one genome has now been characterized for each genus, the genomes of different strains from two taxa were characterized to explore intraspecific variability, and a second taxon has been characterized for the genus Monomorphina to examine intrageneric variability. Overall results showed a large amount of variability among the genomes, though a few trends could be identified both within Euglenaceae and within Euglenophyta. In addition, the intraspecific analysis indicated that the similarity of a genome sequence between strains was taxon dependent, and the intrageneric analysis indicated that the majority of the evolutionary changes within the Euglenaceae occurred intergenerically.     

Molecular Evolution of Drosophila Cuticular Protein Genes

Several multigene families have been described that together encode scores of structural cuticular proteins in Drosophila, although the functional significance of this diversity remains to be explored. Here I investigate the evolutionary histories of several multigene families (CPR, Tweedle, CPLCG, and CPF/CPFL) that vary in age, size, and sequence complexity, using sequenced Drosophila genomes and mosquito outgroups. My objective is to describe the rates and mechanisms of ‘cuticle-ome’ divergence, in order to identify conserved and rapidly evolving elements. I also investigate potential examples of interlocus gene conversion and concerted evolution within these families during Drosophila evolution. The absolute rate of change in gene number (per million years) is an order of magnitude lower for cuticular protein families within Drosophila than it is among Drosophila and the two mosquito taxa, implying that major transitions in the cuticle proteome have occurred at higher taxonomic levels. Several hotspots of intergenic conversion and/or gene turnover were identified, e.g. some gene pairs have independently undergone intergenic conversion within different lineages. Some gene conversion hotspots were characterized by conversion tracts initiating near nucleotide repeats within coding regions, and similar repeats were found within concertedly evolving cuticular protein genes in Anopheles gambiae. Rates of amino-acid substitution were generally severalfold higher along the branch connecting the Sophophora and Drosophila species groups, and 13 genes have Ka/Ks significantly greater than one along this branch, indicating adaptive divergence. Insect cuticular proteins appear to be a source of adaptive evolution within genera and, at higher taxonomic levels, subject to periods of gene-family expansion and contraction followed by quiescence. However, this relative stasis is belied by hotspots of molecular evolution, particularly concerted evolution, during the diversification of Drosophila. The prominent association between interlocus gene conversion and repeats within the coding sequence of interacting genes suggests that the latter promote strand exchange.     

Variation of the Mitochondral Genome in the Evolution of Drosophila

The evidence on mitochondrial genome variation and its role in evolution of the genus Drosophila are reviewed. The mitochondrial genome is represented by a circular double-stranded DNA molecule 16 to 19 kb in length. Mitochondrial genes lack introns and recombination. The entire mitochondrial genome can be arbitrarily divided into three parts: (1) protein-coding genes; (2) genes encoding rRNA and tRNA; and (3) the noncoding regulatory region (A + T region). The selective importance of mutations within different mtDNA regions is therefore unequal. In Drosophila, the content of the A + T pairs in mtDNA is extremely high and a pattern of nucleotide substitution is characterized by a low transition/transversion ratio (and a low threshold of mutation saturation). The deletions and duplications are of common occurrence in the mitochondrial genome. However, this genome lacks such characteristic for the nuclear genome aberrations as inversions and transpositions. The phenomena of introgression and heteroplasmy provide an opportunity to study the adaptive role of the mitochondrial genome and its role in speciation. Analysis of evidence concerning mtDNA variation in different species of the genus Drosophilamade it possible to ascertain data on phylogenetic relationships among species obtained by studying nuclear genome variation. In some species, mtDNA variation may serve as a reliable marker for population differentiation within a species, although evidence on the population dynamics of the mtDNA variation is very scarce.     

Out of the testis,into the ovary: biased outcomes of gene duplication and deletion in Drosophila

Gene turnover is a key source of adaptive variation. Yet most evolutionary studies have focused on gene duplication, dismissing gene deletion as a mechanism that simply eradicates redundancy. Here, I use genome‐scale sequence and multi‐tissue expression data from Drosophila melanogaster and Drosophila pseudoobscura to simultaneously assess the evolutionary outcomes of gene duplication and deletion in Drosophila. I find that gene duplication is more frequent than gene deletion in both species, indicating that it may play a more important role in Drosophila evolution. However, examination of several genic properties reveals that genes likely possess distinct functions after duplication that diverge further before deletion, suggesting that loss of redundancy cannot explain a majority of gene deletion events in Drosophila. Moreover, in addition to providing support for the well‐known “out of the testis” origin of young duplicate genes, analyses of gene expression profiles uncover a preferential bias against deletion of old ovary‐expressed genes. Therefore, I propose a novel “into the ovary” hypothesis for gene deletion in Drosophila, in which gene deletion may promote adaptation by salvaging genes that contribute to the evolution of female reproductive phenotypes. Under this combined “out of the testis, into the ovary” evolutionary model, gene duplication and deletion work in concert to generate and maintain a balanced repertoire of genes that promote sex‐specific adaptation in Drosophila.     

MITOCHONDRIAL–NUCLEAR EPISTASIS AFFECTS FITNESS WITHIN SPECIES BUT DOES NOT CONTRIBUTE TO FIXED INCOMPATIBILITIES BETWEEN SPECIES OF DROSOPHILA

Efficient mitochondrial function requires physical interactions between the proteins encoded by the mitochondrial and nuclear genomes. Coevolution between these genomes may result in the accumulation of incompatibilities between divergent lineages. We test whether mitochondrial–nuclear incompatibilities have accumulated within the Drosophila melanogaster species subgroup by combining divergent mitochondrial and nuclear lineages and quantifying the effects on relative fitness. Precise placement of nine mtDNAs from D. melanogaster, D. simulans, and D. mauritiana into two D. melanogaster nuclear genetic backgrounds reveals significant mitochondrial–nuclear epistasis affecting fitness in females. Combining the mitochondrial genomes with three different D. melanogaster X chromosomes reveals significant epistasis for male fitness between X‐linked and mitochondrial variation. However, we find no evidence that the more than 500 fixed differences between the mitochondrial genomes of D. melanogaster and the D. simulans species complex are incompatible with the D. melanogaster nuclear genome. Rather, the interactions of largest effect occur between mitochondrial and nuclear polymorphisms that segregate within species of the D. melanogaster species subgroup. We propose that a low mitochondrial substitution rate, resulting from a low mutation rate and/or efficient purifying selection, precludes the accumulation of mitochondrial–nuclear incompatibilities among these Drosophila species.     

Conservatism of sites of tRNA loci among the linkage groups of severalDrosophila species

Summary The sites of seven tRNA genes (Arg-2, Lys-2, Ser-2b, Ser-7, Thr-3, Thr-4, Val-3b) were studied by in situ hybridization.¹²⁵I-labeled tRNA probes fromDrosophila melanogaster were hybridized to spreads of polytene chromosomes prepared from fourDrosophila species representing different evolutionary lineages (D. melanogaster, Drosophila hydei, Drosophila pseudoobscura, andDrosophila virilis). Most tRNA loci occurred on homologous chromosomal elements of all four species. In some cases the number of hybridization sites within an element varied and sites on nonhomologous elements were found. It was observed that both tRNA ₂ ^Arg and tRNA ₂ ^Lys hybridized to the same site on homologous elements in several species. These data suggest a limited amount of exchange among different linkage groups during the evolution ofDrosophila species.     

Gene family evolution in green plants with emphasis on the origination and evolution of Arabidopsis thaliana genes

Gene family size variation is an important mechanism that shapes the natural variation for adaptation in various species. Despite its importance, the pattern of gene family size variation in green plants is still not well understood. In particular, the evolutionary pattern of genes and gene families remains unknown in the model plant Arabidopsis thaliana in the context of green plants. In this study, eight representative genomes of green plants are sampled to study gene family evolution and characterize the origination of A. thaliana genes, respectively. Four important insights gained are that: (i) the rate of gene gains and losses is about 0.001359 per gene every million years, similar to the rate in yeast, Drosophila, and mammals; (ii) some gene families evolved rapidly with extreme expansions or contractions, and 2745 gene families present in all the eight species represent the ‘core’ proteome of green plants; (iii) 70% of A. thaliana genes could be traced back to 450 million years ago; and (iv) intriguingly, A. thaliana genes with early origination are under stronger purifying selection and more conserved. In summary, the present study provides genome‐wide insights into evolutionary history and mechanisms of genes and gene families in green plants and especially in A. thaliana.     

Differential impact of retrotransposon populations on the genome of allotetraploid tobacco (<Emphasis Type="Italic">Nicotiana tabacum</Emphasis>)

LTR-retrotransposons contribute substantially to the structural diversity of plant genomes. Recent models of genome evolution suggest that retrotransposon amplification is offset by removal of retrotransposon sequences, leading to a turnover of retrotransposon populations. While bursts of amplification have been documented, it is not known whether removal of retrotransposon sequences occurs continuously, or is triggered by specific stimuli over short evolutionary periods. In this work, we have characterized the evolutionary dynamics of four populations of copia-type retrotransposons in allotetraploid tobacco (Nicotiana tabacum) and its two diploid progenitors Nicotiana sylvestris and Nicotiana tomentosiformis. We have used SSAP (Sequence-Specific Amplification Polymorphism) to evaluate the contribution retrotransposons have made to the diversity of tobacco and its diploid progenitor species, to quantify the contribution each diploid progenitor has made to tobacco's retrotransposon populations, and to estimate losses or amplifications of retrotransposon sequences subsequent to tobacco's formation. Our results show that the tobacco genome derives from a turnover of retrotransposon sequences with removals concomitant with new insertions. We have detected unique behaviour specific to each retrotransposon population, with differences likely reflecting distinct evolutionary histories and activities of particular elements. Our results indicate that the retrotransposon content of a given plant species is strongly influenced by the host evolutionary history, with periods of rapid turnover of retrotransposon sequences stimulated by allopolyploidy.