Emergence of young human genes after a burst of retroposition in primates

The origin of new genes through gene duplication is fundamental to the evolution of lineage- or species-specific phenotypic traits. In this report, we estimate the number of functional retrogenes on the lineage leading to humans generated by the high rate of retroposition (retroduplication) in primates. Extensive comparative sequencing and expression studies coupled with evolutionary analyses and simulations suggest that a significant proportion of recent retrocopies represent bona fide human genes. We estimate that at least one new retrogene per million years emerged on the human lineage during the past ∼63 million years of primate evolution. Detailed analysis of a subset of the data shows that the majority of retrogenes are specifically expressed in testis, whereas their parental genes show broad expression patterns. Consistently, most retrogenes evolved functional roles in spermatogenesis. Proteins encoded by X chromosome−derived retrogenes were strongly preserved by purifying selection following the duplication event, supporting the view that they may act as functional autosomal substitutes during X-inactivation of late spermatogenesis genes. Also, some retrogenes acquired a new or more adapted function driven by positive selection. We conclude that retroduplication significantly contributed to the formation of recent human genes and that most new retrogenes were progressively recruited during primate evolution by natural and/or sexual selection to enhance male germline function.     

Analysis of new retrogenes provides insight into dog adaptive evolution

The origin and subsequent evolution of new genes have been considered as an important source of genetic and phenotypic diversity in organisms. Dog breeds show great phenotypic diversity for morphological, physiological, and behavioral traits. However, the contributions of newly originated retrogenes, which provide important genetic bases for dog species differentiation and adaptive traits, are largely unknown. Here, we analyzed the dog genome to identify new RNA‐based duplications and comprehensively investigated their origin, evolution, functions in adaptive traits, and gene movement processes. First, we totally identified 3,025 retrocopies including 476 intact retrogenes, 2,518 retropseudogenes, and 31 chimerical retrogenes. Second, selective pressure along with ESTs expression analysis showed that most of the intact retrogenes were significantly under stronger purifying selection and subjected to more functional constraints when compared to retropseudogenes. Furthermore, a large number of retrocopies and chimerical retrogenes that occurred approximately 22 million years ago implied a burst of retrotransposition in the dog genome after the divergence time between dog and its closely related species red fox. Interestingly, GO and pathway analyses showed that new retrogenes had expanded in glutathione biosynthetic/metabolic process which likely provided important genetic basis for dogs' adaptation to scavenge human waste dumps. Finally, consistent with the results in human and mouse, a significant excess of functional retrogenes movement on and off the X chromosome in the dog confirmed a general pattern of gene movement process in mammals which was likely driven by natural selection or sexual antagonism. Together, these results increase our understanding that new retrogenes can reshape the dog genome and provide further exploration of the molecular mechanisms underlying the dogs' adaptive evolution.     

Evolutionary patterns of RNA-based duplication in non-mammalian chordates

The role of RNA-based duplication, or retroposition, in the evolution of new gene functions in mammals, plants, and Drosophila has been widely reported. However, little is known about RNA-based duplication in non-mammalian chordates. In this study, we screened ten non-mammalian chordate genomes for retrocopies and investigated their evolutionary patterns. We identified numerous retrocopies in these species. Examination of the age distribution of these retrocopies revealed no burst of young retrocopies in ancient chordate species. Upon comparing these non-mammalian chordate species to the mammalian species, we observed that a larger fraction of the non-mammalian retrocopies was under strong evolutionary constraints than mammalian retrocopies are, as evidenced by signals of purifying selection and expression profiles. For the Western clawed frog, Medaka, and Sea squirt, many retrogenes have evolved gonad and brain expression patterns, similar to what was observed in human. Testing of retrogene movement in the Medaka genome, where the nascent sex chrosomes have been well assembled, did not reveal any significant gene movement. Taken together, our analyses demonstrate that RNA-based duplication generates many functional genes and can make a significant contribution to the evolution of non-mammalian genomes.     

Inter-population Differences in Retrogene Loss and Expression in Humans

