Search

De novo gene birth—the evolution of new protein-coding genes from ancestrally noncoding DNA—is increasingly appreciated as an important source of genetic and phenotypic innovation. However, the frequency and overall biological impact of de novo genes (DNGs) remain controversial. Large-scale surveys of de novo genes are critical to address these issues, but DNG identification represents a persistent challenge due to the lack of standardized protocols and the laborious analyses traditionally used to detect DNGs. Here, we introduced novel approaches to identify de novo genes that rely on Machine Learning Algorithms (MLAs) and are poised to accelerate DNG discovery. We specifically investigated if MLAs developed in one species using known DNGs can accurately predict de novo genes in other genomes. To maximize the applicability of these methods across species, we relied only on DNA and protein sequence features that can be easily obtained from annotation data. Using hundreds of published and newly annotated DNGs from three angiosperms, we trained and tested both Decision Tree (DT) and Neural Network (NN) algorithms. Both MLAs showed high levels of accuracy and recall within-genomes. Although accuracies and recall decreased in cross-species analyses, they remained elevated between evolutionary closely related species. A few training features, including presence of a protein domain and coding probability, held most of the MLAs predictive power. In analyses of all genes from a genome, recall was still elevated. Although false positive rates were relatively high, MLA screenings of whole-genome datasets reduced by up to ten-fold the number of genes to be examined by conventional comparative genomic methods. Thus, a combination of MLAs and traditional strategies can significantly accelerate the accurate discovery of DNG and the annotation in angiosperm genomes.

Background The recurrent evolution of the C4 photosynthetic pathway in angiosperms represents one of the most extraordinary examples of convergent evolution of a complex trait. Comparative genomic analyses have unveiled some of the molecular changes associated with the C4 pathway. For instance, several key enzymes involved in the transition from C3 to C4 photosynthesis have been found to share convergent amino acid replacements along C4 lineages. However, the extent of convergent replacements potentially associated with the emergence of C4 plants remains to be fully assessed. Here, we conducted an organelle-wide analysis to determine if convergent evolution occurred in multiple chloroplast proteins beside the well-known case of the large RuBisCO subunit encoded by the chloroplast gene rbcL. Methods Our study was based on the comparative analysis of 43 C4 and 21 C3 grass species belonging to the PACMAD clade, a focal taxonomic group in many investigations of C4 evolution. We first used protein sequences of 67 orthologous chloroplast genes to build an accurate phylogeny of these species. Then, we inferred amino acid replacements along 13 C4 lineages and 9 C3 lineages using reconstructed protein sequences of their reference branches, corresponding to the branches containing the most recent common ancestors of C4-only clades and C3-only clades. Pairwise comparisons between reference branches allowed us to identify both convergent and non-convergent amino acid replacements between C4:C4, C3:C3 and C3:C4 lineages. Results The reconstructed phylogenetic tree of 64 PACMAD grasses was characterized by strong supports in all nodes used for analyses of convergence. We identified 217 convergent replacements and 201 non-convergent replacements in 45/67 chloroplast proteins in both C4 and C3 reference branches. C4:C4 branches showed higher levels of convergent replacements than C3:C3 and C3:C4 branches. Furthermore, we found that more proteins shared unique convergent replacements in C4 lineages, with both RbcL and RpoC1 (the RNA polymerase beta’ subunit 1) showing a significantly higher convergent/non-convergent replacements ratio in C4 branches. Notably, more C4:C4 reference branches showed higher numbers of convergent vs. non-convergent replacements than C3:C3 and C3:C4 branches. Our results suggest that, in the PACMAD clade, C4 grasses experienced higher levels of molecular convergence than C3 species across multiple chloroplast genes. These findings have important implications for our understanding of the evolution of the C4 photosynthesis pathway.

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 35th biennial Southern Forest Tree Improvement Conference (SFTIC) was held on the 4th and 5th day of June 2019, in Lexington, KY, USA. The theme of the conference was “Genetics and improvement of forest health and productivity”. Plenary speakers detailed significant progress in understanding forest tree genomes, identifying candidate genes for disease resistance and tree architecture traits and utilizing both genomic selection and genetic modification to advance forest health and productivity. Concurrent and general technical sessions provided additional evidence of progress in these areas as well as information on seed zone development for the eastern USA and species prioritization for conservation and restoration programs, in addition to advanced genome-assisted selection and breeding methods and programs in species of commercial interest including pines and eucalyptus. An interactive poster session provided each poster presenter an opportunity to speak in front of a general session audience and to take Q&A in a traditional poster session format. Pre-conference workshops and a post-conference field tour rounded out the 35th SFTIC with 88 registered for the conference, more than 60 participating in workshops and 33 attending the field tour. The conference program and proceedings are accessible through the SFTIC web site at http://www.sftic.org.

