Transposable elements in gene regulation and in the evolution of vertebrate genomes

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state regulatory potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbor distinctive transcription factor binding motifs that are similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we show that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

Untangling the gordian knot: The intertwining interactions between developmental hormone signaling and epigenetic mechanisms in insects

Hormones control developmental and physiological processes, often by regulating the expression of multiple genes simultaneously or sequentially. Crosstalk between hormones and epigenetics is pivotal to dynamically coordinate this process. Hormonal signals can guide the addition and removal of epigenetic marks, steering gene expression. Conversely, DNA methylation, histone modifications and non-coding RNAs can modulate regional chromatin structure and accessibility and regulate the expression of numerous (hormone-related) genes. Here, we provide a review of the interplay between the classical insect hormones, ecdysteroids and juvenile hormones, and epigenetics. We summarize the mode-of-action and roles of these hormones in post-embryonic development, and provide a general overview of epigenetic mechanisms. We then highlight recent advances on the interactions between these hormonal pathways and epigenetics, and their involvement in development. Furthermore, we give an overview of several 'omics techniques employed in the field. Finally, we discuss which questions remain unanswered and possible avenues for future research.

The analysis of complete genome sequence and comparative genomics of Vibrio parahaemolyticus LF1113 in Hainan

Vibrio parahaemolyticus is a Gram-negative, halophilic and polymorphic coccobacillus. It is world-widely distributed and has resulted in great economic losses since its first appearance. In this study, a pathogenic strain was isolated from diseased pearl gentian grouper and identified as V. parahaemolyticus based on the sequencing results of 16S rDNA gene. In order to gain a comprehensive understanding of this isolation, the whole genome sequencing was conducted. Phylogenetic analysis of the complete genomes of 16 Vibrio species showed that LF1113, ATCC17802, ATCC33787, 2210633, FORC 004, and 160807 were the most closely related. Animal experiments demonstrated that the isolated LF1113 strain was pathogenic in a fish model. This study is the first study to describe the complete genome sequence of a V. parahaemolyticus isolate, which infected pearl gentian grouper from an outbreak in a fish factory farm in Hainan. The results will expand our understanding of genetic characteristics, pathogenesis, diagnostics and disease prevention of V. parahaemolyticus, and lay the foundation for further study.

Genome sequencing and analysis of black flounder (Paralichthys orbignyanus) reveals new insights into Pleuronectiformes genomic size and structure

Black flounder (Paralichthys orbignyanus, Pleuronectiformes) is a commercially significant marine fish with promising aquaculture potential in Argentina. Despite extensive studies on Black flounder aquaculture, its limited genetic information available hampers the crucial role genetics plays in the development of this activity. In this study, we first employed Illumina sequencing technology to sequence the entire genome of Black flounder. Utilizing two independent libraries-one from a female and another from a male-with 150 bp paired-end reads, a mean insert length of 350 bp, and over 35 X-fold coverage, we achieved assemblies resulting in a genome size of ~ 538 Mbp. Analysis of the assemblies revealed that more than 98% of the core genes were present, with more than 78% of them having more than 50% coverage. This indicates a somehow complete and accurate genome at the coding sequence level. This genome contains 25,231 protein-coding genes, 445 tRNAs, 3 rRNAs, and more than 1,500 non-coding RNAs of other types. Black flounder, along with pufferfishes, seahorses, pipefishes, and anabantid fish, displays a smaller genome compared to most other teleost groups. In vertebrates, the number of transposable elements (TEs) is often correlated with genome size. However, it remains unclear whether the sizes of introns and exons also play a role in determining genome size. Hence, to elucidate the potential factors contributing to this reduced genome size, we conducted a comparative genomic analysis between Black flounder and other teleost orders to determine if the small genomic size could be explained by repetitive elements or gene features, including the whole genome genes and introns sizes. We show that the smaller genome size of flounders can be attributed to several factors, including changes in the number of repetitive elements, and decreased gene size, particularly due to lower amount of very large and small introns. Thus, these components appear to be involved in the genome reduction in Black flounder. Despite these insights, the full implications and potential benefits of genome reduction in Black flounder for reproduction and aquaculture remain incompletely understood, necessitating further research.

Transposons in the Williams-Beuren Syndrome Critical Region are Associated with Social Behavior in Assistance Dogs

A strong signature of selection in the domestic dog genome is found in a five-megabase region of chromosome six in which four structural variants derived from transposons have previously been associated with human-oriented social behavior, such as attentional bias to social stimuli and social interest in strangers. To explore these genetic associations in more phenotypic detail-as well as their role in training success in a specialized assistance dog program-we genotyped 1001 assistance dogs from Canine Companions for Independence®, including both successful graduates and dogs released from the training program for behaviors incompatible with their working role. We collected phenotypes on each dog using puppy-raiser questionnaires, trainer questionnaires, and both cognitive and behavioral tests. Using Bayesian mixed models, we found strong associations (95% credibility intervals excluding zero) between genotypes and certain behavioral measures, including separation-related problems, aggression when challenged or corrected, and reactivity to other dogs. Furthermore, we found moderate differences in the genotypes of dogs who graduated versus those who did not; insertions in GTF2I showed the strongest association with training success (β = 0.23, CI95% = - 0.04, 0.49), translating to an odds-ratio of 1.25 for one insertion. Our results provide insight into the role of each of these four transposons in canine sociability and may inform breeding and training practices for working dog organizations. Furthermore, the observed importance of the gene GTF2I supports the emerging consensus that variation in GTF2I genotypes and expression have important consequences for social behavior broadly.

Chronically high stress hormone levels dysregulate sperm long noncoding RNAs and their embryonic microinjection alters development and affective behaviours