Gene retroposition leads to considerable genetic variation between individuals. Recent studies revealed the presence of at least 208 retroduplication variations (RDVs), a class of polymorphisms, in which a retrocopy is present or absent from individual genomes. Most of these RDVs resulted from recent retroduplications. In this study, we used the results of Phase 1 from the 1000 Genomes Project to investigate the variation in loss of ancestral (i.e. shared with other primates) retrocopies among different human populations. In addition, we examined retrocopy expression levels using RNA-Seq data derived from the Ilumina BodyMap project, as well as data from lymphoblastoid cell lines provided by the Geuvadis Consortium. We also developed a new approach to detect novel retrocopies absent from the reference human genome. We experimentally confirmed the existence of the detected retrocopies and determined their presence or absence in the human genomes of 17 different populations. Altogether, we were able to detect 193 RDVs; the majority resulted from retrocopy deletion. Most of these RDVs had not been previously reported. We experimentally confirmed the expression of 11 ancestral retrogenes that underwent deletion in certain individuals. The frequency of their deletion, with the exception of one retrogene, is very low. The expression, conservation and low rate of deletion of the remaining 10 retrocopies may suggest some functionality. Aside from the presence or absence of expressed retrocopies, we also searched for differences in retrocopy expression levels between populations, finding 9 retrogenes that undergo statistically significant differential expression.     

Origins of new male germ-line functions from X-derived autosomal retrogenes in the mouse

Recent literature demonstrates that retrogenes tend to leave the X chromosome and integrate onto the autosomes and evolve male-biased expression patterns. Several selection-based evolutionary mechanisms have been proposed to explain this observation. Testing these selection-based models requires examining the evolutionary history and functional properties of new retrogenes, particularly those that show evidence of directional movement between the X and the autosomes (X-related retrogenes). This includes autosomal retrogenes with parental paralogs on the X chromosome (X-derived autosomal retrogenes) and those retrogenes integrated onto the X chromosomes (X-linked retrogenes). In order to understand why retrogenes tend to move nonrandomly in genomes, we examined the expression patterns and evolutionary mechanisms concerning gene pairs having young retrogenes--originating less than 20 MYA (after mouse-rat split). We demonstrate that these X-derived autosomal retrogenes evolved a more restricted male-biased expression pattern: they are expressed exclusively or predominantly in the testis, in particular, during the late stages of spermatogenesis. In contrast, the parental counterparts have relatively broad expression patterns in various tissues and spermatogenetic stages. We further observed that positive selection is targeting these X-derived autosomal retrogenes with novel male-biased expression patterns. This suggests that such retrogenes evolved new male germ-line functions that may be complementary to the functions of the parental paralogs, which themselves contribute little during spermatogenesis. Such evolutionary changes may be beneficial to the populations. Furthermore, most identified X-related retrogenes have recruited novel adjacent sequences as their untranslated regions (UTRs), suggesting that these UTRs, acquired de novo, may play an important role in establishing new regulatory mechanisms to carry out the new male germ-line functions.     

Evolution of Haplotypes at the DRD2 Locus

We present here the first evolutionary perspective on haplotypes at DRD2, the locus for the dopamine D₂ receptor. The dopamine D₂ receptor plays a critical role in the functioning of many neural circuits in the human brain. If functionally relevant variation at the DRD2 locus exists, understanding the evolution of haplotypes on the basis of polymorphic sites encompassing the gene should provide a powerful framework for identifying that variation. Three DRD2 polymorphisms (TaqI “A” and “B” RFLPs and the (CA)_n short tandem repeat polymorphism) encompassing the coding sequences have been studied in 15 populations; these markers are polymorphic in all the populations studied, and they display strong and significant linkage disequilibria with each other. The common haplotypes for the two TaqI RFLPs are separately derived from the ancestral haplotype but predate the spread of modern humans around the world. The knowledge of how the various haplotypes have evolved, the allele frequencies of the haplotypes in human populations, and the physical relationships of the polymorphisms to each other and to the functional parts of the gene should now allow proper design and interpretation of association studies.     

A case-by-case evolutionary analysis of four imprinted retrogenes

Retroposition is a widespread phenomenon resulting in the generation of new genes that are initially related to a parent gene via very high coding sequence similarity. We examine the evolutionary fate of four retrogenes generated by such an event; mouse Inpp5f_v2, Mcts2, Nap1l5, and U2af1-rs1. These genes are all subject to the epigenetic phenomenon of parental imprinting. We first provide new data on the age of these retrogene insertions. Using codon-based models of sequence evolution, we show these retrogenes have diverse evolutionary trajectories, including divergence from the parent coding sequence under positive selection pressure, purifying selection pressure maintaining parent-retrogene similarity, and neutral evolution. Examination of the expression pattern of retrogenes shows an atypical, broad pattern across multiple tissues. Protein 3D structure modeling reveals that a positively selected residue in U2af1-rs1, not shared by its parent, may influence protein conformation. Our case-by-case analysis of the evolution of four imprinted retrogenes reveals that this interesting class of imprinted genes, while similar in regulation and sequence characteristics, follow very varied evolutionary paths.     

