Search

Abstract Background Retrogenes are processed copies of other genes. This duplication mechanism produces a copy of the parental gene that should not contain introns, and usually does not contain cis-regulatory regions. Here, we computationally address the evolutionary origin of promoter and other cis-regulatory regions in retrogenes using a total of 94 Drosophila retroposition events we recently identified. Previous tissue expression data has revealed that a large fraction of these retrogenes are specifically and/or highly expressed in adult testes of Drosophila. Results In this work, we infer that retrogenes do not generally carry regulatory regions from aberrant upstream or normal transcripts of their parental genes, and that expression patterns of neighboring genes are not consistently shared by retrogenes. Additionally, transposable elements do not appear to substantially provide regulatory regions to retrogenes. Interestingly, we find that there is an excess of retrogenes in male testis neighborhoods that is not explained by insertional biases of the retroelement machinery used for retroposition. Conclusion We conclude that retrogenes' regulatory regions mostly do not represent a random set of existing regulatory regions. On the contrary, our conclusion is that selection is likely to have played an important role in the persistence of autosomal testis biased retrogenes. Selection in favor of retrogenes inserted in male testis neighborhoods and at the sequence level to produce testis expression is postulated to have occurred.

Evolutionary conflicts occur when there is antagonistic selection between different individuals of the same or different species, life stages or between levels of biological organization. Remarkably, conflicts can occur within species or within genomes. In the dynamics of evolutionary conflicts, gene duplications can play a major role because they can bring very specific changes to the genome: changes in protein dose, the generation of novel paralogues with different functions or expression patterns or the evolution of small antisense RNAs. As we describe here, by having those effects, gene duplication might spark evolutionary conflict or fuel arms race dynamics that takes place during conflicts. Interestingly, gene duplication can also contribute to the resolution of a within-locus evolutionary conflict by partitioning the functions of the gene that is under an evolutionary trade-off. In this review, we focus on intraspecific conflicts, including sexual conflict and illustrate the various roles of gene duplications with a compilation of examples. These examples reveal the level of complexity and the differences in the patterns of gene duplications within genomes under different conflicts. These examples also reveal the gene ontologies involved in conflict and the genomic location of the elements of the conflict. The examples provide a blueprint for the direct study of these conflicts or the exploration of the presence of similar conflicts in other lineages.

Background: Transposable elements (TEs) are selfish DNA sequences capable of moving and amplifying at the expense of host cells. Despite this, an increasing number of studies have revealed that TE proteins are important contributors to the emergence of novel host proteins through molecular domestication. We previously described seven transposase-derived domesticated genes from the PIF/Harbinger DNA family of TEs in Drosophila and a co-domestication. All PIF TEs known in plants and animals distinguish themselves from other DNA transposons by the presence of two genes. We hypothesize that there should often be co-domestications of the two genes from the same TE because the transposase (gene 1) has been described to be translocated to the nucleus by the MADF protein (gene 2). To provide support for this model of new gene origination, we investigated available insect species genomes for additional evidence of PIF TE domestication events and explored the co-domestication of the MADF protein from the same TE insertion., Results: After the extensive insect species genomes exploration of hits to PIF transposases and analyses of their context and evolution, we present evidence of at least six independent PIF transposable elements proteins domestication events in insects: two co-domestications of both transposase and MADF proteins in Anopheles (Diptera), one transposase-only domestication event and one co-domestication in butterflies and moths (Lepidoptera), and two transposases-only domestication events in cockroaches (Blattodea). The predicted nuclear localization signals for many of those proteins and dicistronic transcription in some instances support the functional associations of co-domesticated transposase and MADF proteins., Conclusions: Our results add to a co-domestication that we previously described in fruit fly genomes and support that new gene origination through domestication of a PIF transposase is frequently accompanied by the co-domestication of a cognate MADF protein in insects, potentially for regulatory functions. We propose a detailed model that predicts that PIF TE protein co-domestication should often occur from the same PIF TE insertion., (© 2022. The Author(s).)

Genes that originate during evolution are an important source of novel biological functions. Retrogenes are functional copies of genes produced by retroduplication and as such are located in different genomic positions. To investigate retroposition patterns and retrogene expression, we computationally identified interchromosomal retroduplication events in nine portions of the phylogenetic history of malaria mosquitoes, making use of species that do or do not have classical sex chromosomes to test the roles of sex-linkage. We found 40 interchromosomal events and a significant excess of retroduplications from the X chromosome to autosomes among a set of young retrogenes. These young retroposition events occurred within the last 100 million years in lineages where all species possessed differentiated sex chromosomes. An analysis of available microarray and RNA-seq expression data for Anopheles gambiae showed that many of the young retrogenes evolved male-biased expression in the reproductive organs. Young autosomal retrogenes with increased meiotic or postmeiotic expression in the testes tend to be male biased. In contrast, older retrogenes, i.e., in lineages with undifferentiated sex chromosomes, do not show this particular chromosomal bias and are enriched for female-biased expression in reproductive organs. Our reverse-transcription PCR data indicates that most of the youngest retrogenes, which originated within the last 47.6 million years in the subgenus Cellia , evolved non-uniform expression patterns across body parts in the males and females of An. coluzzii . Finally, gene annotation revealed that mitochondrial function is a prominent feature of the young autosomal retrogenes. We conclude that mRNA-mediated gene duplication has produced a set of genes that contribute to mosquito reproductive functions and that different biases are revealed after the sex chromosomes evolve. Overall, these results suggest potential roles for the evolution of meiotic sex chromosome inactivation in males and of sexually antagonistic conflict related to mitochondrial energy function as the main selective pressures for X-to-autosome gene reduplication and testis-biased expression in these mosquito lineages.