Previous studies on paternal epigenetic inheritance have shown that sperm RNAs play a role in this type of inheritance. The microinjection of sperm small noncoding RNAs into fertilised mouse oocytes induces reprogramming of the early embryo, which is thought to be responsible for the differences observed in adult phenotype. While sperm long noncoding RNAs (lncRNAs) have also been investigated in a previous study, their microinjection into fertilised oocytes did not yield conclusive results regarding their role in modulating brain development and adult behavioural phenotypes. Therefore, in the current study we sought to investigate this further. We used our previously established paternal corticosterone (stress hormone) model to assess sperm lncRNA expression using CaptureSeq, a sequencing technique that is more sensitive than the ones used in other studies in the field. Paternal corticosterone exposure led to dysregulation of sperm long noncoding RNA expression, which encompassed lncRNAs, circular RNAs and transposable element transcripts. Although they have limited functional annotation, bioinformatic approaches indicated the potential of these lncRNAs in regulating brain development and function. We then separated and isolated the sperm lncRNAs and performed microinjections into fertilised oocytes, to generate embryos with modulated lncRNA populations. We observed that the resulting adult offspring had lower body weight and altered anxiety and affective behavioural responses, demonstrating roles for lncRNAs in modulating development and brain function. This study provides novel insights into the roles of lncRNAs in epigenetic inheritance, including impacts on brain development and behaviours of relevance to affective disorders.

Mechanisms of Rapid Karyotype Evolution in Mammals

Chromosome reshuffling events are often a foundational mechanism by which speciation can occur, giving rise to highly derivative karyotypes even amongst closely related species. Yet, the features that distinguish lineages prone to such rapid chromosome evolution from those that maintain stable karyotypes across evolutionary time are still to be defined. In this review, we summarize lineages prone to rapid karyotypic evolution in the context of Simpson's rates of evolution-tachytelic, horotelic, and bradytelic-and outline the mechanisms proposed to contribute to chromosome rearrangements, their fixation, and their potential impact on speciation events. Furthermore, we discuss relevant genomic features that underpin chromosome variation, including patterns of fusions/fissions, centromere positioning, and epigenetic marks such as DNA methylation. Finally, in the era of telomere-to-telomere genomics, we discuss the value of gapless genome resources to the future of research focused on the plasticity of highly rearranged karyotypes.

When Dad's Stress Gets under Kid's Skin-Impacts of Stress on Germline Cargo and Embryonic Development

Multiple lines of evidence suggest that paternal psychological stress contributes to an increased prevalence of neuropsychiatric and metabolic diseases in the progeny. While altered paternal care certainly plays a role in such transmitted disease risk, molecular factors in the germline might additionally be at play in humans. This is supported by findings on changes to the molecular make up of germ cells and suggests an epigenetic component in transmission. Several rodent studies demonstrate the correlation between paternal stress induced changes in epigenetic modifications and offspring phenotypic alterations, yet some intriguing cases also start to show mechanistic links in between sperm and the early embryo. In this review, we summarise efforts to understand the mechanism of intergenerational transmission from sperm to the early embryo. In particular, we highlight how stress alters epigenetic modifications in sperm and discuss the potential for these modifications to propagate modified molecular trajectories in the early embryo to give rise to aberrant phenotypes in adult offspring.

Transposable elements are associated with the variable response to influenza infection

Influenza A virus (IAV) infections are frequent every year and result in a range of disease severity. Here, we wanted to explore the potential contribution of transposable elements (TEs) to the variable human immune response. Transcriptome profiling in monocyte-derived macrophages from 39 individuals following IAV infection revealed significant inter-individual variation in viral load post-infection. Using transposase-accessible chromatin using sequencing (ATAC-seq), we identified a set of TE families with either enhanced or reduced accessibility upon infection. Of the enhanced families, 15 showed high variability between individuals and had distinct epigenetic profiles. Motif analysis showed an association with known immune regulators (e.g., BATFs, FOSs/JUNs, IRFs, STATs, NFkBs, NFYs, and RELs) in stably enriched families and with other factors in variable families, including KRAB-ZNFs. We showed that TEs and host factors regulating TEs were predictive of viral load post-infection. Our findings shed light on the role TEs and KRAB-ZNFs may play in inter-individual variation in immunity.

Native functions of short tandem repeats

Over a third of the human genome is comprised of repetitive sequences, including more than a million short tandem repeats (STRs). While studies of the pathologic consequences of repeat expansions that cause syndromic human diseases are extensive, the potential native functions of STRs are often ignored. Here, we summarize a growing body of research into the normal biological functions for repetitive elements across the genome, with a particular focus on the roles of STRs in regulating gene expression. We propose reconceptualizing the pathogenic consequences of repeat expansions as aberrancies in normal gene regulation. From this altered viewpoint, we predict that future work will reveal broader roles for STRs in neuronal function and as risk alleles for more common human neurological diseases.

ZFP462 safeguards neural lineage specification by targeting G9A/GLP-mediated heterochromatin to silence enhancers

ZNF462 haploinsufficiency is linked to Weiss-Kruszka syndrome, a genetic disorder characterized by neurodevelopmental defects, including autism. Though conserved in vertebrates and essential for embryonic development, the molecular functions of ZNF462 remain unclear. We identified its murine homologue ZFP462 in a screen for mediators of epigenetic gene silencing. Here we show that ZFP462 safeguards neural lineage specification of mouse embryonic stem cells (ESCs) by targeting the H3K9-specific histone methyltransferase complex G9A/GLP to silence meso-endodermal genes. ZFP462 binds to transposable elements that are potential enhancers harbouring pluripotency and meso-endoderm transcription factor binding sites. Recruiting G9A/GLP, ZFP462 seeds heterochromatin, restricting transcription factor binding. Loss of ZFP462 in ESCs results in increased chromatin accessibility at target sites and ectopic expression of meso-endodermal genes. Taken together, ZFP462 confers lineage and locus specificity to the broadly expressed epigenetic regulator G9A/GLP. Our results suggest that aberrant activation of lineage non-specific genes in the neuronal lineage underlies ZNF462-associated neurodevelopmental pathology.

Closing the Gap: Horizontal Transfer of Mariner Transposons between Rhus Gall Aphids and Other Insects

