Comparative analysis of the Mercenaria mercenaria genome provides insights into the diversity of transposable elements and immune molecules in bivalve mollusks.

BACKGROUND: The hard clam Mercenaria mercenaria is a major marine resource along the Atlantic coasts of North America and has been introduced to other continents for resource restoration or aquaculture activities. Significant mortality events have been reported in the species throughout its native range as a result of diseases (microbial infections, leukemia) and acute environmental stress. In this context, the characterization of the hard clam genome can provide highly needed resources to enable basic (e.g., oncogenesis and cancer transmission, adaptation biology) and applied (clam stock enhancement, genomic selection) sciences. RESULTS: Using a combination of long and short-read sequencing technologies, a 1.86 Gb chromosome-level assembly of the clam genome was generated. The assembly was scaffolded into 19 chromosomes, with an N50 of 83 Mb. Genome annotation yielded 34,728 predicted protein-coding genes, markedly more than the few other members of the Venerida sequenced so far, with coding regions representing only 2% of the assembly. Indeed, more than half of the genome is composed of repeated elements, including transposable elements. Major chromosome rearrangements were detected between this assembly and another recent assembly derived from a genetically segregated clam stock. Comparative analysis of the clam genome allowed the identification of a marked diversification in immune-related proteins, particularly extensive tandem duplications and expansions in tumor necrosis factors (TNFs) and C1q domain-containing proteins, some of which were previously shown to play a role in clam interactions with infectious microbes. The study also generated a comparative repertoire highlighting the diversity and, in some instances, the specificity of LTR-retrotransposons elements, particularly Steamer elements in bivalves. CONCLUSIONS: The diversity of immune molecules in M. mercenaria may allow this species to cope with varying and complex microbial and environmental landscapes. The repertoire of transposable elements identified in this study, particularly Steamer elements, should be a prime target for the investigation of cancer cell development and transmission among bivalve mollusks.

Mollusc genomes reveal variability in patterns of LTR-retrotransposons dynamics.

BACKGROUND: The three superfamilies of Long Terminal Repeat (LTR) retrotransposons are a widespread kind of transposable element and a major factor in eukaryotic genome evolution. In metazoans, recent studies suggested that Copia LTR-retrotransposons display specific dynamic compared to the more abundant and diverse Gypsy elements. Indeed, Copia elements show a relative scarcity and the prevalence of only a few clades in specific hosts. Thus, BEL/Pao seems to be the second most abundant superfamily. However, the generality of these assumptions remains to be assessed. Therefore, we carried out the first large-scale comparative genomic analysis of LTR-retrotransposons in molluscs. The aim of this study was to analyse the diversity, copy numbers, genomic proportions and distribution of LTR-retrotransposons in a large host phylum. RESULTS: We compare nine genomes of molluscs and further added LTR-retrotransposons sequences detected in databases for 47 additional species. We identified 1709 families, which enabled us to define 31 clades. We show that clade richness was highly dependent on the considered superfamily. We found only three Copia clades, including GalEa and Hydra which appear to be widely distributed and highly dominant as they account for 96% of the characterised Copia elements. Among the seven BEL/Pao clades identified, Sparrow and Surcouf are characterised for the first time. We find no BEL or Pao elements, but the rare clades Dan and Flow are present in molluscs. Finally, we characterised 21 Gypsy clades, only five of which had been previously described, the C-clade being the most abundant one. Even if they are found in the same number of host species, Copia elements are clearly less abundant than BEL/Pao elements in copy number or genomic proportions, while Gypsy elements are always the most abundant ones whatever the parameter considered. CONCLUSIONS: Our analysis confirms the contrasting dynamics of Copia and Gypsy elements in metazoans and indicates that BEL/Pao represents the second most abundant superfamily, probably reflecting an intermediate dynamic. Altogether, the data obtained in several taxa highly suggest that these patterns can be generalised for most metazoans. Finally, we highlight the importance of using database information in complement of genome analyses when analyzing transposable element diversity.

