Transposable elements as genetic accelerators of evolution: contribution to genome size, gene regulatory network rewiring and morphological innovation

Genome-wide expansion and reorganization during grass evolution: from 30 Mb chromosomes in rice and Brachypodium to 550 Mb in Avena

BACKGROUND

The BOP (Bambusoideae, Oryzoideae, and Pooideae) clade of the Poaceae has a common ancestor, with similarities to the genomes of rice, Oryza sativa (2n = 24; genome size 389 Mb) and Brachypodium, Brachypodium distachyon (2n = 10; 271 Mb). We exploit chromosome-scale genome assemblies to show the nature of genomic expansion, structural variation, and chromosomal rearrangements from rice and Brachypodium, to diploids in the tribe Aveneae (e.g., Avena longiglumis, 2n = 2x = 14; 3,961 Mb assembled to 3,850 Mb in chromosomes).

RESULTS

Most of the Avena chromosome arms show relatively uniform expansion over the 10-fold to 15-fold genome-size increase. Apart from non-coding sequence diversification and accumulation around the centromeres, blocks of genes are not interspersed with blocks of repeats, even in subterminal regions. As in the tribe Triticeae, blocks of conserved synteny are seen between the analyzed species with chromosome fusion, fission, and nesting (insertion) events showing deep evolutionary conservation of chromosome structure during genomic expansion. Unexpectedly, the terminal gene-rich chromosomal segments (representing about 50 Mb) show translocations between chromosomes during speciation, with homogenization of genome-specific repetitive elements within the tribe Aveneae. Newly-formed intergenomic translocations of similar extent are found in the hexaploid A. sativa.

CONCLUSIONS

The study provides insight into evolutionary mechanisms and speciation in the BOP clade, which is valuable for measurement of biodiversity, development of a clade-wide pangenome, and exploitation of genomic diversity through breeding programs in Poaceae.

Retrotransposons and Telomeres

Transposable elements (TEs) comprise a significant part of eukaryotic genomes being a major source of genome instability and mutagenesis. Cellular defense systems suppress the TE expansion at all stages of their life cycle. Piwi proteins and Piwi-interacting RNAs (piRNAs) are key elements of the anti-transposon defense system, which control TE activity in metazoan gonads preventing inheritable transpositions and developmental defects. In this review, we discuss various regulatory mechanisms by which small RNAs combat TE activity. However, active transposons persist, suggesting these powerful anti-transposon defense mechanisms have a limited capacity. A growing body of evidence suggests that increased TE activity coincides with genome reprogramming and telomere lengthening in different species. In the Drosophila fruit fly, whose telomeres consist only of retrotransposons, a piRNA-mediated mechanism is required for telomere maintenance and their length control. Therefore, the efficacy of protective mechanisms must be finely balanced in order not only to suppress the activity of transposons, but also to maintain the proper length and stability of telomeres. Structural and functional relationship between the telomere homeostasis and LINE1 retrotransposon in human cells indicates a close link between selfish TEs and the vital structure of the genome, telomere. This relationship, which permits the retention of active TEs in the genome, is reportedly a legacy of the retrotransposon origin of telomeres. The maintenance of telomeres and the execution of other crucial roles that TEs acquired during the process of their domestication in the genome serve as a type of payment for such a "service."

Repetitive DNA sequence detection and its role in the human genome

Repetitive DNA sequences playing critical roles in driving evolution, inducing variation, and regulating gene expression. In this review, we summarized the definition, arrangement, and structural characteristics of repeats. Besides, we introduced diverse biological functions of repeats and reviewed existing methods for automatic repeat detection, classification, and masking. Finally, we analyzed the type, structure, and regulation of repeats in the human genome and their role in the induction of complex diseases. We believe that this review will facilitate a comprehensive understanding of repeats and provide guidance for repeat annotation and in-depth exploration of its association with human diseases.

Transposons contribute to the acquisition of cell type-specific cis-elements in the brain