Copy number variants are insertions and deletions of 1 kb or larger in a genome that play an important role in phenotypic changes and human disease. Many software applications have been developed to detect copy number variants using either whole-genome sequencing or whole-exome sequencing data. However, there is poor agreement in the results from these applications. Simulated datasets containing copy number variants allow comprehensive comparisons of the operating characteristics of existing and novel copy number variant detection methods. Several software applications have been developed to simulate copy number variants and other structural variants in whole-genome sequencing data. However, none of the applications reliably simulate copy number variants in whole-exome sequencing data. We have developed and tested SECNVs (Simulator of Exome Copy Number Variants), a fast, robust and customizable software application for simulating copy number variants and whole-exome sequences from a reference genome. SECNVs is easy to install, implements a wide range of commands to customize simulations, can output multiple samples at once, and incorporates a pipeline to output rearranged genomes, short reads and BAM files in a single command. Variants generated by SECNVs are detected with high sensitivity and precision by tools commonly used to detect copy number variants. SECNVs is publicly available at https://github.com/YJulyXing/SECNVs.

Comparative genomics analyses empowered by the wealth of sequenced genomes have revealed numerous instances of horizontal DNA transfers between distantly related species. In eukaryotes, repetitive DNA sequences known as transposable elements (TEs) are especially prone to move across species boundaries. Such horizontal transposon transfers, or HTTs, are relatively common within major eukaryotic kingdoms, including animals, plants, and fungi, while rarely occurring across these kingdoms. Here, we describe the first case of HTT from animals to plants, involving TEs known as Penelope-like elements, or PLEs, a group of retrotransposons closely related to eukaryotic telomerases. Using a combination of in situ hybridization on chromosomes, polymerase chain reaction experiments, and computational analyses we show that the predominant PLE lineage, EN(+)PLEs, is highly diversified in loblolly pine and other conifers, but appears to be absent in other gymnosperms. Phylogenetic analyses of both protein and DNA sequences reveal that conifers EN(+)PLEs, or Dryads, form a monophyletic group clustering within a clade of primarily arthropod elements. Additionally, no EN(+)PLEs were detected in 1,928 genome assemblies from 1,029 nonmetazoan and nonconifer genomes from 14 major eukaryotic lineages. These findings indicate that Dryads emerged following an ancient horizontal transfer of EN(+)PLEs from arthropods to a common ancestor of conifers approximately 340 Ma. This represents one of the oldest known interspecific transmissions of TEs, and the most conspicuous case of DNA transfer between animals and plants.

The evolution of novel protein-coding genes from noncoding regions of the genome is one of the most compelling pieces of evidence for genetic innovations in nature. One popular approach to identify de novo genes is phylostratigraphy, which consists of determining the approximate time of origin (age) of a gene based on its distribution along a species phylogeny. Several studies have revealed significant flaws in determining the age of genes, including de novo genes, using phylostratigraphy alone. However, the rate of false positives in de novo gene surveys, based on phylostratigraphy, remains unknown. Here, I reanalyze the findings from three studies, two of which identified tens to hundreds of rodent-specific de novo genes adopting a phylostratigraphy-centered approach. Most putative de novo genes discovered in these investigations are no longer included in recently updated mouse gene sets. Using a combination of synteny information and sequence similarity searches, I show that ∼60% of the remaining 381 putative de novo genes share homology with genes from other vertebrates, originated through gene duplication, and/or share no synteny information with nonrodent mammals. These results led to an estimated rate of ∼12 de novo genes per million years in mouse. Contrary to a previous study (Wilson BA, Foy SG, Neme R, Masel J. 2017. Young genes are highly disordered as predicted by the preadaptation hypothesis of de novo gene birth. Nat Ecol Evol. 1:0146), I found no evidence supporting the preadaptation hypothesis of de novo gene formation. Nearly half of the de novo genes confirmed in this study are within older genes, indicating that co-option of preexisting regulatory regions and a higher GC content may facilitate the origin of novel genes.

The evolution of novel protein-coding genes from noncoding regions of the genome is one of the most compelling evidence for genetic innovations in nature. One popular approach to identify de novo genes is phylostratigraphy, which consists of determining the approximate time of origin (age) of a gene based on its distribution along a species phylogeny. Several studies have revealed significant flaws in determining the age of genes, including de novo genes, using phylostratigraphy alone. However, the rate of false positives in de novo gene surveys, based on phylostratigraphy, remains unknown. Here, I re-analyze the findings from three studies, two of which identified tens to hundreds of rodent-specific de novo genes adopting a phylostratigraphy-centered approach. Most of the putative de novo genes discovered in these investigations are no longer included in recently updated mouse gene sets. Using a combination of synteny information and sequence similarity searches, I show that about 60% of the remaining 381 putative de novo genes share homology with genes from other vertebrates, originated through gene duplication, and/or share no synteny information with non-rodent mammals. These results led to an estimated rate of ∼12 de novo genes per million year in mouse. Contrary to a previous study (Wilson et al. 2017), I found no evidence supporting the preadaptation hypothesis of de novo gene formation. Nearly half of the de novo genes confirmed in this study are within older genes, indicating that co-option of preexisting regulatory regions and a higher GC content may facilitate the origin of novel genes.