Filamentous ascomycete genomes provide insights into Copia retrotransposon diversity in fungi.

BACKGROUND: The relative scarcity of Copia retrotransposons has been recently characterized in metazoans in comparison with the other superfamilies of LTR elements. Furthermore, Copia retrotransposons have often a particular dynamics that results in a highly predominant single clade of elements within a large host taxon, such as the GalEa-like retrotransposons in crustaceans. Taking advantage of the skyrocketing amount of genomic data available for fungi, we carried out the first large-scale comparative genomic analysis of the Copia clades in filamentous ascomycetes. RESULTS: Screening 30 completely sequenced genomes allowed us to identify more than 2500 Copia copies with conserved LTR, which are distributed in 138 families. Their characterization revealed that fungal Copia diversity is much broader than previously thought with at least 27 clades, 23 of which likely correspond to new ones. While the Copia copy number is low in most species, the two clades GalEa and FunCo1 are widely distributed and highly dominate Copia content as they both account for 80% of the detected sequences. CONCLUSIONS: In Fungi, GalEa retrotransposons are restricted to Pezizomycotina in which they can make up an outstandingly high proportion of the genome (up to 10% in Cenococcum geophilum). At last, we revealed that fungal GalEa elements structurally differ from all other Copia elements with an absence of Primer Binding Site. These elements however harbor a Conserved Hairpin Site which is probably essential for their transposition.

Genomic patterns of divergence in the early and late steps of speciation of the deep-sea vent thermophilic worms of the genus Alvinella.

BACKGROUND: The transient and fragmented nature of the deep-sea hydrothermal environment made of ridge subduction, plate collision and the emergence of new rifts is currently acting to separate of vent populations, promoting local adaptation and contributing to bursts of speciation and species specialization. The tube-dwelling worms Alvinella pompejana called the Pompeii worm and its sister species A. caudata live syntopically on the hottest part of deep-sea hydrothermal chimneys along the East Pacific Rise. They are exposed to extreme thermal and chemical gradients, which vary greatly in space and time, and thus represent ideal candidates for understanding the evolutionary mechanisms at play in the vent fauna evolution. RESULTS: We explored genomic patterns of divergence in the early and late stages of speciation of these emblematic worms using transcriptome assemblies and the first draft genome to better understand the relative role of geographic isolation and habitat preference in their genome evolution. Analyses were conducted on allopatric populations of Alvinella pompejana (early stage of separation) and between A. pompejana and its syntopic species Alvinella caudata (late stage of speciation). We first identified divergent genomic regions and targets of selection as well as their position in the genome over collections of orthologous genes and, then, described the speciation dynamics by documenting the annotation of the most divergent and/or positively selected genes involved in the isolation process. Gene mapping clearly indicated that divergent genes associated with the early stage of speciation, although accounting for nearly 30% of genes, are highly scattered in the genome without any island of divergence and not involved in gamete recognition or mito-nuclear incompatibilities. By contrast, genomes of A. pompejana and A. caudata are clearly separated with nearly all genes (96%) exhibiting high divergence. This congealing effect however seems to be linked to habitat specialization and still allows positive selection on genes involved in gamete recognition, as a possible long-duration process of species reinforcement. CONCLUSION: Our analyses highlight the non-negligible role of natural selection on both the early and late stages of speciation in the iconic thermophilic worms living on the walls of deep-sea hydrothermal chimneys. They shed light on the evolution of gene divergence during the process of speciation and species specialization over a very long period of time.

Chromosomal assembly of the flat oyster (Ostrea edulis L.) genome as a new genetic resource for aquaculture.

