Modern genomes with retro-look: retrotransposed elements, retroposition and the origin of new genes

A retrotransposon gag-like-3 gene RTL3 and SOX-9 co-regulate the expression of COL2A1 in chondrocytes

PURPOSE

Transposable elements are known to remodel gene structure and provide a known source of genetic variation. Retrotransposon gag-like-3 (RTL3) is a mammalian retrotransposon-derived transcript (MART) whose function in the skeletal tissue is unknown. This study aimed to elucidate the biological significance of RTL3 in chondrogenesis and type-II collagen (COL2A1) gene expression in chondrocytes.

MATERIALS AND METHODS

Expression of RTL3, SOX-9 and COL2A1 mRNAs was determined by TaqMan assays and the protein expression by immunoblotting. RTL3 and Sox-9 depletion in human chondrocytes was achieved using validated siRNAs. An RTL3 mutant (∆RTL3) lacking the zinc finger domain was created using in vitro mutagenesis. Forced expression of RTL3, ∆RTL3, and SOX-9 was achieved using CMV promoter containing expression plasmids. CRISPR-Cas9 was utilized to delete Rtl3 and create a stable ATDC5^Rlt3-/- cell line. Matrix deposition and Col2a1 quantification during chondrogenesis were determined by Alcian blue staining and Sircol™ Soluble Collagen Assay, respectively.

RESULTS

RTL3 is not ubiquitously expressed but showed strong expression in cartilage, chondrocytes and synoviocytes but not in muscle, brain, or other tissues analyzed. Loss-of-function and gain-of-function studies demonstrated a critical role of RTL3 in the regulation of SOX-9 and COL2A1 expression and matrix synthesis during chondrogenesis. Both RTL3 and SOX-9 displayed co-regulated expression in chondrocytes. Gene regulatory activity of RTL3 requires the c-terminal CCHC zinc-finger binding domain.

CONCLUSIONS

Our results identify a novel regulatory mechanism of COL2A1 expression in chondrocytes that may help to further understand the skeletal development and the pathogenesis of diseases with altered COL2A1 expression.

MicroRNA retrocopies generated via L1-mediated retrotransposition in placental mammals help to reveal how their parental genes were transcribed

In mammalian genomes, most retrocopies emerged via the L1 retrotransposition machinery. The hallmarks of an L1-mediated retrocopy, i.e., the intronlessness, the presence of a 3′ poly-A tail, and the TSDs at both ends, were frequently used to identify retrotransposition events. However, most previous studies only focused on protein-coding genes as their possible parental sources and thus only a few retrocopies derived from non-coding genes were reported. Remarkably, none of them was from microRNAs. Here in this study, we found several retrocopies generated from the mir-302–367 cluster gene (MIR302CHG), and identified a novel alternatively spliced exon encoding mir-302a. The other recognized microRNA retrotransposition events are primate-specific with mir-373 and mir-498 as their parental genes. The 3′ poly-A tracts of these two retrocopy groups were directly attached to the end of the microRNA precursor homologous regions, which suggests that their parental transcripts might alternatively terminate at the end of mir-373 and mir-498. All the three parental microRNAs are highly expressed in specific tissues with elevated retrotransposon activity, such as the embryonic stem cells and the placenta. This might be the reason that our first microRNA retrocopy findings were derived from these three microRNA genes.

Characterization of testis-specific serine/threonine kinase 1-like (TSSK1-like) gene and expression patterns in diploid and triploid Pacific abalone (Haliotis discus hannai; Gastropoda; Mollusca) males

Testis-specific serine/threonine kinase 1-like (TSSK1-like), which plays important roles in late-phase spermatogenesis and male fertility, was characterized in Pacific abalone Haliotis discus hannai, an important commercial marine gastropod. Further, its expression patterns were assessed in diploid and induced triploid males showing differential degrees of testis maturation. Abalone TSSK1-like shared conserved structural features with mammalian TSSK1s and other potential metazoan orthologs, especially regarding the catalytic STKc domain. Phylogenetically, abalone TSSK1-like displayed a genetic affiliation with its molluscan TSSK1-like orthologs and human TSSK1. Additionally, abalone TSSK1-like gene showed a tetrapartite exon-intron organization, unlike the intronless structure of most amniotic tetrapodian TSSK1s. Molecular phylogenetic analysis in the metazoan lineage suggested a possible revision in the origin of the earliest ancestral TSSK1. Further, abalone TSSK1-like showed testis-predominant expression, which was significantly influenced by both age and seasonal reproductive cycles. Comparative expression analyses between diploid and triploid abalone males suggested that robust TSSK1-like expression occurred primarily at the post-meiotic stage. Additionally, RT-PCR assay indicates that mature abalone sperms retain TSSK1-like transcripts after release. Taken together, this study provides useful insights for further studies to assess male reproduction and sterility and/or partial fertility of induced male triploidy in abalone species.