Gene duplications and gene losses are major determinants of genome evolution and phenotypic diversity. The frequency of gene turnover (gene gains and gene losses combined) is known to vary between organisms. Comparative genomic analyses of gene families can highlight such variation; however, estimates of gene turnover may be biased when using highly fragmented genome assemblies resulting in poor gene annotations. Here, we address potential biases introduced by gene annotation errors in estimates of gene turnover frequencies in a dataset including both well-annotated angiosperm genomes and the incomplete gene sets of four Pinaceae, including two pine species, Norway spruce and Douglas-fir. We show that Pinaceae experienced higher gene turnover rates than angiosperm lineages lacking recent whole-genome duplications. This finding is robust to both known major issues in Pinaceae gene sets: missing gene models and erroneous annotation of pseudogenes. A separate analysis limited to the four Pinaceae gene sets pointed to an accelerated gene turnover rate in pines compared with Norway spruce and Douglas-fir. Our results indicate that gene turnover significantly contributes to genome variation and possibly to speciation in Pinaceae, particularly in pines. Moreover, these findings indicate that reliable estimates of gene turnover frequencies can be discerned in incomplete and potentially inaccurate gene sets. Because gymnosperms are known to exhibit low overall substitution rates compared with angiosperms, our results suggest that the rate of single-base pair mutations is uncoupled from the rate of large DNA duplications and deletions associated with gene turnover in Pinaceae.

Gene duplications and gene losses are major determinants of genome evolution and phenotypic diversity. The frequency of gene turnover (gene gains and gene losses combined) is known to vary between organisms. Comparative genomic analyses of gene families can highlight such variation; however, estimates of gene turnover rates may be biased when using highly fragmented genome assemblies resulting in poor gene annotations. Here, we address potential biases introduced by gene annotation errors in estimates of gene turnover frequencies in a dataset including both well-annotated angiosperm genomes and the incomplete gene sets of four Pinaceae including two pine species, Norway spruce and Douglas-fir. Previous studies have shown low overall substitution rates, but higher levels of adaptive substitutions, in genes from Pinaceae and other gymnosperms compared to angiosperms. Conversely, our analysis suggests that pines experienced higher gene turnover rates than angiosperm lineages lacking recent whole-genome duplications. This finding is robust to both known major issues in Pinaceae gene sets: missing gene models and erroneous annotation of pseudogenes. A separate analysis limited to the four Pinaceae gene sets confirmed an accelerated gene turnover rate in pines compared to Norway spruce and Douglas-fir. Our results indicate that gene turnover significantly contributes to genome variation and possibly to adaptation and speciation in Pinaceae. Moreover, these findings indicate that reliable estimates of gene turnover frequencies can be discerned in incomplete and potentially inaccurate gene sets.

The Marmoset Genome Sequencing and Analysis Consortium.-- Worley, Kim C. et al., We report the whole-genome sequence of the common marmoset (Callithrix jacchus). The 2.26-Gb genome of a female marmoset was assembled using Sanger read data (6×) and a whole-genome shotgun strategy. A first analysis has permitted comparison with the genomes of apes and Old World monkeys and the identification of specific features that might contribute to the unique biology of this diminutive primate, including genetic changes that may influence body size, frequent twinning and chimerism. We observed positive selection in growth hormone/insulin-like growth factor genes (growth pathways), respiratory complex I genes (metabolic pathways), and genes encoding immunobiological factors and proteases (reproductive and immunity pathways). In addition, both protein-coding and microRNA genes related to reproduction exhibited evidence of rapid sequence evolution. This genome sequence for a New World monkey enables increased power for comparative analyses among available primate genomes and facilitates biomedical research application. © 2014 Nature America, Inc., The marmoset genome project was funded by the National Human Genome Research Institute (NHGRI), including from grants U54 HG003273 (R.A. Gibbs) and U54 HG003079 (R.K.W.), with additional support from the US National Institutes of Health (NIH), including from grants R01 DK077639 (S.D.T.), R01 GM59290 (L.B.J. and M.A.B.), HG002385 (E.E.E.) and P51-OD011133 (Southwest NPRC), and support from the National Science Foundation (NSF BCS-0751508 to D.E.W.) and the VEGA grant agency: 1/0719/14 (T.V.) and 1/1085/12 (B.B.). C.C.F. and M.C.R. were supported in part by a Howard Hughes Medical Institute grant to Louisiana State University through the Undergraduate Biological Sciences Education program. J.X. was supported by NHGRI grant K99 HG005846. P.H.G. was supported by the Cullen Foundation. T.M.-B. was supported by European Research Council Starting Grant (260372) and MICINN (Spain) grant BFU2011-28549. B.L.-G. was supported by the Spanish National Institute of Bioinformatics (see URLs).