The European flat oyster (Ostrea edulis L.) is a native bivalve of the European coasts. Harvest of this species has declined during the last decades because of the appearance of two parasites that have led to the collapse of the stocks and the loss of the natural oyster beds. O. edulis has been the subject of numerous studies in population genetics and on the detection of the parasites Bonamia ostreae and Marteilia refringens. These studies investigated immune responses to these parasites at the molecular and cellular levels. Several genetic improvement programs have been initiated especially for parasite resistance. Within the framework of a European project (PERLE 2) that aims to produce genetic lines of O. edulis with hardiness traits (growth, survival, resistance) for the purpose of repopulating natural oyster beds in Brittany and reviving the culture of this species in the foreshore, obtaining a reference genome becomes essential as done recently in many bivalve species of aquaculture interest. Here, we present a chromosome-level genome assembly and annotation for the European flat oyster, generated by combining PacBio, Illumina, 10X linked, and Hi-C sequencing. The finished assembly is 887.2 Mb with a scaffold-N50 of 97.1 Mb scaffolded on the expected 10 pseudochromosomes. Annotation of the genome revealed the presence of 35,962 protein-coding genes. We analyzed in detail the transposable element (TE) diversity in the flat oyster genome, highlighted some specificities in tRNA and miRNA composition, and provided the first insight into the molecular response of O. edulis to M. refringens. This genome provides a reference for genomic studies on O. edulis to better understand its basic physiology and as a useful resource for genetic breeding in support of aquaculture and natural reef restoration.

Eukaryote DIRS1-like retrotransposons: an overview.

BACKGROUND: DIRS1-like elements compose one superfamily of tyrosine recombinase-encoding retrotransposons. They have been previously reported in only a few diverse eukaryote species, describing a patchy distribution, and little is known about their origin and dynamics. Recently, we have shown that these retrotransposons are common among decapods, which calls into question the distribution of DIRS1-like retrotransposons among eukaryotes. RESULTS: To determine the distribution of DIRS1-like retrotransposons, we developed a new computational tool, ReDoSt, which allows us to identify well-conserved DIRS1-like elements. By screening 274 completely sequenced genomes, we identified more than 4000 DIRS1-like copies distributed among 30 diverse species which can be clustered into roughly 300 families. While the diversity in most species appears restricted to a low copy number, a few bursts of transposition are strongly suggested in certain species, such as Danio rerio and Saccoglossus kowalevskii. CONCLUSION: In this study, we report 14 new species and 8 new higher taxa that were not previously known to harbor DIRS1-like retrotransposons. Now reported in 61 species, these elements appear widely distributed among eukaryotes, even if they remain undetected in streptophytes and mammals. Especially in unikonts, a broad range of taxa from Cnidaria to Sauropsida harbors such elements. Both the distribution and the similarities between the DIRS1-like element phylogeny and conventional phylogenies of the host species suggest that DIRS1-like retrotransposons emerged early during the radiation of eukaryotes.

Comparative genomic and transcriptomic analyses of transposable elements in polychaetous annelids highlight LTR retrotransposon diversity and evolution.

BACKGROUND: With the expansion of high throughput sequencing, we now have access to a larger number of genome-wide studies analyzing the Transposable elements (TEs) composition in a wide variety of organisms. However, genomic analyses often remain too limited in number and diversity of species investigated to study in depth the dynamics and evolutionary success of the different types of TEs among metazoans. Therefore, we chose to investigate the use of transcriptomes to describe the diversity of TEs in phylogenetically related species by conducting the first comparative analysis of TEs in two groups of polychaetes and evaluate the diversity of TEs that might impact genomic evolution as a result of their mobility. RESULTS: We present a detailed analysis of TEs distribution in transcriptomes extracted from 15 polychaetes depending on the number of reads used during assembly, and also compare these results with additional TE scans on associated low-coverage genomes. We then characterized the clades defined by 1021 LTR-retrotransposon families identified in 26 species. Clade richness was highly dependent on the considered superfamily. Copia elements appear rare and are equally distributed in only three clades, GalEa, Hydra and CoMol. Among the eight BEL/Pao clades identified in annelids, two small clades within the Sailor lineage are new for science. We characterized 17 Gypsy clades of which only 4 are new; the C-clade largely dominates with a quarter of the families. Finally, all species also expressed for the majority two distinct transcripts encoding PIWI proteins, known to be involved in control of TEs mobilities. CONCLUSIONS: This study shows that the use of transcriptomes assembled from 40 million reads was sufficient to access to the diversity and proportion of the transposable elements compared to those obtained by low coverage sequencing. Among LTR-retrotransposons Gypsy elements were unequivocally dominant but results suggest that the number of Gypsy clades, although high, may be more limited than previously thought in metazoans. For BEL/Pao elements, the organization of clades within the Sailor lineage appears more difficult to establish clearly. The Copia elements remain rare and result from the evolutionary consistent success of the same three clades.