Exaptation at the molecular genetic level

The realization that body parts of animals and plants can be recruited or coopted for novel functions dates back to, or even predates the observations of Darwin. S.J. Gould and E.S. Vrba recognized a mode of evolution of characters that differs from adaptation. The umbrella term aptation was supplemented with the concept of exaptation. Unlike adaptations, which are restricted to features built by selection for their current role, exaptations are features that currently enhance fitness, even though their present role was not a result of natural selection. Exaptations can also arise from nonaptations; these are characters which had previously been evolving neutrally. All nonaptations are potential exaptations. The concept of exaptation was expanded to the molecular genetic level which aided greatly in understanding the enormous potential of neutrally evolving repetitive DNA-including transposed elements, formerly considered junk DNA-for the evolution of genes and genomes. The distinction between adaptations and exaptations is outlined in this review and examples are given. Also elaborated on is the fact that such distinctions are sometimes more difficult to determine; this is a widespread phenomenon in biology, where continua abound and clear borders between states and definitions are rare.

The retrotransposon gag domain containing protein Rgag4 is an Ikaros target in the pituitary

Previous studies have established the common and critical involvement of the zinc finger protein Ikaros in lymphoid and pituitary cell development and expansion. Key to the assembly of several transcriptional networks, we have demonstrated up-regulation of Ikaros and its interacting partner the C-terminal Binding Protein (CtBP) in response to hypoxia. This prompted us to explore common transcriptional targets using a chromatin immunoprecipitate (ChIP) screen of DNA from pituitary corticotroph cells. This strategy yielded a finite list of targets common to both transcription factors that included the metalloprotease ADAMTS10. In this report, we focus on validation of a second candidate target, the retrotransposon gag domain containing protein Rgag4. We identified the ability of Ikaros to bind the Rgag4 promoter, influence its transcriptional activity, and induce endogenous gene expression. Robust expression of Rgag4 was noted in the anterior lobe of the pituitary gland which was diminished in Ikaros knockout mice. Down-regulation of Rgag4 resulted in profound reduction of hormone gene expression with diminished ACTH secretion, recapitulating the effect of Ikaros deficiency in knockout mice. The results introduce Rgag4 to the repertoire of effectors serving to couple the chromatin remodeler Ikaros with the hormonal stress response.

Cold-induced retrotransposition of fish LINEs

Classes of retrotransposons constitute a large portion of metazoan genome. There have been cases reported that genomic abundance of retrotransposons is correlated with the severity of low environmental temperatures. However, the molecular mechanisms underlying such correlation are unknown. We show here by cell transfection assays that retrotransposition (RTP) of a long interspersed nuclear element (LINE) from an Antarctic notothenioid fish Dissostichus mawsoni (dmL1) could be activated by low temperature exposure, causing increased dmL1 copies in the host cell genome. The cold-induced dmL1 propagation was demonstrated to be mediated by the mitogen-activated protein kinases (MAPK)/p38 signaling pathway, which is activated by accumulation of reactive oxygen species (ROS) in cold-stressed conditions. Surprisingly, dmL1 transfected cells showed an increase in the number of viable cells after prolonged cold exposures than non-transfected cells. Features of cold inducibility of dmL1 were recapitulated in LINEs of zebrafish origin both in cultured cell lines and tissues, suggesting existence of a common cold-induced LINE amplification in fishes. The findings reveal an important function of LINEs in temperature adaptation and provid insights into the MAPK/p38 stress responsive pathway that shapes LINE composition in fishes facing cold stresses.

Not so bad after all: retroviruses and long terminal repeat retrotransposons as a source of new genes in vertebrates