Horizontal transfer of transposons (HTT) is an essential source of genomic evolution in eukaryotes. The HTT dynamics are well characterized in eukaryotes, including insects; however, there is a considerable gap in knowledge about HTT regarding many eukaryotes' species. In this study, we analyzed the events of the HTT between Rhus gall aphids (Hemiptera) and other insects. We analyzed the Mariner-like transposable elements (MLEs) belonging to Rhus gall aphids for the possible HT events. The MLEs have a patchy distribution and high similarity over the entire element length with insect MLEs from different orders. We selected representative sequences from the Rhus gall MLEs and identified five events of HT between MLEs of Rhus gall aphids and other insects from five different orders. We also found multiple HTT events among the MLEs of insects from the five orders, demonstrating that these Mariner elements have been involved in recurrent HT between Rhus gall aphids and other insects. Our current study closed the knowledge gap surrounding HTT and reported the events between Rhus gall aphids and other insects for the first time. We believe that this study about HTT events will help us understand the evolution and spread of transposable elements in the genomes of Rhus gall aphids.

Paramecium epigenetics in development and proliferation

The term epigenetics is used for any layer of genetic information aside from the DNA base-sequence information. Mammalian epigenetic research increased our understanding of chromatin dynamics in terms of cytosine methylation and histone modification during differentiation, aging, and disease. Instead, ciliate epigenetics focused more on small RNA-mediated effects. On the one hand, these do concern the transport of RNA from parental to daughter nuclei, representing a regulated transfer of epigenetic information across generations. On the other hand, studies of Paramecium, Tetrahymena, Oxytricha, and Stylonychia revealed an almost unique function of transgenerational RNA. Rather than solely controlling chromatin dynamics, they control sexual progeny's DNA content quantitatively and qualitatively. Thus epigenetics seems to control genetics, at least genetics of the vegetative macronucleus. This combination offers ciliates, in particular, an epigenetically controlled genetic variability. This review summarizes the epigenetic mechanisms that contribute to macronuclear heterogeneity and relates these to nuclear dimorphism. This system's adaptive and evolutionary possibilities raise the critical question of whether such a system is limited to unicellular organisms or binuclear cells. We discuss here the relevance of ciliate genetics and epigenetics to multicellular organisms.

Migrators within migrators: exploring transposable element dynamics in the monarch butterfly, Danaus plexippus

BACKGROUND

Lepidoptera (butterflies and moths) are an important model system in ecology and evolution. A high-quality chromosomal genome assembly is available for the monarch butterfly (Danaus plexippus), but it lacks an in-depth transposable element (TE) annotation, presenting an opportunity to explore monarch TE dynamics and the impact of TEs on shaping the monarch genome.

RESULTS

We find 6.21% of the monarch genome is comprised of TEs, a reduction of 6.85% compared to the original TE annotation performed on the draft genome assembly. Monarch TE content is low compared to two closely related species with available genomes, Danaus chrysippus (33.97% TE) and Danaus melanippus (11.87% TE). The biggest TE contributions to genome size in the monarch are LINEs and Penelope-like elements, and three newly identified families, r2-hero_dPle (LINE), penelope-1_dPle (Penelope-like), and hase2-1_dPle (SINE), collectively contribute 34.92% of total TE content. We find evidence of recent TE activity, with two novel Tc1 families rapidly expanding over recent timescales (tc1-1_dPle, tc1-2_dPle). LINE fragments show signatures of genomic deletions indicating a high rate of TE turnover. We investigate associations between TEs and wing colouration and immune genes and identify a three-fold increase in TE content around immune genes compared to other host genes.

CONCLUSIONS

We provide a detailed TE annotation and analysis for the monarch genome, revealing a considerably smaller TE contribution to genome content compared to two closely related Danaus species with available genome assemblies. We identify highly successful novel DNA TE families rapidly expanding over recent timescales, and ongoing signatures of both TE expansion and removal highlight the dynamic nature of repeat content in the monarch genome. Our findings also suggest that insect immune genes are promising candidates for future interrogation of TE-mediated host adaptation.

Comparative analysis of transposable elements provides insights into genome evolution in the genus Camelus

BACKGROUND

Transposable elements (TEs) are common features in eukaryotic genomes that are known to affect genome evolution critically and to play roles in gene regulation. Vertebrate genomes are dominated by TEs, which can reach copy numbers in the hundreds of thousands. To date, details regarding the presence and characteristics of TEs in camelid genomes have not been made available.

RESULTS