DIRS1-like retrotransposons are widely distributed among Decapoda and are particularly present in hydrothermal vent organisms.

BACKGROUND: Transposable elements are major constituents of eukaryote genomes and have a great impact on genome structure and stability. Considering their mutational abilities, TEs can contribute to the genetic diversity and evolution of organisms. Knowledge of their distribution among several genomes is an essential condition to study their dynamics and to better understand their role in species evolution. DIRS1-like retrotransposons are a particular group of retrotransposons according to their mode of transposition that implies a tyrosine recombinase. To date, they have been described in a restricted number of species in comparison with the LTR retrotransposons. In this paper, we determine the distribution of DIRS1-like elements among 25 decapod species, 10 of them living in hydrothermal vents that correspond to particularly unstable environments. RESULTS: Using PCR approaches, we have identified 15 new DIRS1-like families in 15 diverse decapod species (shrimps, lobsters, crabs and galatheid crabs). Hydrothermal organisms show a particularly great diversity of DIRS1-like elements with 5 families characterized among Alvinocarididae shrimps and 3 in the galatheid crab Munidopsis recta. Phylogenic analyses show that these elements are divergent toward the DIRS1-like families previously described in other crustaceans and arthropods and form a new clade called AlDIRS1. At larger scale, the distribution of DIRS1-like retrotransposons appears more or less patchy depending on the taxa considered. Indeed, a scattered distribution can be observed in the infraorder Brachyura whereas all the species tested in infraorders Caridea and Astacidea harbor some DIRS1-like elements. CONCLUSION: Our results lead to nearly double both the number of DIRS1-like elements described to date, and the number of species known to harbor these ones. In this study, we provide the first degenerate primers designed to look specifically for DIRS1-like retrotransposons. They allowed for revealing for the first time a widespread distribution of these elements among a large phylum, here the order Decapoda. They also suggest some peculiar features of these retrotransposons in hydrothermal organisms where a great diversity of elements is already observed. Finally, this paper constitutes the first essential step which allows for considering further studies based on the dynamics of the DIRS1-like retrotransposons among several genomes.

LTR-retrotransposons in R. exoculata and other crustaceans: the outstanding success of GalEa-like copia elements.

Transposable elements are major constituents of eukaryote genomes and have a great impact on genome structure and stability. They can contribute to the genetic diversity and evolution of organisms. Knowledge of their distribution among several genomes is an essential condition to study their dynamics and to better understand their role in species evolution. LTR-retrotransposons have been reported in many diverse eukaryote species, describing a ubiquitous distribution. Given their abundance, diversity and their extended ranges in C-values, environment and life styles, crustaceans are a great taxon to investigate the genomic component of adaptation and its possible relationships with TEs. However, crustaceans have been greatly underrepresented in transposable element studies. Using both degenerate PCR and in silico approaches, we have identified 35 Copia and 46 Gypsy families in 15 and 18 crustacean species, respectively. In particular, we characterized several full-length elements from the shrimp Rimicaris exoculata that is listed as a model organism from hydrothermal vents. Phylogenic analyses show that Copia and Gypsy retrotransposons likely present two opposite dynamics within crustaceans. The Gypsy elements appear relatively frequent and diverse whereas Copia are much more homogeneous, as 29 of them belong to the single GalEa clade, and species- or lineage-dependent. Our results also support the hypothesis of the Copia retrotransposon scarcity in metazoans compared to Gypsy elements. In such a context, the GalEa-like elements present an outstanding wide distribution among eukaryotes, from fishes to red algae, and can be even highly predominant within a large taxon, such as Malacostraca. Their distribution among crustaceans suggests a dynamics that follows a "domino days spreading" branching process in which successive amplifications may interact positively.