The accumulation and removal of transposable elements (TEs) is a major driver of genome size evolution in eukaryotes. In plants, long terminal repeat (LTR) retrotransposons (LTR-RTs) represent the majority of TEs and form most of the nuclear DNA in large genomes. Unequal recombination (UR) between LTRs leads to removal of intervening sequence and formation of solo-LTRs. UR is a major mechanism of LTR-RT removal in many angiosperms, but our understanding of LTR-RT-associated recombination within the large, LTR-RT-rich genomes of conifers is quite limited. We employ a novel read-based methodology to estimate the relative rates of LTR-RT-associated UR within the genomes of four conifer and seven angiosperm species. We found the lowest rates of UR in the largest genomes studied, conifers and the angiosperm maize. Recombination may also resolve as gene conversion, which does not remove sequence, so we analyzed LTR-RT-associated gene conversion events (GCEs) in Norway spruce and six angiosperms. Opposite the trend for UR, we found the highest rates of GCEs in Norway spruce and maize. Unlike previous work in angiosperms, we found no evidence that rates of UR correlate with retroelement structural features in the conifers, suggesting that another process is suppressing UR in these species. Recent results from diverse eukaryotes indicate that heterochromatin affects the resolution of recombination, by favoring gene conversion over crossing-over, similar to our observation of opposed rates of UR and GCEs. Control of LTR-RT proliferation via formation of heterochromatin would be a likely step toward large genomes in eukaryotes carrying high LTR-RT content.

Establishing the molecular basis of DNA mutations that cause inherited disease is of fundamental importance to understanding the origin, nature, and clinical sequelae of genetic disorders in humans. The majority of disease-associated mutations constitute single-base substitutions and short deletions and/or insertions resulting from DNA replication errors and the repair of damaged bases. However, pathological mutations can also be introduced by nonreciprocal recombination events between paralogous sequences, a phenomenon known as interlocus gene conversion (IGC). IGC events have thus far been linked to pathology in more than 20 human genes. However, the large number of duplicated gene sequences in the human genome implies that many more disease-associated mutations could originate via IGC. Here, we have used a genome-wide computational approach to identify disease-associated mutations derived from IGC events. Our approach revealed hundreds of known pathological mutations that could have been caused by IGC. Further, we identified several dozen high-confidence cases of inherited disease mutations resulting from IGC in ∼1% of all genes analyzed. About half of the donor sequences associated with such mutations are functional paralogous genes, suggesting that epistatic interactions or differential expression patterns will determine the impact upon fitness of specific substitutions between duplicated genes. In addition, we identified thousands of hitherto undescribed and potentially deleterious mutations that could arise via IGC. Our findings reveal the extent of the impact of interlocus gene conversion upon the spectrum of human inherited disease.

Prolactin (PRL) is a multifunctional signaling molecule best known for its role in regulating lactation in mammals. Systemic PRL is produced by the anterior pituitary, but extrapituitary PRL has also been detected in many tissues including the human endometrium. Prolactin is essential for pregnancy in rodents and one of the most dramatically induced genes in the endometrium during human pregnancy. The promoter for human endometrial Prl is located about 5.8 kb upstream of the pituitary promoter and is derived from a transposable element called MER39. Although it has been shown that prolactin is expressed in the pregnant endometrium of a few mammals other than humans, MER39 has been described as primate specific. Thus, in an effort to understand mechanisms of prolactin regulatory evolution, we sought to determine how uterine prolactin is transcribed in species that lack MER39. Using a variety of complementary strategies, including reverse transcriptase-polymerase chain reaction, 5' rapid amplification of cDNA ends, and whole-transcriptome sequencing, we show that endometrial Prl expression is not a shared character of all placental mammals, as it is not expressed in rabbits, pigs, dogs, or armadillos. We show that in primates, mice, and elephants, prolactin mRNA is transcribed in the pregnant endometrium from alternative promoters, different from the pituitary promoter and different from each other. Moreover, we demonstrate that the spider monkey promoter derives from the long terminal repeat (LTR) element MER39 as in humans, the mouse promoter derives from the LTR element MER77, and the elephant promoter derives from the lineage-specific LINE retrotransposon L1-2_LA. We also find surprising variation of transcriptional start sites within these transposable elements and of Prl splice variants, suggesting a high degree of flexibility in the promoter architecture even among closely related species. Finally, the three groups shown here to express endometrial prolactin-the higher primates, the rodents, and the elephant-represent three of the four lineages that showed adaptive evolution of the Prl gene in an earlier study (Wallis M. 2000. Episodic evolution of protein hormones: molecular evolution of pituitary prolactin. J Mol Evol. 50:465-473), which supports our findings and suggests that the selective forces responsible for accelerated Prl evolution were in the endometrium. This is the first reported case of convergent evolution of gene expression through the independent recruitment of different transposable elements, highlighting the importance of transposable elements in gene regulatory, and potentially adaptive, evolution.

Gene conversion between duplicated genes has been implicated in homogenization of gene families and reassortment of variation among paralogs. If conversion is common, this process could lead to errors in gene tree inference and subsequent overestimation of rates of gene duplication. After performing simulations to assess our power to detect gene conversion events, we determined rates of conversion among young, lineage-specific gene duplicates in four mammal species: human, rhesus macaque, mouse, and rat. Gene conversion rates (number of conversion events/number of gene pairs) among young duplicates range from 8.3% in macaque to 18.96% in rat, including a 5% false-positive rate. For all lineages, only 1–3% of the total amount of sequence examined was converted. There is no increase in GC content in conversion tracts compared to flanking regions of the same genes nor in conversion tracts compared to the same region in nonconverted gene-family members, suggesting that ectopic gene conversion does not significantly alter nucleotide composition in these duplicates. While the majority of gene duplicate pairs reside on different chromosomes in mammalian genomes, the majority of gene conversion events occur between duplicates on the same chromosome, even after controlling for divergence between duplicates. Among intrachromosomal duplicates, however, there is no correlation between the probability of conversion and physical distance between duplicates after controlling for divergence. Finally, we use a novel method to show that at most 5–10% of all gene trees involving young duplicates are likely to be incorrect due to gene conversion. We conclude that gene conversion has had only a small effect on mammalian genomes and gene duplicate evolution in general.

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.