Viruses and transposable elements, once considered as purely junk and selfish sequences, have repeatedly been used as a source of novel protein-coding genes during the evolution of most eukaryotic lineages, a phenomenon called 'molecular domestication'. This is exemplified perfectly in mammals and other vertebrates, where many genes derived from long terminal repeat (LTR) retroelements (retroviruses and LTR retrotransposons) have been identified through comparative genomics and functional analyses. In particular, genes derived from gag structural protein and envelope (env) genes, as well as from the integrase-coding and protease-coding sequences, have been identified in humans and other vertebrates. Retroelement-derived genes are involved in many important biological processes including placenta formation, cognitive functions in the brain and immunity against retroelements, as well as in cell proliferation, apoptosis and cancer. These observations support an important role of retroelement-derived genes in the evolution and diversification of the vertebrate lineage.

Selective expression of sense and antisense transcripts of the sushi-ichi-related retrotransposon--derived family during mouse placentogenesis

Background

LTR-retrotransposons became functional neogenes through evolution by acquiring promoter sequences, regulatory elements and sequence modification. Mammalian retrotransposon transcripts (Mart1-9), also called sushi-ichi-related retrotransposon-homolog (SIRH) genes, are a class of Ty3/gypsy LTR-retroelements showing moderate homology to the sushi-ichi LTR-retrotransposon in pufferfish. Rtl1/Mart1 and Peg10/Mart2 expression in mouse placenta and demonstration of their functional roles during placental development exemplifies their importance in cellular processes. In this study, we analyzed all eleven mouse Mart genes from the blastocyst stage and throughout placentogenesis in order to gain information about their expression and regulation.

Results

Quantitative PCR, in situ hybridization (ISH) and immunoblotting showed various expression patterns of the 11 mouse Mart genes through different placental stages. Zcchc5/Mart3, Zcchc16/ Mart4 and Rgag1/Mart9 expression was undetectable. Rtl1/Mart1, Peg10/Mart2, Rgag4/Mart5 – Cxx1a,b,c/Mart8b,c,a gene expression was very low at the blastocyst stage. Later placental stages showed an increase of expression for Rtl1/Mart1, Rgag4/Mart5 – Cxx1a,b,c/Mart8b,c,a, the latter up to 1,489 molecules/ng cDNA at E9.5. From our recently published findings Peg10/Mart2 was the most highly expressed Mart gene. ISH demonstrated sense and antisense transcript co-localization of Rgag4/Mart5 to Cxx1a,b,c/Mart8b,c,a in trophoblast subtypes at the junctional zone, with an accumulation of antisense transcripts in the nuclei. To validate these results, we developed a TAG-aided sense/antisense transcript detection (TASA-TD) method, which verified sense and antisense transcripts for Rtl1/Mart1, Rgag4/Mart5 – Cxx1a,b,c/Mart8b,c,a. Except for Rtl1/Mart1 and Cxx1a,b/Mart8b,c all other Mart genes showed a reduced amount of antisense transcripts. Northern blot and 5′ and 3′ RACE confirmed both sense and antisense transcripts for Ldoc1/Mart7 and Cxx1a,b,c/Mart8b,c,a. Immunoblotting demonstrated a single protein throughout all placental stages for Ldoc1/Mart7, but for Cxx1a,b,c/Mart8b,c,a a switch occurred from a 57 kDa protein at E10.5 and E14.5 to a 25 kDa protein at E16.5 and E18.5.

Conclusions

RNA and protein detection of mouse Mart genes support neo-functionalization of retrotransposons in mammalian genomes. Undetectable expression of Zcchc5/Mart3, Zcchc16/Mart4 and Rgag1/Mart9 indicate no role during mouse placentogenesis. Rgag4/Mart5 to Cxx1a,b,c/Mart8b,c,a gene expression support a role for differentiation from the ectoplacental cone. Mart antisense transcripts and protein alterations predict unique and complex molecular regulation in a time directed manner throughout mouse placentogenesis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12977-015-0138-8) contains supplementary material, which is available to authorized users.

Evolution of testis-specific kinases TSSK1B and TSSK2 in primates