Fluidity of eukaryotic genomes.

The understanding the different kinds of sequences that make up a genome, as well as their proportions in genomes (obtained by the sequencing of the complete genome), has considerably changed our idea of evolution at the genomic level. The former view of a slowly evolving genome has given way to the idea of a genome that can undergo many transformations, on a large or small scale, depending on the evolution of the different types of sequences constituting it. Here we summarise the evolution of these sequences and the impact it can have on the genome. We have focused on micro-transformations, and especially on the impact of transposable elements on genomes.

Survey of genome size in 28 hydrothermal vent species covering 10 families.

Knowledge of genome size is a useful and necessary prerequisite for the development of many genomic resources. To better understand the origins and effects of DNA gains and losses among species, it is important to collect data from a broad taxonomic base, but also from particular ecosystems. Oceanic thermal vents are an interesting model to investigate genome size in very unstable environments. Here we provide data estimated by flow cytometry for 28 vent-living species among the most representative from different hydrothermal vents. We also report the genome size of closely related coastal decapods. Haploid C-values were compared with those previously reported for species from corresponding orders or infraorders. This is the first broad survey of 2C values in vent organisms. Contrary to expectations, it shows that certain hydrothermal vent species have particularly large genomes. The vent squat lobster Munidopsis recta has the largest genome yet reported for any anomuran: 2C=31.1 pg=30.4x10(9) bp. In several groups, such as Brachyura, Phyllodocida, and Veneroida, vent species have genomes that clearly rank at the high end of published values for each group. We also describe the highest DNA content yet recorded for the Brachyura (coastal crabs Xantho pilipes and Necora puber). Finally, analysis of genome size variation across populations revealed unexpected intraspecific variation in the vent shrimp Mirocaris fortunata that could not be attributed simply to ploidy changes.

GalEa retrotransposons from galatheid squat lobsters (Decapoda, Anomura) define a new clade of Ty1/copia-like elements restricted to aquatic species.

Crustacean species have not been examined in great detail for their transposable elements content. Here we focus on galatheid crabs, which are one of the most diverse and widespread taxonomic groups of Decapoda. Ty1/copia retrotransposons are a diverse and taxonomically dispersed group. Using degenerate primers, we isolated several DNA fragments that show homology with Ty1/copia retroelements reverse transcriptase gene. We named the corresponding elements from which they originated GalEa1 to GalEa3 and analyzed one of them further by isolating various clones containing segments of GalEa1. This is the first LTR retrotransposon described in crustacean genome. Nucleotide sequencing of the clones revealed that GalEa1 has LTRs (124 bp) and that the internal sequence (4,421 bp) includes a single large ORF containing gag and pol regions. Further screening identified highly related elements in six of the nine galatheid species studied. By performing BLAST searches on genome databases, we could also identify GalEa-like elements in some fishes and Urochordata genomes. These elements define a new clade of Ty1/copia retrotransposons that differs from all other Ty1/copia elements and that seems to be restricted to aquatic species.

Distribution dynamics of the Tnt1 retrotransposon in tobacco.

Retrotransposons contribute significantly to the size, organization and genetic diversity of plant genomes. Although many retrotransposon families have been reported in plants, to this day, the tobacco Tnt1 retrotransposon remains one of the few elements for which active transposition has been shown. Demonstration that Tnt1 activation can be induced by stress has lent support to the hypothesis that, under adverse conditions, transposition can be an important source of genetic variability. Here, we compared the insertion site preference of a collection of newly transposed and pre-existing Tnt1 copies identified in plants regenerated from protoplasts or tissue culture. We find that newly transposed Tnt1 copies are targeted within or close to host gene coding sequences and that the distribution of pre-existing insertions does not vary significantly from this trend. Therefore, in spite of their potential to disrupt neighboring genes, insertions within or near CDS are not preferentially removed with age. Elimination of Tnt1 insertions within or near coding sequences may be relaxed due to the polyploid nature of the tobacco genome. Tnt1 insertions within or near CDS are thus better tolerated and can putatively contribute to the diversification of tobacco gene function.