We conducted a genome-wide comparative analysis of camelid TEs, focusing on the identification of TEs and elucidation of transposition histories in four species: Camelus dromedarius, C. bactrianus, C. ferus, and Vicugna pacos. Our TE library was created using both de novo structure-based and homology-based searching strategies ( https://github.com/kacst-bioinfo-lab/TE_ideintification_pipeline ). Annotation results indicated a similar proportion of each genomes comprising TEs (35-36%). Class I LTR retrotransposons comprised 16-20% of genomes, and mostly consisted of the endogenous retroviruses (ERVs) groups ERVL, ERVL-MaLR, ERV_classI, and ERV_classII. Non-LTR elements comprised about 12% of genomes and consisted of SINEs (MIRs) and the LINE superfamilies LINE1, LINE2, L3/CR1, and RTE clades. Least represented were the Class II DNA transposons (2%), consisting of hAT-Charlie, TcMar-Tigger, and Helitron elements and comprising about 1-2% of each genome.

CONCLUSIONS

The findings of the present study revealed that the distribution of transposable elements across camelid genomes is approximately similar. This investigation presents a characterization of TE content in four camelid to contribute to developing a better understanding of camelid genome architecture and evolution.

Role of genes related to performance and reproduction of Thoroughbreds in training and breeding - A review

Thoroughbreds have been selected for speed and stamina since the 1700s. This selection resulted in structural and functional system-wide adaptations that enhanced physiological characteristics for outstanding speed of 61-71 kph (38-44 mph) between 1,000 and 3,200 m (5 furlongs - 2 miles). At present, horseracing is still an economically important industrial sector, therefore intensive research is underway to explore genes that allow the utilisation of genetic abilities and are significant in breeding and training. This study aims to provide an overview of genetic research and its applicability related to Thoroughbreds.

Transposable Element Dynamics and Regulation during Zygotic Genome Activation in Mammalian Embryos and Embryonic Stem Cell Model Systems

Transposable elements (TEs) are mobile genetic sequences capable of duplicating and reintegrating at new regions within the genome. A growing body of evidence has demonstrated that these elements play important roles in host genome evolution, despite being traditionally viewed as parasitic elements. To prevent ectopic activation of TE transposition and transcription, they are epigenetically silenced in most somatic tissues. Intriguingly, a specific class of TEs-retrotransposons-is transiently expressed at discrete phases during mammalian development and has been linked to the establishment of totipotency during zygotic genome activation (ZGA). While mechanisms controlling TE regulation in somatic tissues have been extensively studied, the significance underlying the unique transcriptional reactivation of retrotransposons during ZGA is only beginning to be uncovered. In this review, we summarize the expression dynamics of key retrotransposons during ZGA, focusing on findings from in vivo totipotent embryos and in vitro totipotent-like embryonic stem cells (ESCs). We then dissect the functions of retrotransposons and discuss how their transcriptional activities are finetuned during early stages of mammalian development.

A Model-Driven Quantitative Analysis of Retrotransposon Distributions in the Human Genome

Retrotransposons, DNA sequences capable of creating copies of themselves, compose about half of the human genome and played a central role in the evolution of mammals. Their current position in the host genome is the result of the retrotranscription process and of the following host genome evolution. We apply a model from statistical physics to show that the genomic distribution of the two most populated classes of retrotransposons in human deviates from random placement, and that this deviation increases with time. The time dependence suggests a major role of the host genome dynamics in shaping the current retrotransposon distributions. Focusing on a neutral scenario, we show that a simple model based on random placement followed by genome expansion and sequence duplications can reproduce the empirical retrotransposon distributions, even though more complex and possibly selective mechanisms can have contributed. Besides the inherent interest in understanding the origin of current retrotransposon distributions, this work sets a general analytical framework to analyze quantitatively the effects of genome evolutionary dynamics on the distribution of genomic elements.

Horizontal transfer and subsequent explosive expansion of a DNA transposon in sea kraits (Laticauda)

Transposable elements (TEs) are self-replicating genetic sequences and are often described as important 'drivers of evolution'. This driving force is because TEs promote genomic novelty by enabling rearrangement, and through exaptation as coding and regulatory elements. However, most TE insertions potentially lead to neutral or harmful outcomes, therefore host genomes have evolved machinery to suppress TE expansion. Through horizontal transposon transfer (HTT) TEs can colonize new genomes, and since new hosts may not be able to regulate subsequent replication, these TEs may proliferate rapidly. Here, we describe HTT of the Harbinger-Snek DNA transposon into sea kraits (Laticauda), and its subsequent explosive expansion within Laticauda genomes. This HTT occurred following the divergence of Laticauda from terrestrial Australian elapids approximately 15-25 Mya. This has resulted in numerous insertions into introns and regulatory regions, with some insertions into exons which appear to have altered UTRs or added sequence to coding exons. Harbinger-Snek has rapidly expanded to make up 8-12% of Laticauda spp. genomes; this is the fastest known expansion of TEs in amniotes following HTT. Genomic changes caused by this rapid expansion may have contributed to adaptation to the amphibious-marine habitat.

The Structural, Functional and Evolutionary Impact of Transposable Elements in Eukaryotes

Transposable elements (TEs) are nearly ubiquitous in eukaryotes. The increase in genomic data, as well as progress in genome annotation and molecular biology techniques, have revealed the vast number of ways mobile elements have impacted the evolution of eukaryotes. In addition to being the main cause of difference in haploid genome size, TEs have affected the overall organization of genomes by accumulating preferentially in some genomic regions, by causing structural rearrangements or by modifying the recombination rate. Although the vast majority of insertions is neutral or deleterious, TEs have been an important source of evolutionary novelties and have played a determinant role in the evolution of fundamental biological processes. TEs have been recruited in the regulation of host genes and are implicated in the evolution of regulatory networks. They have also served as a source of protein-coding sequences or even entire genes. The impact of TEs on eukaryotic evolution is only now being fully appreciated and the role they may play in a number of biological processes, such as speciation and adaptation, remains to be deciphered.

Complete genome analysis of Vibrio mimicus strain SCCF01, a highly virulent isolate from the freshwater catfish

Vibrio mimicus is a foodborne pathogen, which is widely distributed in the aquatic environment. Moreover, it is often involved in aquatic animal diseases. In recent years, V. mimicus is an emerging pathogen in some species of Siluriformes. The strain SCCF01 was isolated from yellow catfish (Pelteobagrus fulvidraco). In this study, we aimed to perform genomic analysis of V. mimicus strain SCCF01 to identify genetic features and evolutionary relationships. Information on gene function and classification was obtained by functional annotation, and circular graph of strain SCCF01 genome, which was created by Circos v0.64. Information on virulence genes (adhesion, flagellum system, exotoxin, and secretory system, etc.) was obtained by virulence genes annotation. Genome element prediction showed that most of the mobile elements were distributed in chromosome I. Therefore, chromosome I of SCCF01 genome has more plasticity than chromosome II and might be larger in size. Genomic linear relationship between the strain of V. mimicus and strain SCCF01 was analyzed by linear pairwise comparison but was unable to determine the relationship. Gene family analysis predicted that the evolutionary direction of strain SCCF01 was: clinical strain → environmental strain → SCCF01 strain. Phylogenetic analysis showed that the strain SCCF01 was more closely related to environmental strains. According to gene family analysis and phylogenetic analysis, we speculated that strain SCCF01 has probably diverged from environmental strains.

A small targeting domain in Ty1 integrase is sufficient to direct retrotransposon integration upstream of tRNA genes

Integration of transposable elements into the genome is mutagenic. Mechanisms targeting integrations into relatively safe locations, hence minimizing deleterious consequences for cell fitness, have emerged during evolution. In budding yeast, integration of the Ty1 LTR retrotransposon upstream of RNA polymerase III (Pol III)-transcribed genes requires interaction between Ty1 integrase (IN1) and AC40, a subunit common to Pol I and Pol III. Here, we identify the Ty1 targeting domain of IN1 that ensures (i) IN1 binding to Pol I and Pol III through AC40, (ii) IN1 genome-wide recruitment to Pol I- and Pol III-transcribed genes, and (iii) Ty1 integration only at Pol III-transcribed genes, while IN1 recruitment by AC40 is insufficient to target Ty1 integration into Pol I-transcribed genes. Swapping the targeting domains between Ty5 and Ty1 integrases causes Ty5 integration at Pol III-transcribed genes, indicating that the targeting domain of IN1 alone confers Ty1 integration site specificity.

Co-option of the lineage-specific LAVA retrotransposon in the gibbon genome

Co-option of transposable elements (TEs) to become part of existing or new enhancers is an important mechanism for evolution of gene regulation. However, contributions of lineage-specific TE insertions to recent regulatory adaptations remain poorly understood. Gibbons present a suitable model to study these contributions as they have evolved a lineage-specific TE called LAVA (LINE-AluSz-VNTR-Alu _LIKE), which is still active in the gibbon genome. The LAVA retrotransposon is thought to have played a role in the emergence of the highly rearranged structure of the gibbon genome by disrupting transcription of cell cycle genes. In this study, we investigated whether LAVA may have also contributed to the evolution of gene regulation by adopting enhancer function. We characterized fixed and polymorphic LAVA insertions across multiple gibbons and found 96 LAVA elements overlapping enhancer chromatin states. Moreover, LAVA was enriched in multiple transcription factor binding motifs, was bound by an important transcription factor (PU.1), and was associated with higher levels of gene expression in cis We found gibbon-specific signatures of purifying/positive selection at 27 LAVA insertions. Two of these insertions were fixed in the gibbon lineage and overlapped with enhancer chromatin states, representing putative co-opted LAVA enhancers. These putative enhancers were located within genes encoding SETD2 and RAD9A, two proteins that facilitate accurate repair of DNA double-strand breaks and prevent chromosomal rearrangement mutations. Co-option of LAVA in these genes may have influenced regulation of processes that preserve genome integrity. Our findings highlight the importance of considering lineage-specific TEs in studying evolution of gene regulatory elements.

Insights into epigenetic patterns in mammalian early embryos

Mammalian fertilization begins with the fusion of two specialized gametes, followed by major epigenetic remodeling leading to the formation of a totipotent embryo. During the development of the pre-implantation embryo, precise reprogramming progress is a prerequisite for avoiding developmental defects or embryonic lethality, but the underlying molecular mechanisms remain elusive. For the past few years, unprecedented breakthroughs have been made in mapping the regulatory network of dynamic epigenomes during mammalian early embryo development, taking advantage of multiple advances and innovations in low-input genome-wide chromatin analysis technologies. The aim of this review is to highlight the most recent progress in understanding the mechanisms of epigenetic remodeling during early embryogenesis in mammals, including DNA methylation, histone modifications, chromatin accessibility and 3D chromatin organization.

On transposons and totipotency

Our perception of the role of the previously considered 'selfish' or 'junk' DNA has been dramatically altered in the past 20 years or so. A large proportion of this non-coding part of mammalian genomes is repetitive in nature, classified as either satellites or transposons. While repetitive elements can be termed selfish in terms of their amplification, such events have surely been co-opted by the host, suggesting by itself a likely altruistic function for the organism at the subject of such natural selection. Indeed numerous examples of transposons regulating the functional output of the host genome have been documented. Transposons provide a powerful framework for large-scale relatively rapid concerted regulatory activities with the ability to drive evolution. Mammalian totipotency has emerged as one key stage of development in which transposon-mediated regulation of gene expression has taken centre stage in the past few years. During this period, large-scale (epigenetic) reprogramming must be accomplished in order to activate the host genome. In mice and men, one particular element murine endogenous retrovirus with leucine tRNA primer (MERVL) (and its counterpart human ERVL (HERVL)) appears to have acquired roles as a key driving force in this process. Here, I will discuss and interpret the current knowledge and its implications regarding the role of transposons, particularly of long interspersed nuclear elements (LINE-1s) and endogenous retroviruses (ERVs), in the regulation of totipotency. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Transposable elements as a potent source of diverse cis-regulatory sequences in mammalian genomes

Eukaryotic gene regulation is mediated by cis-regulatory elements, which are embedded within the vast non-coding genomic space and recognized by the transcription factors in a sequence- and context-dependent manner. A large proportion of eukaryotic genomes, including at least half of the human genome, are composed of transposable elements (TEs), which in their ancestral form carried their own cis-regulatory sequences able to exploit the host trans environment to promote TE transcription and facilitate transposition. Although not all present-day TE copies have retained this regulatory function, the preexisting regulatory potential of TEs can provide a rich source of cis-regulatory innovation for the host. Here, we review recent evidence documenting diverse contributions of TE sequences to gene regulation by functioning as enhancers, promoters, silencers and boundary elements. We discuss how TE-derived enhancer sequences can rapidly facilitate changes in existing gene regulatory networks and mediate species- and cell-type-specific regulatory innovations, and we postulate a unique contribution of TEs to species-specific gene expression divergence in pluripotency and early embryogenesis. With advances in genome-wide technologies and analyses, systematic investigation of TEs' cis-regulatory potential is now possible and our understanding of the biological impact of genomic TEs is increasing. This article is part of a discussion meeting issue 'Crossroads between transposons and gene regulation'.

Decoding the Role of Satellite DNA in Genome Architecture and Plasticity-An Evolutionary and Clinical Affair

Repetitive DNA is a major organizational component of eukaryotic genomes, being intrinsically related with their architecture and evolution. Tandemly repeated satellite DNAs (satDNAs) can be found clustered in specific heterochromatin-rich chromosomal regions, building vital structures like functional centromeres and also dispersed within euchromatin. Interestingly, despite their association to critical chromosomal structures, satDNAs are widely variable among species due to their high turnover rates. This dynamic behavior has been associated with genome plasticity and chromosome rearrangements, leading to the reshaping of genomes. Here we present the current knowledge regarding satDNAs in the light of new genomic technologies, and the challenges in the study of these sequences. Furthermore, we discuss how these sequences, together with other repeats, influence genome architecture, impacting its evolution and association with disease.

The role of transposable elements activity in aging and their possible involvement in laminopathic diseases

Eukaryotic genomes contain a large number of transposable elements, part of which are still active and able to transpose in the host genome. Mobile element activation is repressed to avoid deleterious effects, such as gene mutations or chromosome rearrangements. Control of transposable elements includes a variety of mechanisms comprising silencing pathways, which are based on the production of small non-coding RNAs. Silencing can occur either through transposable element RNA degradation or through the targeting of DNA sequences by heterochromatin formation and consequent transcriptional inhibition. Since the important role of the heterochromatin silencing, the gradual loss of heterochromatin marks in constitutive heterochromatin regions during the aging process promotes derepression of transposable elements, which is considered a cause of the progressive increase in genomic instability and of the activation of inflammatory responses. This review provides an overview of the effects of heterochromatin loss on the activity of transposable elements during the aging process and the possible impact on genome function. In this context, we discuss the possible role of the nuclear lamina, a major player in heterochromatin dynamics, in the regulation of transposable element activity and potential implications in laminopathic diseases.

Expansion of a core regulon by transposable elements promotes Arabidopsis chemical diversity and pathogen defense

Plants synthesize numerous ecologically specialized, lineage-specific metabolites through biosynthetic gene duplication and functional specialization. However, it remains unclear how duplicated genes are wired into existing regulatory networks. We show that the duplicated gene CYP82C2 has been recruited into the WRKY33 regulon and indole-3-carbonylnitrile (ICN) biosynthetic pathway through exaptation of a retroduplicated LINE retrotransposon (EPCOT3) into an enhancer. The stepwise development of a chromatin-accessible WRKY33-binding site on EPCOT3 has potentiated the regulatory neofunctionalization of CYP82C2 and the evolution of inducible defense metabolite 4-hydroxy-ICN in Arabidopsis thaliana. Although transposable elements (TEs) have long been recognized to have the potential to rewire regulatory networks, these results establish a more complete understanding of how duplicated genes and TEs contribute in concert to chemical diversity and pathogen defense.

Arabidopsis plants can produce 4-hydroxyindole-3-carbonitrile (4OH-ICN) upon pathogen infection. Here, the authors show that EPCOT3, a retrotransposonderived enhancer, mediates WRKY33-binding, pathogen-responsive transcription of CYP82C2, and synthesis of 4OH-ICN.

Expansion of the fatty acyl reductase gene family shaped pheromone communication in Hymenoptera

Fatty acyl reductases (FARs) are involved in the biosynthesis of fatty alcohols that serve a range of biological roles. Insects typically harbor numerous FAR gene family members. While some FARs are involved in pheromone biosynthesis, the biological significance of the large number of FARs in insect genomes remains unclear.

Using bumble bee (Bombini) FAR expression analysis and functional characterization, hymenopteran FAR gene tree reconstruction, and inspection of transposable elements (TEs) in the genomic environment of FARs, we uncovered a massive expansion of the FAR gene family in Hymenoptera, presumably facilitated by TEs. The expansion occurred in the common ancestor of bumble bees and stingless bees (Meliponini). We found that bumble bee FARs from the expanded FAR-A ortholog group contribute to the species-specific pheromone composition. Our results indicate that expansion and functional diversification of the FAR gene family played a key role in the evolution of pheromone communication in Hymenoptera.

Ancient exapted transposable elements promote nuclear enrichment of human long noncoding RNAs

The sequence domains underlying long noncoding RNA (lncRNA) activities, including their characteristic nuclear enrichment, remain largely unknown. It has been proposed that these domains can originate from neofunctionalized fragments of transposable elements (TEs), otherwise known as RIDLs (repeat insertion domains of lncRNA), although just a handful have been identified. It is challenging to distinguish functional RIDL instances against a numerous genomic background of neutrally evolving TEs. We here show evidence that a subset of TE types experience evolutionary selection in the context of lncRNA exons. Together these comprise an enrichment group of 5374 TE fragments in 3566 loci. Their host lncRNAs tend to be functionally validated and associated with disease. This RIDL group was used to explore the relationship between TEs and lncRNA subcellular localization. By using global localization data from 10 human cell lines, we uncover a dose-dependent relationship between nuclear/cytoplasmic distribution and evolutionarily conserved L2b, MIRb, and MIRc elements. This is observed in multiple cell types and is unaffected by confounders of transcript length or expression. Experimental validation with engineered transgenes shows that these TEs drive nuclear enrichment in a natural sequence context. Together these data reveal a role for TEs in regulating the subcellular localization of lncRNAs.

Ancient exapted transposable elements promote nuclear enrichment of human long noncoding RNAs

The sequence domains underlying long noncoding RNA (lncRNA) activities, including their characteristic nuclear enrichment, remain largely unknown. It has been proposed that these domains can originate from neofunctionalized fragments of transposable elements (TEs), otherwise known as RIDLs (repeat insertion domains of lncRNA), although just a handful have been identified. It is challenging to distinguish functional RIDL instances against a numerous genomic background of neutrally evolving TEs. We here show evidence that a subset of TE types experience evolutionary selection in the context of lncRNA exons. Together these comprise an enrichment group of 5374 TE fragments in 3566 loci. Their host lncRNAs tend to be functionally validated and associated with disease. This RIDL group was used to explore the relationship between TEs and lncRNA subcellular localization. By using global localization data from 10 human cell lines, we uncover a dose-dependent relationship between nuclear/cytoplasmic distribution and evolutionarily conserved L2b, MIRb, and MIRc elements. This is observed in multiple cell types and is unaffected by confounders of transcript length or expression. Experimental validation with engineered transgenes shows that these TEs drive nuclear enrichment in a natural sequence context. Together these data reveal a role for TEs in regulating the subcellular localization of lncRNAs.

piRNAs and PIWI proteins: regulators of gene expression in development and stem cells

PIWI proteins and Piwi-interacting RNAs (piRNAs) have established and conserved roles in repressing transposable elements (TEs) in the germline of animals. However, in several biological contexts, a large proportion of piRNAs are not related to TE sequences and, accordingly, functions for piRNAs and PIWI proteins that are independent of TE regulation have been identified. This aspect of piRNA biology is expanding rapidly. Indeed, recent reports have revealed the role of piRNAs in the regulation of endogenous gene expression programs in germ cells, as well as in somatic tissues, challenging dogma in the piRNA field. In this Review, we focus on recent data addressing the biological and developmental functions of piRNAs, highlighting their roles in embryonic patterning, germ cell specification, stem cell biology, neuronal activity and metabolism.

A LINE1-Nucleolin Partnership Regulates Early Development and ESC Identity

Transposable elements represent nearly half of mammalian genomes and are generally described as parasites, or "junk DNA." The LINE1 retrotransposon is the most abundant class and is thought to be deleterious for cells, yet it is paradoxically highly expressed during early development. Here, we report that LINE1 plays essential roles in mouse embryonic stem cells (ESCs) and pre-implantation embryos. In ESCs, LINE1 acts as a nuclear RNA scaffold that recruits Nucleolin and Kap1/Trim28 to repress Dux, the master activator of a transcriptional program specific to the 2-cell embryo. In parallel, LINE1 RNA mediates binding of Nucleolin and Kap1 to rDNA, promoting rRNA synthesis and ESC self-renewal. In embryos, LINE1 RNA is required for Dux silencing, synthesis of rRNA, and exit from the 2-cell stage. The results reveal an essential partnership between LINE1 RNA, Nucleolin, Kap1, and peri-nucleolar chromatin in the regulation of transcription, developmental potency, and ESC self-renewal.

Transposable elements and gene expression during the evolution of amniotes

Background

Transposable elements (TEs) are primarily responsible for the DNA losses and gains in genome sequences that occur over time within and between species. TEs themselves evolve, with clade specific LTR/ERV, LINEs and SINEs responsible for the bulk of species-specific genomic features. Because TEs can contain regulatory motifs, they can be exapted as regulators of gene expression. While TE insertions can provide evolutionary novelty for the regulation of gene expression, their overall impact on the evolution of gene expression is unclear. Previous investigators have shown that tissue specific gene expression in amniotes is more similar across species than within species, supporting the existence of conserved developmental gene regulation. In order to understand how species-specific TE insertions might affect the evolution/conservation of gene expression, we have looked at the association of gene expression in six tissues with TE insertions in six representative amniote genomes.

Results

A novel bootstrapping approach has been used to minimise the conflation of effects of repeat types on gene expression. We compared the expression of orthologs containing recent TE insertions to orthologs that contained older TE insertions, and the expression of non-orthologs containing recent TE insertions to non-orthologs with older TE insertions. Both orthologs and non-orthologs showed significant differences in gene expression associated with TE insertions. TEs were found associated with species-specific changes in gene expression, and the magnitude and direction of expression changes were noteworthy. Overall, orthologs containing species-specific TEs were associated with lower gene expression, while in non-orthologs, non-species specific TEs were associated with higher gene expression. Exceptions were SINE elements in human and chicken, which had an opposite association with gene expression compared to other species.

Conclusions

Our observed species-specific associations of TEs with gene expression support a role for TEs in speciation/response to selection by species. TEs do not exhibit consistent associations with gene expression and observed associations can vary depending on the age of TE insertions. Based on these observations, it would be prudent to refrain from extrapolating these and previously reported associations to distantly related species.

Mobilization of retrotransposons as a cause of chromosomal diversification and rapid speciation: the case for the Antarctic teleost genus Trematomus

Background

The importance of transposable elements (TEs) in the genomic remodeling and chromosomal rearrangements that accompany lineage diversification in vertebrates remains the subject of debate. The major impediment to understanding the roles of TEs in genome evolution is the lack of comparative and integrative analyses on complete taxonomic groups. To help overcome this problem, we have focused on the Antarctic teleost genus Trematomus (Notothenioidei: Nototheniidae), as they experienced rapid speciation accompanied by dramatic chromosomal diversity. Here we apply a multi-strategy approach to determine the role of large-scale TE mobilization in chromosomal diversification within Trematomus species.

Results

Despite the extensive chromosomal rearrangements observed in Trematomus species, our measurements revealed strong interspecific genome size conservation. After identifying the DIRS1, Gypsy and Copia retrotransposon superfamilies in genomes of 13 nototheniid species, we evaluated their diversity, abundance (copy numbers) and chromosomal distribution. Four families of DIRS1, nine of Gypsy, and two of Copia were highly conserved in these genomes; DIRS1 being the most represented within Trematomus genomes. Fluorescence in situ hybridization mapping showed preferential accumulation of DIRS1 in centromeric and pericentromeric regions, both in Trematomus and other nototheniid species, but not in outgroups: species of the Sub-Antarctic notothenioid families Bovichtidae and Eleginopsidae, and the non-notothenioid family Percidae.

Conclusions

In contrast to the outgroups, High-Antarctic notothenioid species, including the genus Trematomus, were subjected to strong environmental stresses involving repeated bouts of warming above the freezing point of seawater and cooling to sub-zero temperatures on the Antarctic continental shelf during the past 40 millions of years (My). As a consequence of these repetitive environmental changes, including thermal shocks; a breakdown of epigenetic regulation that normally represses TE activity may have led to sequential waves of TE activation within their genomes. The predominance of DIRS1 in Trematomus species, their transposition mechanism, and their strategic location in “hot spots” of insertion on chromosomes are likely to have facilitated nonhomologous recombination, thereby increasing genomic rearrangements. The resulting centric and tandem fusions and fissions would favor the rapid lineage diversification, characteristic of the nototheniid adaptive radiation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4714-x) contains supplementary material, which is available to authorized users.