Evolutionary patterns of RNA-based gene duplicates in Caenorhabditis nematodes coincide with their genomic features

ABSTRACT: BACKGROUND: RNA-based gene duplicates (retrocopies) played pivotal roles in many physiological processes. Nowadays, functional retrocopies have been systematically identified in several mammals, fruit flies, plants, zebrafish and other chordates, etc. However, studies about this kind of duplication in Caenorhabditis nematodes have not been reported. FINDINGS: We identified 43, 48, 43, 9, and 42 retrocopies, of which 6, 15, 18, 3, and 13 formed chimeric genes in C. brenneri, C. briggsae, C. elegans, C. japonica, and C. remanei, respectively. At least 5 chimeric types exist in Caenorhabditis species, of which retrocopy recruiting both N and C terminus is the commonest one. Evidences from different analyses demonstrate many retrocopies and almost all chimeric genes may be functional in these species. About half of retrocopies in each species has coordinates in other species, and we suggest that retrocopies in closely related species may be helpful in identifying retrocopies for one certain species. CONCLUSIONS: A number of retrocopies and chimeric genes exist in Caenorhabditis genomes, and some of them may be functional. The evolutionary patterns of these genes may correlate with their genomic features, such as the activity of retroelements, the high rate of mutation and deletion rate, and a large proportion of genes subject to trans-splicing.     

Rapid genome‐wide evolution in Brassica rapa populations following drought revealed by sequencing of ancestral and descendant gene pools