Y‐box proteins are a highly conserved family of gene expression regulatory factors. During gametogenesis they may play a dual role, as both transcriptional activators of germ cell‐specific genes and as translational repressors of stored maternal transcripts. We report the identification of RlYB2, a Y‐box homolog gene specifically expressed in the germ cells of green frogs belonging to Rana esculenta complex, a model system characterized by a hybridogenetic gametogenesis. In developing germ cells of the hybrid R. esculenta, arisen by natural cross of the parental species R. ridibunda and R. lessonae, one set of the parental genomes is excluded and the remaining one, first endoreduplicates and then is transmitted to gametes. In situ hybridizations performed on gonadal tissues showed that the RlYB2 transcript was widely expressed in the ooplasm at early stages of oogenesis in both the parental species and hybrids. Interestingly, a hybridization signal, presumably related to RlYB2‐like nascent transcripts, wa...

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.

A family of interspersed repetitive elements, RlBamHI, with sequence similarity to the transposase of hAT DNA transposons, occurs in genomes of eight western Palaearctic water frog taxa and the brown frog Rana dalmatina, but was not detected in Xenopus laevis or Salamandra salamandra. RlBamHI elements are not tandemly arrayed, are dispersed across all chromosomes although not uniformly distributed, and based on dot-blot hybridizations may constitute as much as approximately 10% of the genomes of Rana lessonae, Rana ridibunda, and Rana perezi, but only approximately 1% of that of Rana saharica. Eleven nucleotide sequences of a 572 bp fragment from the nine taxa are very similar (pairwise differences 0.4–8.1% for nucleotides, 0.6–18.2% for amino acids), and all share a single open reading frame across the whole RlBamHI fragment. The reading frame is maintained despite several indels, most of which are multiples of 3 bp, but a pair of which in one species alters and, after 13 codons, restores the reading frame. It is possible that the reading frame is selectively maintained, suggesting recent or even present transposition capacity. The amino acid sequences encoded by RlBamHI elements, but not the nucleotide sequences themselves, reveal the similarity of RlBamHI to members of the widespread hAT superfamily of DNA transposons. Amino acid sequence comparisons permitted no convincing phylogenetic placement of RlBamHI among 32 representative hAT transposons across organisms, probably reflecting both the genetic divergence of RlBamHI elements and the paucity of aligned residues available for analysis. Among DNA transposons, a genomic fraction of 10% is extraordinarily high. Phylogenetic analyses of Rana RlBamHI sequences match poorly with independent genetic and molecular phylogenies; the elements compared probably are not orthologous, which renders their sequences inadequate for reconstructing organismal phylogenetic histories.

Carbone, Lucia et al.-- This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported licence., Gibbons are small arboreal apes that display an accelerated rate of evolutionary chromosomal rearrangement and occupy a key node in the primate phylogeny between Old World monkeys and great apes. Here we present the assembly and analysis of a northern white-cheeked gibbon (Nomascus leucogenys) genome. We describe the propensity for a gibbon-specific retrotransposon (LAVA) to insert into chromosome segregation genes and alter transcription by providing a premature termination site, suggesting a possible molecular mechanism for the genome plasticity of the gibbon lineage. We further show that the gibbon genera (Nomascus, Hylobates, Hoolock and Symphalangus) experienced a near-instantaneous radiation ∼5 million years ago, coincident with major geographical changes in southeast Asia that caused cycles of habitat compression and expansion. Finally, we identify signatures of positive selection in genes important for forelimb development (TBX5) and connective tissues (COL1A1) that may have been involved in the adaptation of gibbons to their arboreal habitat.

Gene duplication is a major driver of organismal adaptation and evolution and plays an important role in multiple human diseases. Whole-genome analyses have shown similar and high rates of gene duplication across a variety of eukaryotic species. Most of these studies, however, did not address the possible impact of interlocus gene conversion (IGC) on the evolution of gene duplicates. Because IGC homogenizes pairs of duplicates, widespread conversion would cause gene duplication events that happened long ago to appear more recent, resulting in artificially high estimates of duplication rates. Although the majority of genome-wide studies (including in the budding yeast Saccharomyces cerevisiae [Scer]) point to levels of IGC between paralogs ranging from 2% to 18%, Gao and Innan (Gao LZ, Innan H. 2004. Very low gene duplication rate in the yeast genome. Science 306:1367-1370.) found that gene conversion in yeast affected >80% of paralog pairs. If conversion rates really are this high, it would imply that the rate of gene duplication in eukaryotes is much lower than previously reported. In this work, we apply four different methodologies-including one approach that closely mirrors Gao and Innan's method-to estimate the level of IGC in Scer. Our analyses point to a maximum conversion level of 13% between paralogs in this species, in close agreement with most estimates of IGC in eukaryotes. We also show that the exceedingly high levels of conversion found previously derive from application of an accurate method to an inappropriate data set. In conclusion, our work provides the most striking evidence to date supporting the reduced incidence of IGC among Scer paralogs and sets up a framework for future analyses in other eukaryotes.