Mammalian brains have evolved in stages over a long history to acquire higher functions. Recently, several transposable element (TE) families have been shown to evolve into cis-regulatory elements of brain-specific genes. However, it is not fully understood how TEs are important for gene regulatory networks. Here, we performed a single-cell level analysis using public data of scATAC-seq to discover TE-derived cis-elements that are important for specific cell types. Our results suggest that DNA elements derived from TEs, MER130 and MamRep434, can function as transcription factor-binding sites based on their internal motifs for Neurod2 and Lhx2, respectively, especially in glutamatergic neuronal progenitors. Furthermore, MER130- and MamRep434-derived cis-elements were amplified in the ancestors of Amniota and Eutheria, respectively. These results suggest that the acquisition of cis-elements with TEs occurred in different stages during evolution and may contribute to the acquisition of different functions or morphologies in the brain.

Splashed E-box and AP-1 motifs cooperatively drive regeneration response and shape regeneration abilities

Injury triggers a genetic program that induces gene expression for regeneration. Recent studies have identified regeneration-response enhancers (RREs); however, it remains unclear whether a common mechanism operates in these RREs. We identified three RREs from the zebrafish fn1b promoter by searching for conserved sequences within the surrounding genomic regions of regeneration-induced genes and performed a transgenic assay for regeneration response. Two regions contained in the transposons displayed RRE activity when combined with the -0.7 kb fn1b promoter. Another non-transposon element functioned as a stand-alone enhancer in combination with a minimum promoter. By searching for transcription factor-binding motifs and validation by transgenic assays, we revealed that the cooperation of E-box and activator protein 1 motifs is necessary and sufficient for regenerative response. Such RREs respond to variety of tissue injuries, including those in the zebrafish heart and Xenopus limb buds. Our findings suggest that the fidelity of regeneration response is ensured by the two signals evoked by tissue injuries. It is speculated that a large pool of potential enhancers in the genome has helped shape the regenerative capacities during evolution.

Chemotherapy Drugs Act Differently in the Expression and Somatic Mobilization of the mariner Transposable Element in Drosophila simulans

Transposable elements (TEs) are abundant in genomes. Their mobilization can lead to genetic variability that is useful for evolution, but can also have deleterious biological effects. Somatic mobilization (SM) has been linked to degenerative diseases, such as Alzheimer's disease and cancer. We used a Drosophila simulans strain, in which SM can be measured by counting red spots in the eyes, to investigate how chemotherapeutic agents affect expression and SM of the mariner TE. Flies were treated with Cisplatin, Dacarbazine, and Daunorubicin. After acute exposure, relative expression of mariner was quantified by RT-qPCR and oxidative stress was measured by biochemical assays. Exposure to 50 and 100 µg/mL Cisplatin increased mariner expression and ROS levels; catalase activity increased at 100 µg/mL. With chronic exposure, the number of spots also increased, indicating higher mariner SM. Dacarbazine (50 and 100 µg/mL) did not significantly alter mariner expression or mobilization or ROS levels, but decreased catalase activity (100 µg/mL). Daunorubicin (25 and 50 µM) increased mariner expression, but decreased mariner SM. ROS and catalase activity were also reduced. Our data suggest that stress factors may differentially affect the expression and SM of TEs. The increase in mariner transposase gene expression is necessary, but not sufficient for mariner SM.

Evolutionary perspective of Big tau structure: 4a exon variants of MAPT