Recent studies on nuclear-encoded mitochondrial genes (N-mt genes) in Drosophila melanogaster have shown a unique pattern of expression for newly duplicated N-mt genes, with many duplicates having a testis-biased expression and playing an essential role in spermatogenesis. In this study, we investigated a newly duplicated N-mt gene-i.e., Cytochrome c oxidase 4-like ( COX4L )-in order to understand its function and, consequently, the reason behind its retention in the D. melanogaster genome. The COX4L gene is a duplicate of the Cytochrome c oxidase 4 ( COX4 ) gene of OXPHOS complex IV. While the parental COX4 gene has been found in all eukaryotes, including single-cell eukaryotes such as yeast, we show that COX4L is only present in the Brachycera suborder of Diptera; thus, both genes are present in all Drosophila species, but have significantly different patterns of expression: COX4 is highly expressed in all tissues, while COX4L has a testis-specific expression. To understand the function of this new gene, we first knocked down its expression in the D. melanogaster germline using two different RNAi lines driven by the bam-Gal4 driver; second, we created a knockout strain for this gene using CRISPR-Cas9 technology. Our results showed that knockdown and knockout lines of COX4L produce partial sterility and complete sterility in males, respectively, where a lack of sperm individualization was observed in both cases. Male infertility was prevented by driving COX4L -HA in the germline, but not when driving COX4 -HA. In addition, ectopic expression of COX4L in the soma caused embryonic lethality, while overexpression in the germline led to a reduction in male fertility. COX4L -KO mitochondria show reduced membrane potential, providing a plausible explanation for the male sterility observed in these flies. This prominent loss-of-function phenotype, along with its testis-biased expression and its presence in the Drosophila sperm proteome, suggests that COX4L is a paralogous, specialized gene that is assembled in OXPHOS complex IV of male germline cells and/or sperm mitochondria.

Background: The nuclear transport machinery is involved in a well-known male meiotic drive system in Drosophila. Fast gene evolution and gene duplications have been major underlying mechanisms in the evolution of meiotic drive systems, and this might include some nuclear transport genes in Drosophila. So, using a comprehensive, detailed phylogenomic study, we examined 51 insect genomes for the duplication of the same nuclear transport genes., Results: We find that most of the nuclear transport duplications in Drosophila are of a few classes of nuclear transport genes, RNA mediated and fast evolving. We also retrieve many pseudogenes for the Ran gene. Some of the duplicates are relatively young and likely contributing to the turnover expected for genes under strong but changing selective pressures. These duplications are potentially revealing what features of nuclear transport are under selection. Unlike in flies, we find only a few duplications when we study the Drosophila duplicated nuclear transport genes in dipteran species outside of Drosophila, and none in other insects., Conclusions: These findings strengthen the hypothesis that nuclear transport gene duplicates in Drosophila evolve either as drivers or suppressors of meiotic drive systems or as other male-specific adaptations circumscribed to flies and involving a handful of nuclear transport functions., (© 2021. The Author(s).)

Emerging evidence suggests that ribosome heterogeneity may have important functional consequences in the translation of specific mRNAs within different cell types and under various conditions. Ribosome heterogeneity comes in many forms, including post-translational modification of ribosome proteins (RPs), absence of specific RPs and inclusion of different RP paralogs. The Drosophila genome encodes two RpS5 paralogs: RpS5a and RpS5b. While RpS5a is ubiquitously expressed, RpS5b exhibits enriched expression in the reproductive system. Deletion of RpS5b results in female sterility marked by developmental arrest of egg chambers at stages 7-8, disruption of vitellogenesis and posterior follicle cell (PFC) hyperplasia. While transgenic rescue experiments suggest functional redundancy between RpS5a and RpS5b, molecular, biochemical and ribo-seq experiments indicate that RpS5b mutants display increased rRNA transcription and RP production, accompanied by increased protein synthesis. Loss of RpS5b results in microtubule-based defects and in mislocalization of Delta and Mindbomb1, leading to failure of Notch pathway activation in PFCs. Together, our results indicate that germ cell-specific expression of RpS5b promotes proper egg chamber development by ensuring the homeostasis of functional ribosomes., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2021. Published by The Company of Biologists Ltd.)

Linear chromosomes shorten in every round of replication. In Drosophila, telomere-specialized long interspersed retrotransposable elements (LINEs) belonging to the jockey clade offset this shortening by forming head-to-tail arrays at Drosophila telomere ends. As such, these telomeric LINEs have been considered adaptive symbionts of the genome, protecting it from premature decay, particularly as Drosophila lacks a conventional telomerase holoenzyme. However, as reviewed here, recent work reveals a high degree of variation and turnover in the telomere-specialized LINE lineages across Drosophila. There appears to be no absolute requirement for LINE activity to maintain telomeres in flies, hence the suggestion that the telomere-specialized LINEs may instead be neutral or in conflict with the host, rather than adaptive., (© 2019 WILEY Periodicals, Inc.)

Here we use a chromosome-level genome assembly of a prairie rattlesnake ( Crotalus viridis ), together with Hi-C, RNA-seq, and whole-genome resequencing data, to study key features of genome biology and evolution in reptiles. We identify the rattlesnake Z Chromosome, including the recombining pseudoautosomal region, and find evidence for partial dosage compensation driven by an evolutionary accumulation of a female-biased up-regulation mechanism. Comparative analyses with other amniotes provide new insight into the origins, structure, and function of reptile microchromosomes, which we demonstrate have markedly different structure and function compared to macrochromosomes. Snake microchromosomes are also enriched for venom genes, which we show have evolved through multiple tandem duplication events in multiple gene families. By overlaying chromatin structure information and gene expression data, we find evidence for venom gene-specific chromatin contact domains and identify how chromatin structure guides precise expression of multiple venom gene families. Further, we find evidence for venom gland-specific transcription factor activity and characterize a complement of mechanisms underlying venom production and regulation. Our findings reveal novel and fundamental features of reptile genome biology, provide insight into the regulation of snake venom, and broadly highlight the biological insight enabled by chromosome-level genome assemblies., (© 2019 Schield et al.; Published by Cold Spring Harbor Laboratory Press.)