There is increasing evidence that evolution can occur rapidly in response to selection. Recent advances in sequencing suggest the possibility of documenting genetic changes as they occur in populations, thus uncovering the genetic basis of evolution, particularly if samples are available from both before and after selection. Here, we had a unique opportunity to directly assess genetic changes in natural populations following an evolutionary response to a fluctuation in climate. We analysed genome‐wide differences between ancestors and descendants of natural populations of Brassica rapa plants from two locations that rapidly evolved changes in multiple phenotypic traits, including flowering time, following a multiyear late‐season drought in California. These ancestor‐descendant comparisons revealed evolutionary shifts in allele frequencies in many genes. Some genes showing evolutionary shifts have functions related to drought stress and flowering time, consistent with an adaptive response to selection. Loci differentiated between ancestors and descendants (F_ST outliers) were generally different from those showing signatures of selection based on site frequency spectrum analysis (Tajima's D), indicating that the loci that evolved in response to the recent drought and those under historical selection were generally distinct. Very few genes showed similar evolutionary responses between two geographically distinct populations, suggesting independent genetic trajectories of evolution yielding parallel phenotypic changes. The results show that selection can result in rapid genome‐wide evolutionary shifts in allele frequencies in natural populations, and highlight the usefulness of combining resurrection experiments in natural populations with genomics for studying the genetic basis of adaptive evolution.     

Loci under selection during multiple range expansions of an invasive plant are mostly population specific,but patterns are associated with climate

Identifying the genes underlying rapid evolutionary changes, describing their function and ascertaining the environmental pressures that determine fitness are the central elements needed for understanding of evolutionary processes and phenotypic changes that improve the fitness of populations. It has been hypothesized that rapid adaptive changes in new environments may contribute to the rapid spread and success of invasive plants and animals. As yet, studies of adaptation during invasion are scarce, as is knowledge of the genes underlying adaptation, especially in multiple replicated invasions. Here, we quantified how genotype frequencies change during invasions, resulting in rapid evolution of naturalized populations. We used six fully replicated common garden experiments in Brazil where Pinus taeda (loblolly pine) was introduced at the same time, in the same numbers, from the same seed sources, and has formed naturalized populations expanding outward from the plantations. We used a combination of nonparametric, population genetics and multivariate statistics to detect changes in genotype frequencies along each of the six naturalization gradients and their association with climate as well as shifts in allele frequencies compared to the source populations. Results show 25 genes with significant shifts in genotype frequencies. Six genes had shifts in more than one population. Climate explained 25% of the variation in the groups of genes under selection across all locations, but specific genes under strong selection during invasions did not show climate‐related convergence. In conclusion, we detected rapid evolutionary changes during invasive range expansions, but the particular gene‐level patterns of evolution may be population specific.     

Tigers of Sundarbans in India: Is the Population a Separate Conservation Unit?

The Sundarbans tiger inhabits a unique mangrove habitat and are morphologically distinct from the recognized tiger subspecies in terms of skull morphometrics and body size. Thus, there is an urgent need to assess their ecological and genetic distinctiveness and determine if Sundarbans tigers should be defined and managed as separate conservation unit. We utilized nine microsatellites and 3 kb from four mitochondrial DNA (mtDNA) genes to estimate genetic variability, population structure, demographic parameters and visualize historic and contemporary connectivity among tiger populations from Sundarbans and mainland India. We also evaluated the traits that determine exchangeability or adaptive differences among tiger populations. Data from both markers suggest that Sundarbans tiger is not a separate tiger subspecies and should be regarded as Bengal tiger (P. t. tigris) subspecies. Maximum likelihood phylogenetic analyses of the mtDNA data revealed reciprocal monophyly. Genetic differentiation was found stronger for mtDNA than nuclear DNA. Microsatellite markers indicated low genetic variation in Sundarbans tigers (He= 0.58) as compared to other mainland populations, such as northern and Peninsular (Hebetween 0.67- 0.70). Molecular data supports migration between mainland and Sundarbans populations until very recent times. We attribute this reduction in gene flow to accelerated fragmentation and habitat alteration in the landscape over the past few centuries. Demographic analyses suggest that Sundarbans tigers have diverged recently from peninsular tiger population within last 2000 years. Sundarbans tigers are the most divergent group of Bengal tigers, and ecologically non-exchangeable with other tiger populations, and thus should be managed as a separate “evolutionarily significant unit” (ESU) following the adaptive evolutionary conservation (AEC) concept.     

Phenotypic plasticity in gene expression contributes to divergence of locally adapted populations of Fundulus heteroclitus

We examine the interaction between phenotypic plasticity and evolutionary adaptation using muscle gene expression levels among populations of the fish Fundulus heteroclitus acclimated to three temperatures. Our analysis reveals shared patterns of phenotypic plasticity due to thermal acclimation as well as non‐neutral patterns of variation among populations adapted to different thermal environments. For the majority of significant differences in gene expression levels, phenotypic plasticity and adaptation operate on different suites of genes. The subset of genes that demonstrate both adaptive differences and phenotypic plasticity, however, exhibit countergradient variation of expression. Thus, expression differences among populations counteract environmental effects, reducing the phenotypic differentiation between populations. Finally, gene‐by‐environment interactions among genes with non‐neutral patterns of expression suggest that the penetrance of adaptive variation depends on the environmental conditions experienced by the individual.     

A spatially aware likelihood test to detect sweeps from haplotype distributions

The inference of positive selection in genomes is a problem of great interest in evolutionary genomics. By identifying putative regions of the genome that contain adaptive mutations, we are able to learn about the biology of organisms and their evolutionary history. Here we introduce a composite likelihood method that identifies recently completed or ongoing positive selection by searching for extreme distortions in the spatial distribution of the haplotype frequency spectrum along the genome relative to the genome-wide expectation taken as neutrality. Furthermore, the method simultaneously infers two parameters of the sweep: the number of sweeping haplotypes and the “width” of the sweep, which is related to the strength and timing of selection. We demonstrate that this method outperforms the leading haplotype-based selection statistics, though strong signals in low-recombination regions merit extra scrutiny. As a positive control, we apply it to two well-studied human populations from the 1000 Genomes Project and examine haplotype frequency spectrum patterns at the LCT and MHC loci. We also apply it to a data set of brown rats sampled in NYC and identify genes related to olfactory perception. To facilitate use of this method, we have implemented it in user-friendly open source software.     

Transposition of reversed Ac element ends generates novel chimeric genes in maize

The maize Activator/Dissociation (Ac/Ds) elements are members of the hAT (hobo, Ac, and Tam3) superfamily of type II (DNA) transposons that transpose through a “cut-and-paste” mechanism. Previously, we reported that a pair of Ac ends in reversed orientation is capable of undergoing alternative transposition reactions that can generate large-scale chromosomal rearrangements, including deletions and inversions. We show here that rearrangements induced by reversed Ac ends transposition can join the coding and regulatory sequences of two linked paralogous genes to generate a series of chimeric genes, some of which are functional. To our knowledge, this is the first report demonstrating that alternative transposition reactions can recombine gene segments, leading to the creation of new genes.