The role of transposable elements in functional evolution of amphioxus genome: the case of opsin gene family

Transposable elements (TEs) are able to jump to new locations (transposition) in the genome, usually after replication. They constitute the so-called selfish or junk DNA and take over large proportions of some genomes. Due to their ability to move around they can change the DNA landscape of genomes and are therefore a rich source of innovation in genes and gene regulation. Surge of sequence data in the past years has significantly facilitated large scale comparative studies. Cephalochordates have been regarded as a useful proxy to ancestral chordate condition partially due to the comparatively slow evolutionary rate at morphological and genomic level. In this study, we used opsin gene family from three Branchiostoma species as a window into cephalochordate genome evolution. We compared opsin complements in terms of family size, gene structure and sequence allowing us to identify gene duplication and gene loss events. Furthermore, analysis of the opsin containing genomic loci showed that they are populated by TEs. In summary, we provide evidence of the way transposable elements may have contributed to the evolution of opsin gene family and to the shaping of cephalochordate genomes in general.

Identifying co-opted transposable elements using comparative epigenomics

The human genome gives rise to different epigenomic landscapes that define each cell type and can be deregulated in disease. Recent efforts by ENCODE, the NIH Roadmap and the International Human Epigenome Consortium (IHEC) have made significant advances towards assembling reference epigenomic maps of various tissues. Notably, these projects have found that approximately 80% of human DNA was biochemically active in at least one epigenomic assay while only approximately 10% of the sequence displayed signs of purifying selection. Given that transposable elements (TEs) make up at least 50% of the human genome and can be actively transcribed or act as regulatory elements either for their own purposes or be co-opted for the benefit of their host; we are interested in exploring their overall contribution to the "functional" genome. Traditional methods used to identify functional DNA have relied on comparative genomics, conservation analysis and low throughput validation assays. To discover co-opted TEs, and distinguish them from noisy genomic elements, we argue that comparative epigenomic methods will also be important.