In many organisms, the specification of cell fate and the formation of embryonic axes depend on a proper distribution of maternal mRNAs during oogenesis. Asymmetrically localized determinants are required both for embryonic axes and germline determination in anuran amphibians. As a model system of these processes, we have used a species complex of the genus Pelophylax (Rana), characterized by a hybridogenetic reproduction that involves events of genome exclusion and endoreduplication during meiosis in both sexes. With the aim of characterizing the still largely unknown molecular events regulating Pelophylax gametogenesis, we have isolated in this animal model homologues of the deleted in AZoospermia-like (DAZl) and pumilio gene families (named RlDazl and RlPum1, respectively), which encode posttranscriptional regulators. Expression pattern analysis of these genes showed that RlDazl is exclusively expressed in gonadal tissues, whereas RlPum1 is expressed in both somatic tissues and gonads. In situ hybridization carried out on gonads revealed that the two transcripts were asymmetrically localized along the animal–vegetal (A–V) axis of oocytes. In particular, the RlDazl transcript progressively collected to the vegetal pole during oogenesis, whereas the RlPum1 mRNA was preferentially enriched at the animal hemisphere. In adult testes, RlDazl and RlPum1 were expressed in specific phases of spermatogenetic divisions as shown by immunostaining with anti-H3 phosphohistone antibody. Our results indicate that RlDazl and RlPum1 represent two early indicators of oocyte polarity in this hybridogenetic vertebrate model. Additionally, RlDazl share with vertebrate DAZ- like genes a germ cell-specific expression pattern. J. Exp. Zool. (Mol. Dev. Evol.) 316:330–338, 2011. © 2011 Wiley-Liss, Inc.

5 páginas, 5 figuras, 2 tablas.-- This paper is distributed under the terms of the Creative Commons Attribution-Non-Commercial-Share Alike Licence, and is freely available to all readers atwww.nature.com/nature.-- et al., ‘Orang-utan’ is derived from a Malay term meaning ‘man of the forest’ and aptly describes the southeast Asian great apes native to Sumatra and Borneo. The orang-utan species, Pongo abelii (Sumatran) and Pongo pygmaeus (Bornean), are the most phylogenetically distant great apes from humans, thereby providing an informative perspective on hominid evolution. Here we present a Sumatran orang-utan draft genome assembly and short read sequence data from five Sumatran and five Bornean orang-utan genomes. Our analyses reveal that, compared to other primates, the orang-utan genome has many unique features. Structural evolution of the orang-utan genome has proceeded much more slowly than other great apes, evidenced by fewer rearrangements, less segmental duplication, a lower rate of gene family turnover and surprisingly quiescent Alu repeats, which have played a major role in restructuring other primate genomes. We also describe a primate polymorphic neocentromere, found in both Pongo species, emphasizing the gradual evolution of orang-utan genome structure. Orang-utans have extremely low energy usage for a eutherian mammal1, far lower than their hominid relatives. Adding their genome to the repertoire of sequenced primates illuminates new signals of positive selection in several pathways including glycolipid metabolism. From the population perspective, both Pongo species are deeply diverse; however, Sumatran individuals possess greater diversity than their Bornean counterparts, and more species-specific variation. Our estimate of Bornean/Sumatran speciation time, 400,000 years ago, is more recent than most previous studies and underscores the complexity of the orang-utan speciation process. Despite a smaller modern census population size, the Sumatran effective population size (Ne) expanded exponentially relative to the ancestral Ne after the split, while Bornean Ne declined over the same period. Overall, the resources and analyses presented here offer new opportunities in evolutionary genomics, insights into hominid biology, and an extensive database of variation for conservation efforts., The orang-utan genome project was funded by the National Human Genome Research Institute (NHGRI), including grantsU54HG003079 (R.K.W.) and U54 HG003273 (R.A.G), with further support from National Institutes of Health R01 GM59290 (M.A.B.), PO1 AG022064 (M.A.B.), HG002385 (E.E.E.) and HG002238 (W.M.), National Science Foundation DBI-0644111 (A.S. and B.B.), David and Lucile Packard Foundation(A.S., V.T. and T.V.),CornellUniversity Provost’s Fellowship (A.L.M.), UK Medical Research Council (C.P.P., G.L., S.M. and A.H.), Marie Curie Fellowship (T.M.-B.), Ministerio de Ciencia e Innovación-Spain (MCI-Spain) and Fundación M. Botín (V.Q., X.S.P., G.R.O. and C.L.-O.), MCI-Spain BFU2006-15413-C02-01 and BFU2009-13409-C02-02 (A.N.), Spanish National Institute for Bioinformatics (INAB) and Fundação para a Ciência e a Tecnologia (Portugal), SFRH/BPD/26384/2006(R.F.) and SFRH/BD/15856/2005 (O.F.), PRIN and CEGBA (M.R., N.A. and G.D.V.), and the Commission of the European Communities IRG-224885 (T.V.), IRG-231025 (B.B.).