The MAPT gene encoding the microtubule-associated protein tau can generate multiple isoforms by alternative splicing giving rise to proteins which are differentially expressed in specific areas of the nervous system and at different developmental stages. Tau plays important roles in modulating microtubule dynamics, axonal transport, synaptic plasticity, and DNA repair, and has also been associated with neurodegenerative diseases (tauopathies) including Alzheimer's disease and frontotemporal dementia. A unique high-molecular-weight isoform of tau, originally found to be expressed in the peripheral nervous system and projecting neurons, has been termed Big tau and has been shown to uniquely contain the large exon 4a that significantly increases the size and 3D structure of tau. With little progress since the original discovery of Big tau, more than 25 years ago, we have now completed a comprehensive comparative study to analyze the structure of the MAPT gene against available databases with respect to the composition of the tau exons as they evolved from early vertebrates to primates and human. We focused the analysis on the evolution of the 4a exon variants and their homology relative to humans. We discovered that the 4a exon defining Big tau appears to be present early in vertebrate evolution as a large insert that dramatically changed the size of the tau protein with low sequence conservation despite a stable size range of about 250aa, and in some species a larger 4a-L exon of 355aa. We suggest that 4a exon variants evolved independently in different species by an exonization process using new alternative splicing to address the growing complexities of the evolving nervous systems. Thus, the appearance of a significantly larger isoform of tau independently repeated itself multiple times during evolution, accentuating the need across vertebrate species for an elongated domain that likely endows Big tau with novel physiological functions as well as properties related to neurodegeneration.

Eusocial Transition in Blattodea: Transposable Elements and Shifts of Gene Expression

(1) Unravelling the molecular basis underlying major evolutionary transitions can shed light on how complex phenotypes arise. The evolution of eusociality, a major evolutionary transition, has been demonstrated to be accompanied by enhanced gene regulation. Numerous pieces of evidence suggest the major impact of transposon insertion on gene regulation and its role in adaptive evolution. Transposons have been shown to be play a role in gene duplication involved in the eusocial transition in termites. However, evidence of the molecular basis underlying the eusocial transition in Blattodea remains scarce. Could transposons have facilitated the eusocial transition in termites through shifts of gene expression? (2) Using available cockroach and termite genomes and transcriptomes, we investigated if transposons insert more frequently in genes with differential expression in queens and workers and if those genes could be linked to specific functions essential for eusocial transition. (3) The insertion rate of transposons differs among differentially expressed genes and displays opposite trends between termites and cockroaches. The functions of termite transposon-rich queen- and worker-biased genes are related to reproduction and ageing and behaviour and gene expression, respectively. (4) Our study provides further evidence on the role of transposons in the evolution of eusociality, potentially through shifts in gene expression.

The Brassicaceae genome resource (TBGR): A comprehensive genome platform for Brassicaceae plants