Comparative genomic analyses highlight the contribution of pseudogenized protein-coding genes to human lincRNAs

BACKGROUND

The regulatory roles of long intergenic noncoding RNAs (lincRNAs) in humans have been revealed through the use of advanced sequencing technology. Recently, three possible scenarios of lincRNA origins have been proposed: de novo origination from intergenic regions, duplication from other long noncoding RNAs, and pseudogenization from protein-coding genes. The first two scenarios are largely studied and supported, yet few studies focused on the evolution from pseudogenized protein-coding sequence to lincRNA. Due to the non-mutually exclusive nature of these three scenarios and the need of systematic investigation of lincRNA origination, we conducted a comparative genomics study to investigate the evolution of human lincRNAs.

RESULTS

Combining with syntenic analysis and stringent Blastn e-value cutoff, we found that the majority of lincRNAs are aligned to intergenic regions of other species. Interestingly, 193 human lincRNAs could have protein-coding orthologs in at least two of nine vertebrates. Transposable elements in these conserved regions in human genome are much less than expectation. Moreover, 19% of these lincRNAs have overlaps with or are close to pseudogenes in the human genome.

CONCLUSIONS

We suggest that a notable portion of lincRNAs could be derived from pseudogenized protein-coding genes. Furthermore, based on our computational analysis, we hypothesize that a subset of these lincRNAs could have potential to regulate their paralogs by functioning as competing endogenous RNAs. Our results provide evolutionary evidence of the relationship between human lincRNAs and protein-coding genes.