Nonallelic gene conversion has been proposed as a major force in homogenizing the sequences of paralogous genes. In this work, we investigate the extent and characteristics of gene conversion among gene families in nine species of the genus Drosophila. We carried out a genome-wide study of 2855 gene families (including 17,742 genes) and determined that conversion events involved 2628 genes. The proportion of converted genes ranged across species from 1 to 9% when paralogs of all ages were included. Although higher levels of gene conversion were found among young gene duplicates, at most 1–2% of the coding sequences of these duplicates were affected by conversion. Using a second approach relying on gene family size changes and gene-tree/species-tree reconciliation methods, we estimate that only 1–15% of gene trees are misled by gene conversion, depending on the lineage considered. Several features of paralogous genes correlate with gene conversion, such as intra-/interchromosomal location, level of nucleotide divergence, and GC content, although we found no definitive evidence for biased substitution patterns. After considering species-specific differences in the age and distance between paralogs, we found a highly significant difference in the amount of gene conversion among species. In particular, members of the melanogaster group showed the lowest proportion of converted genes. Our data therefore suggest underlying differences in the mechanistic basis of gene conversion among species.

Duplicate genes act as a source of genetic material from which new functions arise. They exist in large numbers in every sequenced eukaryotic genome and may be responsible for many differences in phenotypes between species. However, recent work searching for the targets of positive selection in humans has largely ignored duplicated genes due to complications in orthology assignment. Here we find that a high proportion of young gene duplicates in the human, macaque, mouse, and rat genomes have experienced adaptive natural selection. Approximately 10% of all lineage-specific duplicates show evidence for positive selection on their protein sequences, larger than any reported amount of selection among single-copy genes in these lineages using similar methods. We also find that newly duplicated genes that have been transposed to new chromosomal locations are significantly more likely to have undergone positive selection than the ancestral copy. Human-specific duplicates evolving under adaptive natural selection include a surprising number of genes involved in neuronal and cognitive functions. Our results imply that genome scans for selection that ignore duplicated loci are missing a large fraction of all adaptive substitutions. The results are also in agreement with the classical model of evolution by gene duplication, supporting a common role for neofunctionalization in the long-term maintenance of gene duplicates.

Genes of the Schlafen family, first discovered in mouse, are expressed in hematopoietic cells and are involved in immune processes. Previous results showed that they are candidate genes for two major phenomena: meiotic drive and embryonic lethality (DDK syndrome). However, these genes remain poorly understood, mostly due to the limitations imposed by their similarity, close location and the potential functional redundancy of the gene family members. Here we use genomic and phylogenetic studies to investigate the evolution and role of this family of genes. Our results show that the Schlafen family is widely distributed in mammals, where we recognize four major clades that experienced lineage-specific expansions or contractions in various orders, including primates and rodents. In addition, we identified members of the Schlafen family in Chondrichthyes and Amphibia, indicating an ancient origin of these genes. We find evidence that positive selection has acted on many Schlafen genes. Moreover, our analyses indicate that a member of the Schlafen family was horizontally transferred from murine rodents to orthopoxviruses, where it is hypothesized to play a role in allowing the virus to survive host immune defense mechanisms. The functional relevance of the viral Schlafen sequences is further underscored by our finding that they are evolving under purifying selection. This is of particular importance, since orthopoxviruses infect mammals and include variola, the causative agent of smallpox, and monkeypox, an emerging virus of great concern for human health.

Previous studies have shown that recombination between allelic sequences can cause likelihood-based methods for detecting positive selection to produce many false-positive results. In this article, we use simulations to study the impact of nonallelic gene conversion on the specificity of PAML to detect positive selection among gene duplicates. Our results show that, as expected, gene conversion leads to higher rates of false-positive results, although only moderately. These rates increase with the genetic distance between sequences, the length of converted tracts, and when no outgroup sequences are included in the analysis. We also find that branch-site models will incorrectly identify unconverted sequences as the targets of positive selection when their close paralogs are converted. Bayesian prediction of sites undergoing adaptive evolution implemented in PAML is affected by conversion, albeit in a less straightforward way. Our work suggests that particular attention should be devoted to the evolutionary analysis of recent duplicates that may have experienced gene conversion because they may provide false signals of positive selection. Fortunately, these results also imply that those cases most susceptible to false-positive results--i.e., high divergence between paralogs, long conversion tracts--are also the cases where detecting gene conversion is the easiest.

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.

Plants use light to optimize growth and development. The photoreceptor phytochrome A (phyA) mediates various far-red light–induced responses. We show that Arabidopsis FHY3 and FAR1 , which encode two proteins related to Mutator -like transposases, act together to modulate phyA signaling by directly activating the transcription of FHY1 and FHL , whose products are essential for light-induced phyA nuclear accumulation and subsequent light responses. FHY3 and FAR1 have separable DNA binding and transcriptional activation domains that are highly conserved in Mutator -like transposases. Further, expression of FHY3 and FAR1 is negatively regulated by phyA signaling. We propose that FHY3 and FAR1 represent transcription factors that have been co-opted from an ancient Mutator -like transposase(s) to modulate phyA-signaling homeostasis in higher plants.