Transposable elements (TEs) are selfish genetic units that typically encode proteins that enable their proliferation in the genome and spread across individual hosts. Here we review a growing number of studies that suggest that TE proteins have often been co-opted or 'domesticated' by their host as adaptations to a variety of evolutionary conflicts. In particular, TE-derived proteins have been recurrently repurposed as part of defense systems that protect prokaryotes and eukaryotes against the proliferation of infectious or invasive agents, including viruses and TEs themselves. We argue that the domestication of TE proteins may often be the only evolutionary path toward the mitigation of the cost incurred by their own selfish activities., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

Gene duplication is a major driver of organismal evolution. Gene retroposition is a mechanism of gene duplication whereby a gene's transcript is used as a template to generate retroposed gene copies, or retrocopies. Intriguingly, the formation of retrocopies depends upon the enzymatic machinery encoded by retrotransposable elements, genomic parasites occurring in the majority of eukaryotes. Most retrocopies are depleted of the regulatory regions found upstream of their parental genes; therefore, they were initially considered transcriptionally incompetent gene copies, or retropseudogenes. However, examples of functional retrocopies, or retrogenes, have accumulated since the 1980s. Here, we review what we have learned about retrocopies in animals, plants and other eukaryotic organisms, with a particular emphasis on comparative and population genomic analyses complemented with transcriptomic datasets. In addition, these data have provided information about the dynamics of the different "life cycle" stages of retrocopies (i.e., polymorphic retrocopy number variants, fixed retropseudogenes and retrogenes) and have provided key insights into the retroduplication mechanisms, the patterns and evolutionary forces at work during the fixation process and the biological function of retrogenes. Functional genomic and transcriptomic data have also revealed that many retropseudogenes are transcriptionally active and a biological role has been experimentally determined for many. Finally, we have learned that not only non-long terminal repeat retroelements but also long terminal repeat retroelements play a role in the emergence of retrocopies across eukaryotes. This body of work has shown that mRNA-mediated duplication represents a widespread phenomenon that produces an array of new genes that contribute to organismal diversity and adaptation., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Several studies in Drosophila have shown a paucity of male-biased genes (i.e., genes that express higher in males than in females) on the X chromosome. Dosage compensation (DC) is a regulatory mechanism of gene expression triggered in males that hypertranscribes the X-linked genes to the level of transcription in females. There are currently two different hypotheses about the effects of DC on the distribution of male-biased genes: (1) it might limit male-expression level, or (2) it might interfere with the male upregulation of gene expression. Here, we used previously published gene expression datasets to reevaluate both hypotheses and introduce a mutually exclusive prediction that helped us to reject the hypothesis that the paucity of male-biased genes in the X chromosome is due to a limit in the male-expression level. Our analysis also uncovers unanticipated details about how DC interferes with the genomic distribution of both, male-biased and female-biased genes. We suggest that DC actually interferes with female downregulation of gene expression and not male upregulation, as previously suggested.