Genome-wide characterization of non-reference transposable element insertion polymorphisms reveals genetic diversity in tropical and temperate maize

Background

Maize was originally domesticated in a tropical environment but is now widely cultivated at temperate latitudes. Temperate and tropical maize populations have diverged both genotypically and phenotypically. Tropical maize lines grown in temperate environments usually exhibit delayed flowering, pollination, and seed set, which reduces their grain yield relative to temperate adapted maize lines. One potential mechanism by which temperate maize may have adapted to a new environment is novel transposable element insertions, which can influence gene regulation. Recent advances in sequencing technology have made it possible to study variation in transposon content and insertion location in large sets of maize lines.

Results

In total, 274,408 non-redundant TEs (NRTEs) were identified using resequencing data generated from 83 maize inbred lines. The locations of DNA TEs and copia-superfamily retrotransposons showed significant positive correlations with gene density and genetic recombination rates, whereas gypsy-superfamily retrotransposons showed a negative correlation with these two parameters. Compared to tropical maize, temperate maize had fewer unique NRTEs but higher insertion frequency, lower background recombination rates, and higher linkage disequilibrium, with more NRTEs close to flowering and stress-related genes in the genome. Association mapping demonstrated that the presence/absence of 48 NRTEs was associated with flowering time and that expression of neighboring genes differed between haplotypes where a NRTE was present or absent.