The mammalian centromere-associated protein B (CENP-B) shares significant sequence similarity with 3 proteins in fission yeast (Abp1, Cbh1, and Cbh2) that also bind centromeres and have essential function for chromosome segregation and centromeric heterochromatin formation. Each of these proteins displays extensive sequence similarity with pogo-like transposases, which have been previously identified in the genomes of various insects and vertebrates, in the protozoan Entamoeba and in plants. Based on this distribution, it has been proposed that the mammalian and fission yeast centromeric proteins are derived from "domesticated" pogo-like transposons. Here we took advantage of the vast amount of sequence information that has become recently available for a wide range of fungal and animal species to investigate the origin of the mammalian CENP-B and yeast CENP-B-like genes. A highly conserved ortholog of CENP-B was detected in 31 species of mammals, including opossum and platypus, but was absent from all nonmammalian species represented in the databases. Similarly, no ortholog of the fission yeast centromeric proteins was identified in any of the various fungal genomes currently available. In contrast, we discovered a plethora of novel pogo-like transposons in diverse invertebrates and vertebrates and in several filamentous fungi. Phylogenetic analysis revealed that the mammalian and fission yeast CENP-B proteins fall into 2 distinct monophyletic clades, each of which includes a different set of pogo-like transposons. These results are most parsimoniously explained by independent domestication events of pogo-like transposases into centromeric proteins in the mammalian and fission yeast lineages, a case of "convergent domestication." These findings highlight the propensity of transposases to give rise to new host proteins and the potential of transposons as sources of genetic innovation.

European water frog hybrids Rana esculenta (R. ridibunda×R. lessonae) reproduce hemiclonally, by hybridogenesis: in the germ line they exclude the genome of one parental species and produce haploid gametes with an unrecombined genome of the other parental species. In the widespread L-E population system, both sexes of hybrids (E) coexist with R. lessonae (L). They exclude the lessonae genome and produce ridibunda gametes. In the R-E system, hybrid males coexist with R. ridibunda (R); they exclude either their ridibunda or their lessonae genome and produce sperm with a lessonae or with a ridibunda genome or a mixture of both kinds of sperm. We examined 13 male offspring, 12 of which were from crosses between L-E system and R-E system frogs. All were somatically hybrid. With one exception, they excluded the lessonae genome in the germ line and subsequently endoreduplicated the ridibunda genome. Spermatogonial metaphases contained a haploid or a diploid number of ridibunda chromosomes, identified through in situ hybridization to a satellite DNA marker, and by spermatocyte I metaphases containing a haploid number of ridibunda bivalents. The exception, an F1 hybrid between L-E system R. lessonae and R-E system R. ridibunda, was not hybridogenetic, showed no genome exclusion, and evidenced a disturbed gametogenesis resulting from the combination of two heterospecific genomes. None of the hybridogenetic hybrids showed any cell lines excluding the ridibunda genome, the pattern most frequent in hybrids of the R-E system, unique to that system, and essential for its persistence. A particular combination of R-E system lessonae and R-E system ridibunda genomes seems necessary to induce the R-E system type of hemiclonal gametogenesis.

Germline cell fate decisions are primarily controlled at the post-transcriptional level with DEAD-box RNA helicases playing a crucial role in germline development. In this study, we report the identification of two DEAD-box vasa/PL10 orthologues (RlVlg and RlPL10) in a species complex of the genus Rana, characterized by hybridogenetic reproduction, an enigmatic process that involves the exclusion of an individual genome, and endoreduplication events. Both genes were expressed during the early stages of gametogenesis of R. ridibunda, R. lessonae, and their natural hybrid R. esculenta. RlVlg expression was germline specific. On the other hand, RlPL10 was also expressed in somatic tissues, although only at low levels. The two genes were expressed in different phases of mitotic and meiotic spermatogenetic divisions as demonstrated by immunostaining with an anti-H3 phosphohistone antibody. The data indicate that RlVlg and RlPL10 may represent useful markers for dissecting the molecular aspects of genome exclusion and endoreduplication of the hybridogenetic gametogenesis.

In the green frog hybrids there exists an exclusive cytogenetic reproductive mechanism, referred to as hybridogenesis, whose evolutionary value is still under debate. As a contribution to the knowledge of the possible link between reproductive strategies and phylogenetic relationships, the analysis of repetitive DNA elements in the genomes of the Balkan green frogs, Rana epeirotica and R. shqiperica was carried out to compare their molecular DNA organization with those already known for R. ridibunda and R. lessonae. The Balkan species were chosen, because they give rise to spontaneous, although infrequent and not hybridogenetic hybrids with R. ridibunda. The elements of the RrSl family, a centromeric satellite DNA previously isolated in R. ridibunda, are clustered at centromeric regions also in both Balkan species, although the molecular organization of the RrSl sequences proved to be most similar in R. ridibunda and R. epeirotica. Furthermore, a repetitive DNA family (R1BamHI) is interspersed th...

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.