Cactophilic Drosophila species provide a valuable model to study gene-environment interactions and ecological adaptation. Drosophila buzzatii and Drosophila mojavensis are two cactophilic species that belong to the repleta group, but have very different geographical distributions and primary host plants. To investigate the genomic basis of ecological adaptation, we sequenced the genome and developmental transcriptome of D. buzzatii and compared its gene content with that of D. mojavensis and two other noncactophilic Drosophila species in the same subgenus. The newly sequenced D. buzzatii genome (161.5 Mb) comprises 826 scaffolds (>3 kb) and contains 13,657 annotated protein-coding genes. Using RNA sequencing data of five life-stages we found expression of 15,026 genes, 80% protein-coding genes, and 20% noncoding RNA genes. In total, we detected 1,294 genes putatively under positive selection. Interestingly, among genes under positive selection in the D. mojavensis lineage, there is an excess of genes involved in metabolism of heterocyclic compounds that are abundant in Stenocereus cacti and toxic to nonresident Drosophila species. We found 117 orphan genes in the shared D. buzzatii-D. mojavensis lineage. In addition, gene duplication analysis identified lineage-specific expanded families with functional annotations associated with proteolysis, zinc ion binding, chitin binding, sensory perception, ethanol tolerance, immunity, physiology, and reproduction. In summary, we identified genetic signatures of adaptation in the shared D. buzzatii-D. mojavensis lineage, and in the two separate D. buzzatii and D. mojavensis lineages. Many of the novel lineage-specific genomic features are promising candidates for explaining the adaptation of these species to their distinct ecological niches., (© The Author(s) 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Retrogenes are functional processed copies of genes that originate via the retrotranscription of an mRNA intermediate and often exhibit testis-specific expression. Although this expression pattern appears to be favored by selection, the origin of such expression bias remains unexplained. Here, we study the regulation of two young testis-specific Drosophila retrogenes, Dntf-2r and Pros28.1A, using genetic transformation and the enhanced green fluorescent protein reporter gene in Drosophila melanogaster. We show that two different short (<24 bp) regions upstream of the transcription start sites (TSSs) act as testis-specific regulatory motifs in these genes. The Dntf-2r regulatory region is similar to the known β2 tubulin 14-bp testis motif (β2-tubulin gene upstream element 1 [β2-UE1]). Comparative sequence analyses reveal that this motif was already present before the Dntf-2r insertion and was likely driving the transcription of a noncoding RNA. We also show that the β2-UE1 occurs in the regulatory regions of other testis-specific retrogenes, and is functional in either orientation. In contrast, the Pros28.1A testes regulatory region in D. melanogaster appears to be novel. Only Pros28.1B, an older paralog of the Pros28.1 gene family, seems to carry a similar regulatory sequence. It is unclear how the Pros28.1A regulatory region was acquired in D. melanogaster, but it might have evolved de novo from within a region that may have been preprimed for testes expression. We conclude that relocation is critical for the evolutionary origin of male germline-specific cis-regulatory regions of retrogenes because expression depends on either the site of the retrogene insertion or the sequence changes close to the TSS thereafter. As a consequence we infer that positive selection will play a role in the evolution of these regulatory regions and can often act from the moment of the retrocopy insertion., (© The Author 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Chromosomal inversions are usually portrayed as simple two-breakpoint rearrangements changing gene order but not gene number or structure. However, increasing evidence suggests that inversion breakpoints may often have a complex structure and entail gene duplications with potential functional consequences. Here, we used a combination of different techniques to investigate the breakpoint structure and the functional consequences of a complex rearrangement fixed in Drosophila buzzatii and comprising two tandemly arranged inversions sharing the middle breakpoint: 2m and 2n. By comparing the sequence in the breakpoint regions between D. buzzatii (inverted chromosome) and D. mojavensis (noninverted chromosome), we corroborate the breakpoint reuse at the molecular level and infer that inversion 2m was associated with a duplication of a ~13 kb segment and likely generated by staggered breaks plus repair by nonhomologous end joining. The duplicated segment contained the gene CG4673, involved in nuclear transport, and its two nested genes CG5071 and CG5079. Interestingly, we found that other than the inversion and the associated duplication, both breakpoints suffered additional rearrangements, that is, the proximal breakpoint experienced a microinversion event associated at both ends with a 121-bp long duplication that contains a promoter. As a consequence of all these different rearrangements, CG5079 has been lost from the genome, CG5071 is now a single copy nonnested gene, and CG4673 has a transcript ~9 kb shorter and seems to have acquired a more complex gene regulation. Our results illustrate the complex effects of chromosomal rearrangements and highlight the need of complementing genomic approaches with detailed sequence-level and functional analyses of breakpoint regions if we are to fully understand genome structure, function, and evolutionary dynamics.

We look at sex-limited chromosome (Y or W) evolution with particular emphasis on the importance of palindromes. Y chromosome palindromes consist of inverted duplicates that allow for local recombination in an otherwise nonrecombining chromosome. Since palindromes enable intrachromosomal gene conversion that can help eliminate deleterious mutations, they are often highlighted as mechanisms to protect against Y degeneration. However, the adaptive significance of recombination resides in its ability to decouple the evolutionary fates of linked mutations, leading to both a decrease in degeneration rate and an increase in adaptation rate. Our paper emphasizes the latter, that palindromes may exist to accelerate adaptation by increasing the potential targets and fixation rates of incoming beneficial mutations. This hypothesis helps reconcile two enigmatic features of the "palindromes as protectors" view: (1) genes that are not located in palindromes have been retained under purifying selection for tens of millions of years, and (2) under models that only consider deleterious mutations, gene conversion benefits duplicate gene maintenance but not initial fixation. We conclude by looking at ways to test the hypothesis that palindromes enhance the rate of adaptive evolution of Y-linked genes and whether this effect can be extended to palindromes on other chromosomes.

Gene duplication is mainly recognized by its primary role in the origin of new genes and functions. However, the idea that gene duplication can be a central player in resolving sexual genetic conflicts through its potential to generate sex-biased and sex-specifically expressed genes, has been almost entirely overlooked. We review recent data and theory that support gene duplication as a theoretically predicted and experimentally supported means of resolving intralocus sexual antagonism. We believe that this role is probably the consequence of sexual conflict for housekeeping genes that are required in males and females, and which are expressed in sexually dimorphic tissues (i.e. where sexually antagonistic selection is exerted). We think that these genes cannot evolve tissue-specific expression unless they duplicate., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

In species that have two sexes, a single genome encodes two morphs, as each sex can be thought of as a distinct morph. This means that the same set of genes are differentially expressed in the different sexes. Many questions emanate from this statement. What proportion of genes contributes to sexual dimorphism? How do they contribute to sexual dimorphism? How is sex-biased expression achieved? Which sex and what tissues contribute the most to sex-biased expression? Do sex-biased genes have the same evolutionary patterns as nonbiased genes? We review the current data on sex-biased expression in species with heteromorphic sex chromosomes and comment on the most important hypotheses suggested to explain the origin, evolution, and distribution patterns of sex-biased genes. In this perspective we emphasize how gene duplication serves as an important molecular mechanism to resolve genomic clashes and genetic conflicts by generating sex-biased genes, often sex-specific genes, and contributes greatly to the underlying genetic basis of sexual dimorphism.

Gene duplication is a major force driving genome evolution, and examples of this mode of evolution and of the functions of duplicated genes are needed to reveal general patterns. Here, our study focuses on a particular retrogene (i.e., CG9573) that originated about 5-13 million years ago that we have named Drcd-1 related. It originated in Drosophila through retroposition of the parental gene Required for cell differentiation 1 of Drosophila (Drcd-1; CG14213), which is a known transcription cofactor. Drcd-1r is only present in D. melanogaster, D. simulans, D. sechellia, and D. mauritiana. Drcd-1r is an X to autosome retroposition event. Many retrogenes are X to autosome copies and it has been shown that positive selection underlies this bias. We sought to understand Drcd-1r mode of evolution and function to contribute to the understanding of the selective pressures acting on X to autosome retrogenes. Drcd-1r overlaps with another gene, it is within the 3' UTR of the gene CG13102 and is encoded in the opposite orientation. We have studied the characteristics of the transcripts and quantified expression of CG13102 and Drcd-1r in wild-type flies. We found that Drcd-1r is transcribed specifically in testes. We also studied the molecular evolution of Drcd-1r and Drcd-1 and found that the parental gene has evolved under very strong purifying selection but the retrogene has evolved very rapidly (Ka/Ks ~1) under both positive and purifying selection, as revealed using divergence and polymorphism data. These results indicate that Drcd-1r has a novel function in the Drosophila testes. To further explore Drcd-1r function we used a strain containing a P element inserted in the region where CG13102 and Drcd-1r are located that shows recessive male sterility. Analysis of this strain reveals the difficulties that can be encountered in studying the functions of genes with overlapping transcripts. Avenues for studying of the function of this gene are proposed.

The analyses of gene duplications by retroposition have revealed an excess of male-biased duplicates generated from X chromosome to autosomes in flies and mammals. Investigating these genes is of primary importance in understanding sexual dimorphism and genome evolution. In a particular instance in Drosophila, X-linked nuclear transport genes (Ntf-2 and ran) have given rise to autosomal retroposed copies three independent times (along the lineages leading to Drosophila melanogaster, D. ananassae, and D. grimshawi). Here we explore in further detail the expression and the mode of evolution of these Drosophila Ntf-2- and ran-derived retrogenes. Five of the six retrogenes show male-biased expression. The ran-like gene of D. melanogaster and D. simulans has undergone recurrent positive selection. Similarly, in D. ananassae and D. atripex, the Ntf-2 and ran retrogenes show evidence of past positive selection. The data suggest that strong selection is acting on the origin and evolution of these retrogenes. Avoiding male meiotic X inactivation, increasing level of expression of X-linked genes in male testes, and/or sexual antagonism might explain the recurrent duplication of retrogenes from X to autosomes. Interestingly, the ran-like in D. yakuba has mostly pseudogenized alleles. Disablement of the ran-like gene in D. yakuba indicates turnover of these duplicates. We discuss the possibility that Dntf-2r and ran-like might be involved in genomic conflicts during spermatogenesis.

Gene duplication is probably the most important mechanism for generating new gene functions. However, gene duplication has been overlooked as a potentially effective way to resolve genetic conflicts. Here, we analyze the entire set of Drosophila melanogaster nuclearly encoded mitochondrial duplicate genes and show that both RNA- and DNA-mediated mitochondrial gene duplications exhibit an unexpectedly high rate of relocation (change in location between parental and duplicated gene) as well as an extreme tendency to avoid the X chromosome. These trends are likely related to our observation that relocated genes tend to have testis-specific expression. We also infer that these trends hold across the entire Drosophila genus. Importantly, analyses of gene ontology and functional interaction networks show that there is an overrepresentation of energy production-related functions in these mitochondrial duplicates. We discuss different hypotheses to explain our results and conclude that our findings substantiate the hypothesis that gene duplication for male germline function is likely a mechanism to resolve intralocus sexually antagonistic conflicts that we propose are common in testis. In the case of nuclearly encoded mitochondrial duplicates, our hypothesis is that past sexually antagonistic conflict related to mitochondrial energy function in Drosophila was resolved by gene duplication.

Proteasomes are large multisubunit complexes responsible for regulated protein degradation. Made of a core particle (20S) and regulatory caps (19S), proteasomal proteins are encoded by at least 33 genes, of which 12 have been shown to have testis-specific isoforms in Drosophila melanogaster. Pros28.1A (also known as Prosalpha4T1), a young retroduplicate copy of Pros28.1 (also known as Prosalpha4), is one of these isoforms. It is present in the D. melanogaster subgroup and was previously shown to be testis-specific in D. melanogaster. Here, we show its testis-specific transcription in all D. melanogaster subgroup species. Due to this conserved pattern of expression in the species harboring this insertion, we initially expected that a regulatory region common to these species evolved prior to the speciation event. We determined that the region driving testis expression in D. melanogaster is not far from the coding region (within 272 bp upstream of the ATG). However, different Transcription Start Sites (TSSs) are used in D. melanogaster and D. simulans, and a "broad" transcription start site is used in D. yakuba. These results suggest one of the following scenarios: (1) there is a conserved motif in the 5' region of the gene that can be used as an upstream or downstream element or at different distance depending on the species; (2) different species evolved diverse regulatory sequences for the same pattern of expression (i.e., "TSS turnover"); or (3) the transcription start site can be broad or narrow depending on the species. This work reveals the difficulties of studying gene regulation in one species and extrapolating those findings to close relatives.

Retrogenes are processed copies of genes that are inserted into new genomic regions and that acquire new regulatory elements from the sequences in their surroundings. Here we use a comparative approach of phylogenetic footprinting and a non-comparative approach of measuring motif over-representation in retrogenes in order to describe putative elements present in cis-regulatory regions of 94 retrogenes recently described in Drosophila. The detailed examination of the motifs found in the core promoter regions of retrogenes reveals an abundance of the DNA replication-related element (DRE), the Initiator (Inr), and a new over-represented motif that we call the GCT motif. Parental genes also show an abundance of DRE and Inr motifs, but these do not seem to have been carried over with retrogenes. In particular, we also examined motifs upstream of retrogenes expressed in adult testis and were able to identify 6 additional over-represented motifs. Comparative analyses provide data on the conservation and origin of some of these motifs and reveal 15 additional conserved motifs in these retrogenes. Some of those conserved motifs are sequences bound by known transcription factors, while others are novel motifs. In this report we provide the first genome-wide data on which specific cis-regulatory regions can be recruited by retrogenes after they are inserted into new coding regions in the genome. Future experiments are needed to determine the function and role of the new elements presented here.

The P instability factor or PIF superfamily of DNA transposons constitutes an important group of transposable elements (TEs) in plants, but it is still poorly characterized in metazoans. Taking advantage of the availability of draft genome sequences for twelve Drosophila species, we discovered 4 different lineages of Drosophila PIF-like transposons, named DPLT1-4. These lineages have experienced a complex evolutionary history during the Drosophila radiation, involving differential amplification and retention among species and probable events of horizontal transmission. Like previously described plant and animal PIF transposons, full-length DPLTs encode a putative transposase as well as a second predicted protein containing a Myb/SANT domain. In DPLTs, this domain is most closely related to the MADF DNA-binding domain found in several Drosophila transcription factors. In addition, we identified 7 distinct genes distributed across the Drosophila genus that encode proteins related to PIF transposases, but lack the hallmarks of transposons. Instead, these sequences show features of functional genes, such as an intact coding region evolving under purifying selection, the presence of orthologs in at least 2 Drosophila species, and the conservation of intron/exon structure across orthologs. We also provide evidence that most of these genes are transcribed and that some are developmentally regulated. Together the data indicate that these genes derived from PIF-transposons that have been "domesticated" to serve cellular functions. In one instance the recruitment of the transposase gene was accompanied by the co-recruitment of the adjacent second PIF gene, which raises the hypothesis that both proteins now function in the same pathway. The second PIF gene has retained the capacity to encode a protein with an intact MADF domain, suggesting that it may function as a transcription factor. We conclude that PIF transposons are common in the Drosophila lineage and have been a recurrent source of new genes during Drosophila evolution.

The strand-biased mutation spectrum in vertebrate mitochondrial genomes results in an AC-rich L-strand and a GT-rich H-strand. Because the L-strand is the sense strand of 12 protein-coding genes out of the 13, the third codon position is overall strongly AC-biased. The wobble site of the anticodon of the 22 mitochondrial tRNAs is either U or G to pair with the most abundant synonymous codon, with only one exception. The wobble site of Met-tRNA is C instead of U, forming the Watson-Crick match with AUG instead of AUA, the latter being much more frequent than the former. This has been attributed to a compromise between translation initiation and elongation; i.e., AUG is not only a methionine codon, but also an initiation codon, and an anticodon matching AUG will increase the initiation rate. However, such an anticodon would impose selection against the use of AUA codons because AUA needs to be wobble-translated. According to this translation conflict hypothesis, AUA should be used relatively less frequently compared to UUA in the UUR codon family. A comprehensive analysis of mitochondrial genomes from a variety of vertebrate species revealed a general deficiency of AUA codons relative to UUA codons. In contrast, urochordate mitochondrial genomes with two tRNA(Met) genes with CAU and UAU anticodons exhibit increased AUA codon usage. Furthermore, six bivalve mitochondrial genomes with both of their tRNA-Met genes with a CAU anticodon have reduced AUA usage relative to three other bivalve mitochondrial genomes with one of their two tRNA-Met genes having a CAU anticodon and the other having a UAU anticodon. We conclude that the translation conflict hypothesis is empirically supported, and our results highlight the fine details of selection in shaping molecular evolution.

Background: Processed copies of genes (retrogenes) are duplicate genes that originated through the reverse-transcription of a host transcript and insertion in the genome. This type of gene duplication, as any other, could be a source of new genes and functions. Using whole genome sequence data for 12 Drosophila species, we dated the origin of 94 retroposition events that gave rise to candidate functional genes in D. melanogaster., Results: Based on this analysis, we infer that functional retrogenes have emerged at a fairly constant rate of 0.5 genes per million years per lineage over the last approximately 63 million years of Drosophila evolution. The number of functional retrogenes and the rate at which they are recruited in the D. melanogaster lineage are of the same order of magnitude as those estimated in the human lineage, despite the higher deletion bias in the Drosophila genome. However, unlike primates, the rate of retroposition in Drosophila seems to be fairly constant and no burst of retroposition can be inferred from our analyses. In addition, our data also support an important role for retrogenes as a source of lineage-specific male functions, in agreement with previous hypotheses. Finally, we identified three cases of functional retrogenes in D. melanogaster that have been independently retroposed and recruited in parallel as new genes in other Drosophila lineages., Conclusion: Together, these results indicate that retroposition is a persistent mechanism and a recurrent pathway for the emergence of new genes in Drosophila.

This is the first detailed study of the evolution, phylogenetic distribution, and transcription of one young retroposed gene, CG13732, and its parental gene CG15645, whose functions are unknown. CG13732 is a recognizable retroposed copy of CG15645 retaining the signals of this process. We name the parental gene Cervantes and the retrogene Quijote. To determine when this duplication occurred and the phylogenetic distribution of Quijote, we employed polymerase chain reaction, Southern blotting, and the available information on sequenced Drosophila genomes. Interestingly, these analyses revealed that Quijote is present only in 4 species of Drosophila (Drosophila melanogaster, Drosophila simulans, Drosophila sechellia, and Drosophila mauritiana) and that retroposed copies of Cervantes have also originated in the lineages leading to Drosophila yakuba and Drosophila erecta independently in the 3 instances. We name the new retrogene in the D. yakuba lineage Rocinante and the new retrogene in the D. erecta lineage Sancho. In this work, we present data on Quijote and its parental gene Cervantes. Polymorphism analysis of the derived gene and divergence data for both parental and derived genes were used to determine that both genes likely produce functional proteins and that they are changing at a fast rate (KA/KS approximately 0.38). The negative value of H of Fay and Wu in the non-African sample reveals an excess of derived variants at high frequency. This could be explained either by positive selection in the region or by demographic effects. The comparative expression pattern shows that both genes express in the same adult tissues (male and female germ line) in D. melanogaster. Quijote is also expressed in male and female in D. simulans, D. sechellia, and D. mauritiana. We argue that the fast rate of evolution of these genes could be related to their putative germ line function and are further studying the independent recruitment of Cervantes-derived retrogenes in multiple lineages.

The position of a gene in the genome may have important consequences for its function. Therefore, when a new duplicate gene arises, its location may be critical in determining its fate. Our recent work in humans, mouse, and Drosophila provided a test by studying the patterns of duplication in sex chromosome evolution. We revealed a bias in the generation and recruitment of new gene copies involving the X chromosome that has been shaped largely by selection for male germline functions. The gene movement patterns we observed reflect an ongoing process as some of the new genes are very young while others were present before the divergence of humans and mouse. This suggests a continuing redistribution of male-related genes to achieve a more efficient allocation of male functions. This notion should be further tested in organisms employing other sex determination systems or in organisms differing in germline sex chromosome inactivation. It is likely that the selective forces that were detected in these studies are also acting on other types of duplicate genes. As a result, future work elucidating sex chromosome differentiation by other mutational mechanisms will shed light on this important process.

Mammalian sex chromosomes have undergone profound changes since evolving from ancestral autosomes. By examining retroposed genes in the human and mouse genomes, we demonstrate that, during evolution, the mammalian X chromosome has generated and recruited a disproportionately high number of functional retroposed genes, whereas the autosomes experienced lower gene turnover. Most autosomal copies originating from X-linked genes exhibited testis-biased expression. Such export is incompatible with mutational bias and is likely driven by natural selection to attain male germline function. However, the excess recruitment is consistent with a combination of both natural selection and mutational bias.

Genome data have revealed great variation in the numbers of genes in different organisms, which indicates that there is a fundamental process of genome evolution: the origin of new genes. However, there has been little opportunity to explore how genes with new functions originate and evolve. The study of ancient genes has highlighted the antiquity and general importance of some mechanisms of gene origination, and recent observations of young genes at early stages in their evolution have unveiled unexpected molecular and evolutionary processes.

Exon shuffling is an essential molecular mechanism for the formation of new genes. Many cases of exon shuffling have been reported in vertebrate genes. These discoveries revealed the importance of exon shuffling in the origin of new genes. However, only a few cases of exon shuffling were reported from plants and invertebrates, which gave rise to the assertion that the intron-mediated recombination mechanism originated very recently. We focused on the origin of new genes by exon shuffling and retroposition. We will first summarize our experimental work, which revealed four new genes in Drosophila, plants, and humans. These genes are 10(6) to 10(8) million years old. The recency of these genes allows us to directly examine the origin and evolution of genes in detail. These observations show firstly the importance of exon shuffling and retroposition in the rapid creation of new gene structures. They also show that the resultant chimerical structures appearing as mosaic proteins or as retroposed coding structures with novel regulatory systems, often confer novel functions. Furthermore, these newly created genes appear to have been governed by positive Darwinian selection throughout their history, with rapid changes of amino acid sequence and gene structure in very short periods of evolution. We further analyzed the distribution of intron phases in three non-vertebrate species, Drosophila melanogaster, Caenorhabditis elegans, and Arabidosis thaliana, as inferred from their genome sequences. As in the case of vertebrate genes, we found that intron phases in these species are unevenly distributed with an excess of phase zero introns and a significant excess of symmetric exons. Both findings are consistent with the requirements for the molecular process of exon shuffling. Thus, these non-vertebrate genomes may have also been strongly impacted by exon shuffling in general.

A direct approach to investigating new gene origination is to examine recently evolved genes. We report a new gene in the Drosophila melanogaster subgroup, Drosophila nuclear transport factor-2-related (Dntf-2r). Its sequence features and phylogenetic distribution indicate that Dntf-2r is a retroposed functional gene and originated in the common ancestor of D. melanogaster, D. simulans, D. sechellia, and D. mauritiana, within the past 3-12 million years (MY). Dntf-2r evolved more rapidly than the parental gene, under positive Darwinian selection as revealed by the McDonald-Kreitman test and other evolutionary analyses. Comparative expression analysis shows that Dntf-2r is male specific whereas the parental gene, Dntf-2, is widely expressed in D. melanogaster. In agreement with its new expression pattern, the Dntf-2r putative promoter sequence is similar to the late testis promoter of beta2-tubulin. We discuss the possibility that the action of positive selection in Dntf-2r is related to its putative male-specific functions.

New genes that originated by various molecular mechanisms are an essential component in understanding the evolution of genetic systems. We investigated the pattern of origin of the genes created by retroposition in Drosophila. We surveyed the whole Drosophila melanogaster genome for such new retrogenes and experimentally analyzed their functionality and evolutionary process. These retrogenes, functional as revealed by the analysis of expression, substitution, and population genetics, show a surprisingly asymmetric pattern in their origin. There is a significant excess of retrogenes that originate from the X chromosome and retropose to autosomes; new genes retroposed from autosomes are scarce. Further, we found that most of these X-derived autosomal retrogenes had evolved a testis expression pattern. These observations may be explained by natural selection favoring those new retrogenes that moved to autosomes and avoided the spermatogenesis X inactivation, and suggest the important role of genome position for the origin of new genes.

As it is the case for non-coding regions, the coding regions of organisms can be expanded or shrunk during evolutionary processes. However, the dynamics of coding regions are expected to be more correlated with functional complexity and diversity than are the dynamics of non-coding regions. Hence, it is interesting to investigate the increase of diversity in coding regions--the origin and evolution of new genes - because this provides a new component to the genetic variation underlying the diversity of living organisms. Here, we examine what is known about the mechanisms responsible for the increase in gene number. Every mechanism affects genomes in a distinct way and to a different extent and it appears that certain organisms favor particular mechanisms. The detail of some interesting gene acquisitions reveals the extreme dynamism of genomes. Finally, we discuss what is known about the fate of new genes and conclude that many of the acquisitions are likely to have been driven by natural selection; they increase functional complexity, diversity, and/or adaptation of species. Despite this, the correlation between complexity of life and gene number is low and closely related species (with very similar life histories) can have very different number of genes. We call this phenomenon the G-value paradox.

Processed genes are created by retroposition from messenger RNA of expressed genes. The estimated amount of processed copies of genes in the human genome is 10,000-14,000. Some of these might be pseudogenes with the expected pattern for nonfunctional sequences, but some others might be an important source of new genes. We have studied the evolution of a Phosphoglycerate mutase processed gene (PGAM3) described in humans and believed to be a pseudogene. We sequenced PGAM3 in chimpanzee and macaque and obtained polymorphism data for human coding region. We found evidence that PGAM3 likely produces a functional protein, as an example of addressing functionality for human processed pseudogenes. First, the open reading frame was intact despite many deletions that occurred in the 3' untranslated region. Second, it appears that the gene is expressed. Finally, interspecies and intraspecies variation for PGAM3 was not consistent with the neutral model proposed for pseudogenes, suggesting that a new functional primate gene has originated. Amino acid divergence was significantly higher than synonymous divergence in PGAM3 lineage, supporting positive selection acting in this gene. This role of selection was further supported by the excess of rare alleles in a population genetic analysis. PGAM3 is located in a region of very low recombination; therefore, it is conceivable that the rapid fixation events in this newly arising gene may have contributed to a selective sweep of variation in the region.

It has recently been discovered that the Alcohol dehydrogenase and Alcohol dehydrogenase-related genes of Drosophila melanogaster and closely related species constitute a single transcription unit and that the Alcohol dehydrogenase-related gene is exclusively expressed from a dicistronic mRNA. Here, we show that in Drosophila lebanonensis, subgenus Scaptodrosophila, Adhr: is also transcribed as a dicistronic transcript with Adh Using degenerate primers designed on the sequence of the known Adhr proteins, we have been able to amplify and sequence a partial sequence of Adhr: in species representative of the whole subgenus Drosophila. This has allowed the study of the organization and expression of Adhr: in Drosophila buzzatii. We find that in D. buzzatii Adhr is transcribed as a monocistronic transcript. Adh and Adhr are believed to originate by duplication, and our data suggest that the cotranscription of these two genes was the primitive state, and that their independent transcription in the subgenus Drosophila is derived. We can rationalize the D. buzzatii condition as being correlated with the two genes evolving independent transcriptional control. However, why these two genes with clear divergence in the functions of their proteins should remain cotranscribed in groups as divergent as the subgenus Sophophora and the subgenus Scaptodrosophila remains a mystery.

The molecular organization of a 1.944-Mb chromosomal region of Drosophila melanogaster around the Adh locus has been analyzed in two repleta group species: D. repleta and D. buzzatii. The extensive genetic and molecular information about this region in D. melanogaster makes it a prime choice for comparative studies of genomic organization among distantly related species. A set of 26 P1 phages from D. melanogaster were successfully hybridized using fluorescence in-situ hybridization (FISH) to the salivary gland chromosomes of both repleta group species. The results show that the Adh region is distributed in D. repleta and D. buzatii over six distant sites of chromosome 3, homologous to chromosomal arm 2L of D. melanogaster (Muller's element B). This observation implies a density of 2.57 fixed breakpoints per Mb in the Adh region and suggests a considerable reorganization of this chromosomal element via the fixation of paracentric inversions. Nevertheless, breakpoint density in the Adh region is three times lower than that estimated for D. repleta chromosome 2, homologous to D. melanogaster 3R (Muller's element E). Differences in the rate of evolution among chromosomal elements are seemingly persistent in the Drosophila genus over long phylogenetic distances.

A simple way to think of evolutionary trade-offs is to suppose genetic effects of opposed direction that give rise to antagonistic pleiotropy. Maintenance of additive genetic variability for fitness related characters, in association with negative correlations between these characters, may result. In the cactophilic species Drosophila buzzatii, there is evidence that second-chromosome polymorphic inversions affect size-related traits. Because a trade-off between body size and larval developmental time has been reported in Drosophila, we study here whether or not these inversions also affect larva-adult viability and developmental time. In particular, we expect that polymorphic inversions make a statistically significant contribution to the genetic correlation between body size (as measured by thorax length) and larval developmental time. This contribution is expected to be in the direction predicted by the trade-off, namely, those flies whose karyotypes cause them to be genetically larger should also have a longer developmental time than flies with other karyotypes. Using two different experimental approaches, a statistically significant contribution of the second-chromosome inversions to the phenotypic variances of body size and developmental time in D. buzzatii was found. Further, these inversions make a positive contribution to the total genetic correlation between the traits, as expected by the suggested trade-off. The data do not provide evidence as to whether the genetic correlation is due to antagonistic pleiotropic gene action or to gametic disequilibrium of linked genes that affect one or both traits. The results do suggest, however, a possible explanation for the maintenance of inversion polymorphism in this species., (© 1998 The Society for the Study of Evolution.)

A theoretical analysis of the effects of inversions on recombination and gene flux between arrangements caused by gene conversion and crossing over was carried out. Two different mathematical models of recombination were used: the Poisson model (without interference) and the Counting model (with interference). The main results are as follows. (1) Recombination and gene flux are highly site-dependent both inside and outside the inverted regions. (2) Crossing over overwhelms gene conversion as a cause of gene flux in large inversions, while conversion becomes relatively significant in short inversions and in regions around the breakpoints. (3) Under the Counting model the recombination rate between two markers depends strongly on the position of the markers along the inverted segment. Two equally spaced markers in the central part of the inverted segment have less recombination than if they are in a more extreme position. (4) Inversions affect recombination rates in the univerted regions of the chromosome. Recombination increases in the distal segment and decreases in the proximal segment. These results provide an explanation for a number of observations reported in the literature. Because inversions are ubiquitous in the evolutionary history of many Drosophila species, the effects of inversions on recombination are expected to influence DNA variation patterns.

DNA sequence variation studies report the transfer of small segments of DNA among different sequences caused by gene conversion events. Here, we provide an algorithm to detect gene conversion tracts and a statistical model to estimate the number and the length distribution of conversion tracts for population DNA sequence data. Two length distributions are defined in the model: (1) that of the observed tract lengths and (2) that of the true tract lengths. If the latter follows a geometric distribution, the relationship between both distributions depends on two basic parameters: psi, which measures the probability of detecting a converted site, and phi, the parameter of the geometric distribution, from which the average true tract length, 1/(1-phi), can be estimated. Expressions are provided for estimating phi by the method of the moments and that of the maximum likelihood. The robustness of the model is examined by computer simulation. The present methods have been applied to the published rp49 sequences of Drosophila subobscura. Maximum likelihood estimate of phi for this data set is 0.9918, which represents an average conversion tract length of 122 bp. Only a small percentage of extant conversion events is detected.

The total gametic disequilibrium between two loci linked to polymorphic inversions can be partitioned into two types of components: within and between chromosome arrangements. The within components depend on the gametic disequilibrium within each chromosome arrangement. The between components depend on the locus-inversion disequilibria. This partitioning has practical applications and is indispensable for studying the dynamics of these systems because inversions greatly reduce recombination in the heterokaryotypes while allowing free, and sometimes different, recombination in each of the homokaryotypes. We provide equations for the per generation change of the various disequilibria for systems with two and three chromosome arrangements, and the general recursive equations predicting the disequilibria after any number of generations for the case of two arrangements. Simulation studies were carried out using different values of the recombination parameters and all possible initial conditions. The results show a complex convergence to linkage equilibrium in inversion systems. The various disequilibria can have local maxima and minima while approaching equilibrium and, moreover, their dynamics cannot be described, in general, using a single parameter, i.e. an effective recombination rate. We conclude that the effects of inversions on gametic disequilibria must be carefully considered when dealing with disequilibria in inversion systems. The formulae provided in this paper can be used for such purpose.