The Brassicaceae is an important plant family. We built a user-friendly, web-based, comparative, and functional genomic database, The Brassicaceae Genome Resource (TBGR, http://www.tbgr.org.cn), based on 82 released genomes from 27 Brassicaceae species. The TBGR database contains a large number of important functional genes, including 4,096 glucosinolate genes, 6,625 auxin genes, 13,805 flowering genes, 36,632 resistance genes, 1,939 anthocyanin genes, and 1,231 m6A genes. A total of 1,174,049 specific guide sequences for clustered regularly interspaced short palindromic repeats and 5,856,479 transposable elements were detected in Brassicaceae. TBGR also provides information on synteny, duplication, and orthologs for 27 Brassicaceae species. The TBGR database contains 1,183,851 gene annotations obtained using the TrEMBL, Swiss-Prot, Nr, GO, and Pfam databases. The BLAST, Synteny, Primer Design, Seq_fetch, and JBrowse tools are provided to help users perform comparative genomic analyses. All the genome assemblies, gene models, annotations, and bioinformatics results can be easily downloaded from the TBGR database. We plan to improve and continuously update the database with newly assembled genomes and comparative genomic studies. We expect the TBGR database to become a key resource for the study of the Brassicaceae.

ExplorATE: A new pipeline to explore active transposable elements from RNA-seq data

MOTIVATION

Transposable elements (TEs) are ubiquitous in genomes and many remain active. TEs comprise an important fraction of the transcriptomes with potential effects on the host genome, either by generating deleterious mutations or promoting evolutionary novelties. However, their functional study is limited by the difficulty in their identification and quantification, particularly in non-model organisms.

RESULTS

We developed a new pipeline (ExplorATE or Explore Active Transposable Elements) implemented in R and bash that allows the quantification of active TEs in both model and non-model organisms. ExplorATE creates TE-specific indexes and uses the Selective Alignment (SA) to filter out co-transcribed transposons within genes based on alignment scores. Moreover, our software incorporates a Wicker-like criteria to refine a set of target TEs and avoid spurious mapping. Based on simulated and real data, we show that the SA strategy adopted by ExplorATE achieved better estimates of non-co-transcribed elements than other available alignment-based or mapping-based software. ExplorATE results showed high congruence with alignment-based tools with and without a reference genome, yet ExplorATE required less execution time. Likewise, ExplorATE expands and complements most previous TE analyses by incorporating the co-transcription and multi-mapping effects during quantification, and provides a seamless integration with other downstream tools within the R environment.

AVAILABILITY

Source code is available at https://github.com/FemeniasM/ExplorATEproject and https://github.com/FemeniasM/ExplorATE_shell_script.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Somatic Mobilization: High Somatic Insertion Rate of mariner Transposable Element in Drosophila simulans

Although transposable elements (TEs) are usually silent in somatic tissues, they are sometimes mobilized in the soma and can potentially have biological consequences. The mariner element is one of the TEs involved in somatic mobilization (SM) in Drosophila and has a high rate of somatic excision. It is also known that temperature is an important factor in the increase of the mariner element SM in the fly. However, it is important to emphasize that excision is only one step of TE transposition, and the final step in this process is insertion. In the present study, we used an assay based on sequencing of the mariner flanking region and developed a pipeline to identify novel mariner insertions in Drosophila simulans at 20 and 28 °C. We found that flies carrying two mariner copies (one autonomous and one non-autonomous) had an average of 236.4 (±99.3) to 279 (±107.7) new somatic insertions at 20 °C and an average of 172.7 (±95.3) to 252.6 (±67.3) at 28 °C. In addition, we detected fragments containing mariner and others without mariner in the same regions with low-coverage long-read sequencing, indicating the process of excision and insertion. In conclusion, a low number of autonomous copies of the mariner transposon can promote a high rate of new somatic insertions during the developmental stages of Drosophila. Additionally, the developed method seems to be sensitive and adequate for the verification and estimation of somatic insertion.

The Evolution of Cytogenetic Traits in Cuscuta (Convolvulaceae), the Genus With the Most Diverse Chromosomes in Angiosperms

Karyotypes are characterized by traits such as chromosome number, which can change through whole-genome duplication and dysploidy. In the parasitic plant genus Cuscuta (Convolvulaceae), chromosome numbers vary more than 18-fold. In addition, species of this group show the highest diversity in terms of genome size among angiosperms, as well as a wide variation in the number and distribution of 5S and 35S ribosomal DNA (rDNA) sites. To understand its karyotypic evolution, ancestral character state reconstructions were performed for chromosome number, genome size, and position of 5S and 35S rDNA sites. Previous cytogenetic data were reviewed and complemented with original chromosome counts, genome size estimates, and rDNA distribution assessed via fluorescence in situ hybridization (FISH), for two, seven, and 10 species, respectively. Starting from an ancestral chromosome number of x = 15, duplications were inferred as the prevalent evolutionary process. However, in holocentric clade (subgenus Cuscuta), dysploidy was identified as the main evolutionary mechanism, typical of holocentric karyotypes. The ancestral genome size of Cuscuta was inferred as approximately 1C = 12 Gbp, with an average genome size of 1C = 2.8 Gbp. This indicates an expansion of the genome size relative to other Convolvulaceae, which may be linked to the parasitic lifestyle of Cuscuta. Finally, the position of rDNA sites varied mostly in species with multiple sites in the same karyotype. This feature may be related to the amplification of rDNA sites in association to other repeats present in the heterochromatin. The data suggest that different mechanisms acted in different subgenera, generating the exceptional diversity of karyotypes in Cuscuta.

Impact of transposable elements on the evolution of complex living systems and their epigenetic control

Transposable elements (TEs) contribute to genomic innovations, as well as genome instability, across a wide variety of species. Popular designations such as 'selfish DNA' and 'junk DNA,' common in the 1980s, may be either inaccurate or misleading, while a more enlightened view of the TE-host relationship covers a range from parasitism to mutualism. Both plant and animal hosts have evolved epigenetic mechanisms to reduce the impact of TEs, both by directly silencing them and by reducing their ability to transpose in the genome. However, TEs have also been co-opted by both plant and animal genomes to perform a variety of physiological functions, ranging from TE-derived proteins acting directly in normal biological functions to innovations in transcription factor activity and also influencing gene expression. Their presence, in fact, can affect a range of features at genome, phenotype, and population levels. The impact TEs have had on evolution is multifaceted, and many aspects still remain unexplored. In this review, the epigenetic control of TEs is contextualized according to the evolution of complex living systems.

The Evolutionary Volte-Face of Transposable Elements: From Harmful Jumping Genes to Major Drivers of Genetic Innovation

Transposable elements (TEs) are self-replicating DNA elements that constitute major fractions of eukaryote genomes. Their ability to transpose can modify the genome structure with potentially deleterious effects. To repress TE activity, host cells have developed numerous strategies, including epigenetic pathways, such as DNA methylation or histone modifications. Although TE neo-insertions are mostly deleterious or neutral, they can become advantageous for the host under specific circumstances. The phenomenon leading to the appropriation of TE-derived sequences by the host is known as TE exaptation or co-option. TE exaptation can be of different natures, through the production of coding or non-coding DNA sequences with ultimately an adaptive benefit for the host. In this review, we first give new insights into the silencing pathways controlling TE activity. We then discuss a model to explain how, under specific environmental conditions, TEs are unleashed, leading to a TE burst and neo-insertions, with potential benefits for the host. Finally, we review our current knowledge of coding and non-coding TE exaptation by providing several examples in various organisms and describing a method to identify TE co-option events.

Development of infrastructure for a systemic multidisciplinary approach to study aging in retired sled dogs

Canines represent a valuable model for mammalian aging studies as large animals with short lifespans, allowing longitudinal analyses within a reasonable time frame. Moreover, they develop a spectrum of aging-related diseases resembling that of humans, are exposed to similar environments, and have been reasonably well studied in terms of physiology and genetics. To overcome substantial variables that complicate studies of privately-owned household dogs, we have focused on a more uniform population composed of retired Alaskan sled dogs that shared similar lifestyles, including exposure to natural stresses, and are less prone to breed-specific biases than a pure breed population. To reduce variability even further, we have collected a population of 103 retired (8-11 years-old) sled dogs from multiple North American kennels in a specialized research facility named Vaika. Vaika dogs are maintained under standardized conditions with professional veterinary care and participate in a multidisciplinary program to assess the longitudinal dynamics of aging. The established Vaika infrastructure enables periodic gathering of quantitative data reflecting physical, physiological, immunological, neurological, and cognitive decline, as well as monitoring of aging-associated genetic and epigenetic alterations occurring in somatic cells. In addition, we assess the development of age-related diseases such as arthritis and cancer. In-depth data analysis, including artificial intelligence-based approaches, will build a comprehensive, integrated model of canine aging and potentially identify aging biomarkers that will allow use of this model for future testing of antiaging therapies.

B2 SINE Copies Serve as a Transposable Boundary of DNA Methylation and Histone Modifications in the Mouse

More than one million copies of short interspersed elements (SINEs), a class of retrotransposons, are present in the mammalian genomes, particularly within gene-rich genomic regions. Evidence has accumulated that ancient SINE sequences have acquired new binding sites for transcription factors (TFs) through multiple mutations following retrotransposition, and as a result have rewired the host regulatory network during the course of evolution. However, it remains unclear whether currently active SINEs contribute to the expansion of TF binding sites. To study the mobility, expression, and function of SINE copies, we first identified about 2,000 insertional polymorphisms of SINE B1 and B2 families within Mus musculus. Using a novel RNA sequencing method designated as melRNA-seq, we detected the expression of SINEs in male germ cells at both the subfamily and genomic copy levels: the vast majority of B1 RNAs originated from evolutionarily young subfamilies, whereas B2 RNAs originated from both young and old subfamilies. DNA methylation and chromatin immunoprecipitation-sequencing (ChIP-seq) analyses in liver revealed that polymorphic B2 insertions served as a boundary element inhibiting the expansion of DNA hypomethylated and histone hyperacetylated regions, and decreased the expression of neighboring genes. Moreover, genomic B2 copies were enriched at the boundary of various histone modifications, and chromatin insulator protein, CCCTC-binding factor, a well-known chromatin boundary protein, bound to >100 polymorphic and >10,000 non-polymorphic B2 insertions. These results suggest that the currently active B2 copies are mobile boundary elements that can modulate chromatin modifications and gene expression, and are likely involved in epigenomic and phenotypic diversification of the mouse species.

More than causing (epi)genomic instability: emerging physiological implications of transposable element modulation

Transposable elements (TEs) initially attracted attention because they comprise a major portion of the genomic sequences in plants and animals. TEs may jump around the genome and disrupt both coding genes as well as regulatory sequences to cause disease. Host cells have therefore evolved various epigenetic and functional RNA-mediated mechanisms to mitigate the disruption of genomic integrity by TEs. TE associated sequences therefore acquire the tendencies of attracting various epigenetic modifiers to induce epigenetic alterations that may spread to the neighboring genes. In addition to posting threats for (epi)genome integrity, emerging evidence suggested the physiological importance of endogenous TEs either as cis-acting control elements for controlling gene regulation or as TE-containing functional transcripts that modulate the transcriptome of the host cells. Recent advances in long-reads sequence analysis technologies, bioinformatics and genetic editing tools have enabled the profiling, precise annotation and functional characterization of TEs despite their challenging repetitive nature. The importance of specific TEs in preimplantation embryonic development, germ cell differentiation and meiosis, cell fate determination and in driving species specific differences in mammals will be discussed.

Retrotransposons as Drivers of Mammalian Brain Evolution

Retrotransposons, a large and diverse class of transposable elements that are still active in humans, represent a remarkable force of genomic innovation underlying mammalian evolution. Among the features distinguishing mammals from all other vertebrates, the presence of a neocortex with a peculiar neuronal organization, composition and connectivity is perhaps the one that, by affecting the cognitive abilities of mammals, contributed mostly to their evolutionary success. Among mammals, hominids and especially humans display an extraordinarily expanded cortical volume, an enrichment of the repertoire of neural cell types and more elaborate patterns of neuronal connectivity. Retrotransposon-derived sequences have recently been implicated in multiple layers of gene regulation in the brain, from transcriptional and post-transcriptional control to both local and large-scale three-dimensional chromatin organization. Accordingly, an increasing variety of neurodevelopmental and neurodegenerative conditions are being recognized to be associated with retrotransposon dysregulation. We review here a large body of recent studies lending support to the idea that retrotransposon-dependent evolutionary novelties were crucial for the emergence of mammalian, primate and human peculiarities of brain morphology and function.

KRAB-ZFP Transcriptional Regulators Acting as Oncogenes and Tumor Suppressors: An Overview

Krüppel-associated box zinc finger proteins (KRAB-ZFPs) constitute the largest family of transcriptional factors exerting co-repressor functions in mammalian cells. In general, KRAB-ZFPs have a dual structure. They may bind to specific DNA sequences via zinc finger motifs and recruit a repressive complex through the KRAB domain. Such a complex mediates histone deacetylation, trimethylation of histone 3 at lysine 9 (H3K9me3), and subsequent heterochromatization. Nevertheless, apart from their repressive role, KRAB-ZFPs may also co-activate gene transcription, likely through interaction with other factors implicated in transcriptional control. KRAB-ZFPs play essential roles in various biological processes, including development, imprinting, retroelement silencing, and carcinogenesis. Cancer cells possess multiple genomic, epigenomic, and transcriptomic aberrations. A growing number of data indicates that the expression of many KRAB-ZFPs is altered in several tumor types, in which they may act as oncogenes or tumor suppressors. Hereby, we review the available literature describing the oncogenic and suppressive roles of various KRAB-ZFPs in cancer. We focused on their association with the clinicopathological features and treatment response, as well as their influence on the cancer cell phenotype. Moreover, we summarized the identified upstream and downstream molecular mechanisms that may govern the functioning of KRAB-ZFPs in a cancer setting.

Transposable Elements and Stress in Vertebrates: An Overview

Since their identification as genomic regulatory elements, Transposable Elements (TEs) were considered, at first, molecular parasites and later as an important source of genetic diversity and regulatory innovations. In vertebrates in particular, TEs have been recognized as playing an important role in major evolutionary transitions and biodiversity. Moreover, in the last decade, a significant number of papers has been published highlighting a correlation between TE activity and exposition to environmental stresses and dietary factors. In this review we present an overview of the impact of TEs in vertebrate genomes, report the silencing mechanisms adopted by host genomes to regulate TE activity, and finally we explore the effects of environmental and dietary factor exposures on TE activity in mammals, which is the most studied group among vertebrates. The studies here reported evidence that several factors can induce changes in the epigenetic status of TEs and silencing mechanisms leading to their activation with consequent effects on the host genome. The study of TE can represent a future challenge for research for developing effective markers able to detect precocious epigenetic changes and prevent human diseases.

Comprehensive genomic analysis reveals dynamic evolution of endogenous retroviruses that code for retroviral-like protein domains

Background

Endogenous retroviruses (ERVs) are remnants of ancient retroviral infections of mammalian germline cells. A large proportion of ERVs lose their open reading frames (ORFs), while others retain them and become exapted by the host species. However, it remains unclear what proportion of ERVs possess ORFs (ERV-ORFs), become transcribed, and serve as candidates for co-opted genes.

Results

We investigated characteristics of 176,401 ERV-ORFs containing retroviral-like protein domains (gag, pro, pol, and env) in 19 mammalian genomes. The fractions of ERVs possessing ORFs were overall small (~ 0.15%) although they varied depending on domain types as well as species. The observed divergence of ERV-ORF from their consensus sequences showed bimodal distributions, suggesting that a large proportion of ERV-ORFs either recently, or anciently, inserted themselves into mammalian genomes. Alternatively, very few ERVs lacking ORFs were found to exhibit similar divergence patterns. To identify candidates for ERV-derived genes, we estimated the ratio of non-synonymous to synonymous substitution rates (dN/dS) for ERV-ORFs in human and non-human mammalian pairs, and found that approximately 42% of the ERV-ORFs showed dN/dS < 1. Further, using functional genomics data including transcriptome sequencing, we determined that approximately 9.7% of these selected ERV-ORFs exhibited transcriptional potential.

Conclusions

These results suggest that purifying selection operates on a certain portion of ERV-ORFs, some of which may correspond to uncharacterized functional genes hidden within mammalian genomes. Together, our analyses suggest that more ERV-ORFs may be co-opted in a host-species specific manner than we currently know, which are likely to have contributed to mammalian evolution and diversification.

Comparative genomic analyses illuminate the distinct evolution of megabats within Chiroptera

The revision of the sub-order Microchiroptera is one of the most intriguing outcomes in recent mammalian molecular phylogeny. The unexpected sister–taxon relationship between rhinolophoid microbats and megabats, with the exclusion of other microbats, suggests that megabats arose in a relatively short period of time from a microbat-like ancestor. In order to understand the genetic mechanism underlying adaptive evolution in megabats, we determined the whole-genome sequences of two rousette megabats, Leschenault’s rousette (Rousettus leschenaultia) and the Egyptian fruit bat (R. aegyptiacus). The sequences were compared with those of 22 other mammals, including nine bats, available in the database. We identified that megabat genomes are distinct in that they have extremely low activity of SINE retrotranspositions, expansion of two chemosensory gene families, including the trace amine receptor (TAAR) and olfactory receptor (OR), and elevation of the dN/dS ratio in genes for immunity and protein catabolism. The adaptive signatures discovered in the genomes of megabats may provide crucial insight into their distinct evolution, including key processes such as virus resistance, loss of echolocation, and frugivorous feeding.