High polymorphism of TPE repeats within natural populations of Drosophila melanogaster: a gradient of the 5TPE hobo element in Western Europe.

To find out whether the polymorphism of TPE repeats of the hobo transposable element observed in some populations results from polymorphism within flies or from variability between flies, or both, we carried out isofemale line analyses of 25 populations. We found that polymorphic populations result from the presence of polymorphic flies combined with interfly variability within these populations. The fact that populations display different levels of polymorphism, i.e., different types of element and different frequencies of polymorphic flies, can be used to differentiate between qualitatively identical populations. This showed that the geographical structuring previously observed is reinforced and, in particular, that the western European populations, which have 3TPE and 5TPE elements, display a centrifugal decrease in the frequency of 5TPE hobo elements which start in western France. This gradient supports the hypothesis of a dynamic invasion by this type of elements: a total invasion by 3TPE elements, followed by further invasions involving other types of hobo elements. Moreover, the analysis of numerous sequences in current populations revealed the existence of seven types of never-previously described hobo elements with regard to TPE repeats. This diversity, which contrasts with the conservation of other parts of the element, highlights the high mutation rate of the S region.

High mutation rate of TPE repeats: a microsatellite in the putative transposase of the hobo element in Drosophila melanogaster.

The hobo transposable element contains a polymorphic microsatellite sequence located in its coding region, the TPE repeats. Previous surveys of natural populations of Drosophila melanogaster have detected at least seven different hobo transposons. These natural populations are geographically structured with regard to TPE polymorphism, and a scenario has been proposed for the invasion process. Natural populations have recently been completely invaded by hobo elements with three TPE repeats. New elements then appeared by mutation, triggering a new stage of invasion by other elements. Since TPE polymorphism appeared over a short period of time, we focused on estimating the mutation rate of these TPE repeats. We used transgenic lines harboring three TPE and/or five TPE hobo elements that had been evolving for at least 16 generations to search for a new TPE repeat polymorphism. We detected three mutants, with four, seven, and eight TPE repeats, respectively. The estimated mutation rate of the TPE repeats is therefore higher than that of neutral microsatellites in D. melanogaster (4.2 x 10-4 versus 6.5 x 10-6). The role of the transposition mechanism and the particular structure of the TPE repeats of the hobo element in this increase in the mutation rate are discussed.

Behavior of the hobo transposable element with regard to TPE repeats in transgenic lines of Drosophila melanogaster.

The hobo transposable element of Drosophila melanogaster is known to induce a hybrid dysgenesis syndrome. Moreover it displays a polymorphism of a microsatellite in its coding region: TPE repeats. In European populations, surveys of the distribution of hobo elements with regard to TPE repeats revealed that the 5TPE element is distributed along a frequency gradient, and it is even more frequent than the 3TPE element in Western populations. This suggests that the invasive ability of the hobo elements could be related to the number of TPE repeats they contain. To test this hypothesis we monitored the evolution of 16 lines derived from five initial independent transgenic lines bearing the 3TPE element and/or the 5TPE element. Four lines bearing 5TPE elements and four bearing 3TPE elements were used as a noncompetitive genetic background to compare the evolution of the 5TPE element to that of the 3TPE element. Eight lines bearing both elements provided a competitive genetic context to study potential interactions between these two elements. We studied genetic and molecular aspects of the first 20 generations. At the molecular level, we showed that the 5TPE element is able to spread within the genome at least as efficiently as the 3TPE element. Surprisingly, at the genetic level we found that the 5TPE element is less active than the 3TPE element, and moreover may be able to regulate the activity of the 3TPE element. Our findings suggest that the invasive potential of the 5TPE element could be due not only to its intrinsic transposition capacity but also to a regulatory potential.