The testis-specific serine/threonine protein kinases TSSK1 and TSSK2 are known to be essential for male fertility, in mice. The enzymes are present in elongating spermatids, and targeted deletion of the two genes Tssk1 and Tssk2 results in dysregulation of spermiogenesis. The mouse genes are genetically closely linked, forming a Tssk1-Tssk2 tandem. In human, TSSK1 is present in the form of a pseudogene, TSSK1A, which is linked to an intact TSSK2 gene, and in the form of an intact gene, TSSK1B, which is not genetically linked to TSSK2. Studies on conservation of genes and gene function between mouse and human are relevant, to be able to use mouse models for studies on human infertility, and to evaluate possible targets for non-hormonal contraception targeting the male. Therefore, we have performed a detailed analysis of the evolution of genes encoding TSSK1 and TSSK2 among mammals, in particular among primates. This study includes functional analysis of replacement mutation K27R in TSSK2, which is frequently observed among humans. In primates, the kinase domains of TSSK1B and TSSK2 have evolved under negative selection, reflecting the importance to maintain their kinase activity. Positive selection was observed for the C-terminal domain of TSSK1B, which indicates that TSSK1B and TSSK2 may perform at least partly differential functions.

Expansion of CORE-SINEs in the genome of the Tasmanian devil

Background

The genome of the carnivorous marsupial, the Tasmanian devil (Sarcophilus harrisii, Order: Dasyuromorphia), was sequenced in the hopes of finding a cure for or gaining a better understanding of the contagious devil facial tumor disease that is threatening the species’ survival. To better understand the Tasmanian devil genome, we screened it for transposable elements and investigated the dynamics of short interspersed element (SINE) retroposons.

Results

The temporal history of Tasmanian devil SINEs, elucidated using a transposition in transposition analysis, indicates that WSINE1, a CORE-SINE present in around 200,000 copies, is the most recently active element. Moreover, we discovered a new subtype of WSINE1 (WSINE1b) that comprises at least 90% of all Tasmanian devil WSINE1s. The frequencies of WSINE1 subtypes differ in the genomes of two of the other Australian marsupial orders. A co-segregation analysis indicated that at least 66 subfamilies of WSINE1 evolved during the evolution of Dasyuromorphia. Using a substitution rate derived from WSINE1 insertions, the ages of the subfamilies were estimated and correlated with a newly established phylogeny of Dasyuromorphia. Phylogenetic analyses and divergence time estimates of mitochondrial genome data indicate a rapid radiation of the Tasmanian devil and the closest relative the quolls (Dasyurus) around 14 million years ago.

Conclusions

The radiation and abundance of CORE-SINEs in marsupial genomes indicates that they may be a major player in the evolution of marsupials. It is evident that the early phases of evolution of the carnivorous marsupial order Dasyuromorphia was characterized by a burst of SINE activity. A correlation between a speciation event and a major burst of retroposon activity is for the first time shown in a marsupial genome.

Primate genome gain and loss: a bone dysplasia, muscular dystrophy, and bone cancer syndrome resulting from mutated retroviral-derived MTAP transcripts

Diaphyseal medullary stenosis with malignant fibrous histiocytoma (DMS-MFH) is an autosomal-dominant syndrome characterized by bone dysplasia, myopathy, and bone cancer. We previously mapped the DMS-MFH tumor-suppressing-gene locus to chromosomal region 9p21-22 but failed to identify mutations in known genes in this region. We now demonstrate that DMS-MFH results from mutations in the most proximal of three previously uncharacterized terminal exons of the gene encoding methylthioadenosine phosphorylase, MTAP. Intriguingly, two of these MTAP exons arose from early and independent retroviral-integration events in primate genomes at least 40 million years ago, and since then, their genomic integration has gained a functional role. MTAP is a ubiquitously expressed homotrimeric-subunit enzyme critical to polyamine metabolism and adenine and methionine salvage pathways and was believed to be encoded as a single transcript from the eight previously described exons. Six distinct retroviral-sequence-containing MTAP isoforms, each of which can physically interact with archetype MTAP, have been identified. The disease-causing mutations occur within one of these retroviral-derived exons and result in exon skipping and dysregulated alternative splicing of all MTAP isoforms. Our results identify a gene involved in the development of bone sarcoma, provide evidence of the primate-specific evolution of certain parts of an existing gene, and demonstrate that mutations in parts of this gene can result in human disease despite its relatively recent origin.

Widespread occurrence of power-law distributions in inter-repeat distances shaped by genome dynamics