Conclusions

This study suggests that NRTEs may have played an important role in creating the variation in gene regulation that enabled the rapid adaptation of maize to diverse environments.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4103-x) contains supplementary material, which is available to authorized users.

Transposable Element Domestication As an Adaptation to Evolutionary Conflicts

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

Splendor and misery of adaptation, or the importance of neutral null for understanding evolution

The study of any biological features, including genomic sequences, typically revolves around the question: what is this for? However, population genetic theory, combined with the data of comparative genomics, clearly indicates that such a "pan-adaptationist" approach is a fallacy. The proper question is: how has this sequence evolved? And the proper null hypothesis posits that it is a result of neutral evolution: that is, it survives by sheer chance provided that it is not deleterious enough to be efficiently purged by purifying selection. To claim adaptation, the neutral null has to be falsified. The adaptationist fallacy can be costly, inducing biologists to relentlessly seek function where there is none.

RAG Recombinase as a Selective Pressure for Genome Evolution

The RAG recombinase is a domesticated transposable element co-opted in jawed vertebrates to drive the process of the so-called V(D)J recombination, which is the hallmark of the adaptive immune system to produce antigen receptors. RAG targets, namely, the Recombination Signal Sequences (RSS), are rather long and degenerated sequences, which highlights the ability of the recombinase to interact with a wide range of target sequences, including outside of antigen receptor loci. The recognition of such cryptic targets by the recombinase threatens genome integrity by promoting aberrant DNA recombination, as observed in lymphoid malignancies. Genomes evolution resulting from RAG acquisition is an ongoing discussion, in particular regarding the counter-selection of sequences resembling the RSS and the modifications of epigenetic regulation at these potential cryptic sites. Here, we describe a new bioinformatics tool to map potential RAG targets in all jawed vertebrates. We show that our REcombination Classifier (REC) outperforms the currently available tool and is suitable for full genomes scans from species other than human and mouse. Using the REC, we document a reduction in density of potential RAG targets at the transcription start sites of genes co-expressed with the rag genes and marked with high levels of the trimethylation of the lysine 4 of the histone 3 (H3K4me3), which correlates with the retention of functional RAG activity after the horizontal transfer.