Repetitive DNA sequences derived from transposable elements (TE) are distributed in a non-random way, co-clustering with other classes of repeat elements, genes and other genomic components. In a previous work we reported power-law-like size distributions (linearity in log-log scale) in the spatial arrangement of Alu and LINE1 elements in the human genome. Here we investigate the large-scale features of the spatial arrangement of all principal classes of TEs in 14 genomes from phylogenetically distant organisms by studying the size distribution of inter-repeat distances. Power-law-like size distributions are found to be widespread, extending up to several orders of magnitude. In order to understand the emergence of this distributional pattern, we introduce an evolutionary scenario, which includes (i) Insertions of DNA segments (e.g., more recent repeats) into the considered sequence and (ii) Eliminations of members of the studied TE family. In the proposed model we also incorporate the potential for transposition events (characteristic of the DNA transposons' life-cycle) and segmental duplications. Simulations reproduce the main features of the observed size distributions. Furthermore, we investigate the effects of various genomic features on the presence and extent of power-law size distributions including TE class and age, mode of parental TE transmission, GC content, deletion and recombination rates in the studied genomic region, etc. Our observations corroborate the hypothesis that insertions of genomic material and eliminations of repeats are at the basis of power-laws in inter-repeat distances. The existence of these power-laws could facilitate the formation of the recently proposed "fractal globule" for the confined chromatin organization.

Do repeated arrays of regulatory small-RNA genes elicit genomic imprinting?: Concurrent emergence of large clusters of small non-coding RNAs and genomic imprinting at four evolutionarily distinct eutherian chromosomal loci

The basic premise of the host-defense theory is that genomic imprinting, the parent-of-origin expression of a subset of mammalian genes, derives from mechanisms originally dedicated to silencing repeated and retroviral-like sequences that deeply colonized mammalian genomes. We propose that large clusters of tandemly-repeated C/D-box small nucleolar RNAs (snoRNAs) or microRNAs represent a novel category of sequences recognized as "genomic parasites", contributing to the emergence of genomic imprinting in a subset of chromosomal regions that contain them. Such a view is supported by evidence derived from studies of the imprinted snoRNA- and/or miRNA-encoding Dlk1-Dio3, Snurf-Snrpn, Sfbmt2, and C19MC domains. While adding a new piece to the challenging puzzle of mammalian genome history, this hypothesis also reinforces the notion that dissecting the features and molecular mechanisms that discriminate between "foreign" and "endogenous" sequences is of crucial importance in the field of mammalian epigenetics.

The Gypsy Database (GyDB) of mobile genetic elements: release 2.0

This article introduces the second release of the Gypsy Database of Mobile Genetic Elements (GyDB 2.0): a research project devoted to the evolutionary dynamics of viruses and transposable elements based on their phylogenetic classification (per lineage and protein domain). The Gypsy Database (GyDB) is a long-term project that is continuously progressing, and that owing to the high molecular diversity of mobile elements requires to be completed in several stages. GyDB 2.0 has been powered with a wiki to allow other researchers participate in the project. The current database stage and scope are long terminal repeats (LTR) retroelements and relatives. GyDB 2.0 is an update based on the analysis of Ty3/Gypsy, Retroviridae, Ty1/Copia and Bel/Pao LTR retroelements and the Caulimoviridae pararetroviruses of plants. Among other features, in terms of the aforementioned topics, this update adds: (i) a variety of descriptions and reviews distributed in multiple web pages; (ii) protein-based phylogenies, where phylogenetic levels are assigned to distinct classified elements; (iii) a collection of multiple alignments, lineage-specific hidden Markov models and consensus sequences, called GyDB collection; (iv) updated RefSeq databases and BLAST and HMM servers to facilitate sequence characterization of new LTR retroelement and caulimovirus queries; and (v) a bibliographic server. GyDB 2.0 is available at http://gydb.org.

Polypteridae (Actinopterygii: Cladistia) and DANA-SINEs insertions

SINE sequences are interspersed throughout virtually all eukaryotic genomes and greatly outnumber the other repetitive elements. These sequences are of increasing interest for phylogenetic studies because of their diagnostic power for establishing common ancestry among taxa, once properly characterized. We identified and characterized a peculiar family of composite tRNA-derived short interspersed SINEs, DANA-SINEs, associated with mutational activities in Danio rerio, in a group of species belonging to one of the most basal bony fish families, the Polypteridae, in order to investigate their own inner specific phylogenetic relationships. DANA sequences were identified, sequenced and then localized, by means of fluorescent in situ hybridization (FISH), in six Polypteridae species (Polypterus delhezi, P. ornatipinnis, P. palmas, P. buettikoferi P. senegalus and Erpetoichthys calabaricus) After cloning, the sequences obtained were aligned for phylogenetic analysis, comparing them with three Dipnoan lungfish species (Protopterus annectens, P. aethiopicus, Lepidosiren paradoxa), and Lethenteron reissneri (Petromyzontidae)was used as outgroup. The obtained overlapping MP, ML and NJ tree clustered together the species belonging to the two taxonomically different Osteichthyans groups: the Polypteridae, by one side, and the Protopteridae by the other, with the monotypic genus Erpetoichthys more distantly related to the Polypterus genus comprising three distinct groups: P. palmas and P. buettikoferi, P. delhezi and P. ornatipinnis and P. senegalus. In situ hybridization with DANA probes marked along the whole chromosome arms in the metaphases of all the Polypteridae species examined.

The primate-specific microRNA gene cluster (C19MC) is imprinted in the placenta

Imprinted genes play crucial roles in mammalian development and disruption of their expression is associated with many human disorders including tumourigenesis; yet, the actual number of imprinted genes in the human genome remains a matter of debate. Here, we report on the unexpected finding that the chromosome 19 microRNA cluster (C19MC), the largest human microRNA gene cluster discovered so far, is regulated by genomic imprinting with only the paternally inherited allele being expressed in the placenta. DNA methylation profiling identified a differentially methylated region (C19MC-DMR1) that overlaps an upstream CpG-rich promoter region associated with short tandem repeats. It displays a maternal-specific methylation imprint acquired in oocytes and generates a complex population of large, compartimentalized non-coding RNA (ncRNA) species retained in close proximity to the C19MC transcription site. This occurs adjacent to, but not within, a poorly characterized nuclear Alu-rich domain. Interestingly, C19MC maps near another imprinted gene, the maternally expressed ZNF331 gene, and therefore may define a novel, previously unrecognized large imprinted primate-specific chromosomal domain. Altogether, our study adds C19MC to the growing list of imprinted repeated small RNA gene clusters and further strengthens the potential involvement of small ncRNAs in the function and/or the evolution of imprinted gene networks.

Genomic organization of eukaryotic tRNAs

BACKGROUND

Surprisingly little is known about the organization and distribution of tRNA genes and tRNA-related sequences on a genome-wide scale. While tRNA gene complements are usually reported in passing as part of genome annotation efforts, and peculiar features such as the tandem arrangements of tRNA gene in Entamoeba histolytica have been described in some detail, systematic comparative studies are rare and mostly restricted to bacteria. We therefore set out to survey the genomic arrangement of tRNA genes and pseudogenes in a wide range of eukaryotes to identify common patterns and taxon-specific peculiarities.

RESULTS

In line with previous reports, we find that tRNA complements evolve rapidly and tRNA gene and pseudogene locations are subject to rapid turnover. At phylum level, the distributions of the number of tRNA genes and pseudogenes numbers are very broad, with standard deviations on the order of the mean. Even among closely related species we observe dramatic changes in local organization. For instance, 65% and 87% of the tRNA genes and pseudogenes are located in genomic clusters in zebrafish and stickleback, resp., while such arrangements are relatively rare in the other three sequenced teleost fish genomes. Among basal metazoa, Trichoplax adherens has hardly any duplicated tRNA gene, while the sea anemone Nematostella vectensis boasts more than 17000 tRNA genes and pseudogenes. Dramatic variations are observed even within the eutherian mammals. Higher primates, for instance, have 616 +/- 120 tRNA genes and pseudogenes of which 17% to 36% are arranged in clusters, while the genome of the bushbaby Otolemur garnetti has 45225 tRNA genes and pseudogenes of which only 5.6% appear in clusters. In contrast, the distribution is surprisingly uniform across plant genomes. Consistent with this variability, syntenic conservation of tRNA genes and pseudogenes is also poor in general, with turn-over rates comparable to those of unconstrained sequence elements. Despite this large variation in abundance in Eukarya we observe a significant correlation between the number of tRNA genes, tRNA pseudogenes, and genome size.

CONCLUSIONS

The genomic organization of tRNA genes and pseudogenes shows complex lineage-specific patterns characterized by an extensive variability that is in striking contrast to the extreme levels of sequence-conservation of the tRNAs themselves. The comprehensive analysis of the genomic organization of tRNA genes and pseudogenes in Eukarya provides a basis for further studies into the interplay of tRNA gene arrangements and genome organization in general.

An overview of the processes shaping protein evolution

From a comparatively small number of protein structural domains a staggering array of structural variants has evolved which has, in turn, facilitated an expanse of functional derivatives. Herein I review the primary mechanisms which have contributed to the vastness of our existing, and expanding, protein repertoires.

Deconstructing the dogma: a new view of the evolution and genetic programming of complex organisms

Since the birth of molecular biology it has been generally assumed that most genetic information is transacted by proteins, and that RNA plays an intermediary role. This led to the subsidiary assumption that the vast tracts of noncoding sequences in the genomes of higher organisms are largely nonfunctional, despite the fact that they are transcribed. These assumptions have since become articles of faith, but they are not necessarily correct. I propose an alternative evolutionary history whereby developmental and cognitive complexity has arisen by constructing sophisticated RNA-based regulatory networks that interact with generic effector complexes to control gene expression patterns and the epigenetic trajectories of differentiation and development. Environmental information can also be conveyed into this regulatory system via RNA editing, especially in the brain. Moreover, the observations that RNA-directed epigenetic changes can be inherited raises the intriguing question: has evolution learnt how to learn?

Large-scale discovery of insertion hotspots and preferential integration sites of human transposed elements

Throughout evolution, eukaryotic genomes have been invaded by transposable elements (TEs). Little is known about the factors leading to genomic proliferation of TEs, their preferred integration sites and the molecular mechanisms underlying their insertion. We analyzed hundreds of thousands nested TEs in the human genome, i.e. insertions of TEs into existing ones. We first discovered that most TEs insert within specific ‘hotspots’ along the targeted TE. In particular, retrotransposed Alu elements contain a non-canonical single nucleotide hotspot for insertion of other Alu sequences. We next devised a method for identification of integration sequence motifs of inserted TEs that are conserved within the targeted TEs. This method revealed novel sequences motifs characterizing insertions of various important TE families: Alu, hAT, ERV1 and MaLR. Finally, we performed a global assessment to determine the extent to which young TEs tend to nest within older transposed elements and identified a 4-fold higher tendency of TEs to insert into existing TEs than to insert within non-TE intergenic regions. Our analysis demonstrates that TEs are highly biased to insert within certain TEs, in specific orientations and within specific targeted TE positions. TE nesting events also reveal new characteristics of the molecular mechanisms underlying transposition.

A global view of genomic information--moving beyond the gene and the master regulator

The current view of gene regulation in complex organisms holds that gene expression is largely controlled by the combinatoric actions of transcription factors and other regulatory proteins, some of which powerfully influence cell type. Recent large-scale studies have confirmed that cellular differentiation involves many different regulatory factors. However, other studies indicate that the genome is pervasively transcribed to produce a variety of short and long non-protein-coding RNAs, including those derived from retrotransposed sequences, which also play important roles in the epigenetic regulation of gene expression. The evidence suggests that ontogenesis requires interplay between state-specific regulatory proteins, multitasked effector complexes and target-specific RNAs that recruit these complexes to their sites of action. Moreover, the semi-continuous nature of the transcriptome prompts the reassessment of 'genes' as discrete entities and indicates that the mammalian genome might be more accurately viewed as islands of protein-coding information in a sea of cis- and trans-acting regulatory sequences.

Evolutionarily stable association of intronic snoRNAs and microRNAs with their host genes

Small nucleolar RNAs (snoRNAs) and microRNAs (miRNAs) are integral to a range of processes, including ribosome biogenesis and gene regulation. Some are intron encoded, and this organization may facilitate coordinated coexpression of host gene and RNA. However, snoRNAs and miRNAs are known to be mobile, so intron-RNA associations may not be evolutionarily stable. We have used genome alignments across 11 mammals plus chicken to examine positional orthology of snoRNAs and miRNAs and report that 21% of annotated snoRNAs and 11% of miRNAs are positionally conserved across mammals. Among RNAs traceable to the bird–mammal common ancestor, 98% of snoRNAs and 76% of miRNAs are intronic. Comparison of the most evolutionarily stable mammalian intronic snoRNAs with those positionally conserved among primates reveals that the former are more overrepresented among host genes involved in translation or ribosome biogenesis and are more broadly and highly expressed. This stability is likely attributable to a requirement for overlap between host gene and intronic snoRNA expression profiles, consistent with an ancestral role in ribosome biogenesis. In contrast, whereas miRNA positional conservation is comparable to that observed for snoRNAs, intronic miRNAs show no obvious association with host genes of a particular functional category, and no statistically significant differences in host gene expression are found between those traceable to mammalian or primate ancestors. Our results indicate evolutionarily stable associations of numerous intronic snoRNAs and miRNAs and their host genes, with probable continued diversification of snoRNA function from an ancestral role in ribosome biogenesis.

When is it time for reverse transcription to start and go?

Upon cell infection by a retrovirus, the viral DNA polymerase, called reverse transcriptase (RT), copies the genomic RNA to generate the proviral DNA flanked by two long terminal repeats (LTR). A discovery twenty years ago demonstrated that the structural viral nucleocapsid protein (NC) encoded by Gag is an essential cofactor of reverse transcription, chaperoning RT during viral DNA synthesis. However, it is only recently that NC was found to exert a control on the timing of reverse transcription, in a spatio-temporal manner. This brief review summarizes findings on the timing of reverse transcription in wild type HIV-1 and in nucleopcapsid (NC) mutants where virions contain a large amount of newly made viral DNA. This brief review also proposes some explanations of how NC may control late reverse transcription during Gag assembly in virus producer cells.

Retrotransposons revisited: the restraint and rehabilitation of parasites

Retrotransposons, mainly LINEs, SINEs, and endogenous retroviruses, make up roughly 40% of the mammalian genome and have played an important role in genome evolution. Their prevalence in genomes reflects a delicate balance between their further expansion and the restraint imposed by the host. In any human genome only a small number of LINE1s (L1s) are active, moving their own and SINE sequences into new genomic locations and occasionally causing disease. Recent insights and new technologies promise answers to fundamental questions about the biology of transposable elements.

Evaluating the protein coding potential of exonized transposable element sequences

Background

Transposable element (TE) sequences, once thought to be merely selfish or parasitic members of the genomic community, have been shown to contribute a wide variety of functional sequences to their host genomes. Analysis of complete genome sequences have turned up numerous cases where TE sequences have been incorporated as exons into mRNAs, and it is widely assumed that such 'exonized' TEs encode protein sequences. However, the extent to which TE-derived sequences actually encode proteins is unknown and a matter of some controversy. We have tried to address this outstanding issue from two perspectives: i-by evaluating ascertainment biases related to the search methods used to uncover TE-derived protein coding sequences (CDS) and ii-through a probabilistic codon-frequency based analysis of the protein coding potential of TE-derived exons.

Results

We compared the ability of three classes of sequence similarity search methods to detect TE-derived sequences among data sets of experimentally characterized proteins: 1-a profile-based hidden Markov model (HMM) approach, 2-BLAST methods and 3-RepeatMasker. Profile based methods are more sensitive and more selective than the other methods evaluated. However, the application of profile-based search methods to the detection of TE-derived sequences among well-curated experimentally characterized protein data sets did not turn up many more cases than had been previously detected and nowhere near as many cases as recent genome-wide searches have. We observed that the different search methods used were complementary in the sense that they yielded largely non-overlapping sets of hits and differed in their ability to recover known cases of TE-derived CDS. The probabilistic analysis of TE-derived exon sequences indicates that these sequences have low protein coding potential on average. In particular, non-autonomous TEs that do not encode protein sequences, such as Alu elements, are frequently exonized but unlikely to encode protein sequences.

Conclusion

The exaptation of the numerous TE sequences found in exons as bona fide protein coding sequences may prove to be far less common than has been suggested by the analysis of complete genomes. We hypothesize that many exonized TE sequences actually function as post-transcriptional regulators of gene expression, rather than coding sequences, which may act through a variety of double stranded RNA related regulatory pathways. Indeed, their relatively high copy numbers and similarity to sequences dispersed throughout the genome suggests that exonized TE sequences could serve as master regulators with a wide scope of regulatory influence.

Reviewers:

This article was reviewed by Itai Yanai, Kateryna D. Makova, Melissa Wilson (nominated by Kateryna D. Makova) and Cedric Feschotte (nominated by John M. Logsdon Jr.).