Transposable elements as tools for reshaping the genome: it is a huge world after all!

High-resolution chromosome-level genome of Scylla paramamosain provides molecular insights into adaptive evolution in crabs

BACKGROUND

Evolutionary adaptation drives organismal adjustments to environmental pressures, exemplified in the diverse morphological and ecological adaptations seen in Decapoda crustaceans, particularly brachyuran crabs. Crabs thrive in diverse ecosystems, from coral reefs to hydrothermal vents and terrestrial habitats. Despite their ecological importance, the genetic mechanisms underpinning their developmental processes, reproductive strategies, and nutrient acquisition remain poorly understood.

RESULTS

Here, we report a comprehensive genomic analysis of the green mud crab Scylla paramamosain using ultralong sequencing technologies, achieving a high-quality chromosome-level assembly. The refined 1.21 Gb genome, with an impressive contig N50 of 11.45 Mb, offers a valuable genomic resource. The genome exhibits 33,662 protein-coding genes, enriched in various pathways related to development and environmental adaptation. Gene family analysis shows expansion in development-related pathways and contraction in metabolic pathways, indicating niche adaptations. Notably, investigation into Hox gene regulation sheds light on their role in pleopod development, with the Abd-A gene identified as a linchpin. Post-transcriptional regulation involving novel-miR1317 negatively regulates Abd-A levels. Furthermore, the potential role of fru gene in ovarian development and the identification of novel-miR35 as a regulator of Spfru2 add complexity to gene regulatory networks. Comparative functional analysis across Decapoda species reveals neo-functionalization of the elovl6 gene in the synthesis of long-chain polyunsaturated fatty acids (LC-PUFA), suggesting its importance in environmental adaptation.

CONCLUSIONS

Our findings shed light on various aspects of crab biology, including genome sequencing, assembly, and annotation, as well as gene family expansion, contraction, and regulatory mechanisms governing crucial developmental processes such as metamorphosis, reproductive strategies, and fatty acid metabolism.

Genome sequencing provides insights into the evolution of gene families encoding plant cell wall-degrading enzymes in longhorned beetles

With more than 36,000 species, the longhorned beetles (family Cerambycidae) are a mega-diverse lineage of mostly xylophagous insects, all of which are represented by the sole sequenced genome of the Asian longhorned beetle (Anoplophora glabripennis; Lamiinae). Their successful radiation has been linked to their ability to degrade plant cell wall components using a range of so-called plant cell wall-degrading enzymes (PCWDEs). Our previous analysis of larval gut transcriptomes demonstrated that cerambycid beetles horizontally acquired genes encoding PCWDEs from various microbial donors; these genes evolved through multiple duplication events to form gene families. To gain further insights into the evolution of these gene families during the Cerambycidae radiation, we assembled draft genomes for four beetle species belonging to three subfamilies using long-read nanopore sequencing. All the PCWDE-encoding genes we annotated from the corresponding larval gut transcriptomes were present in these draft genomes. We confirmed that the newly discovered horizontally acquired glycoside hydrolase family 7 (GH7), subfamily 26 of GH43 (GH43_26), and GH53 (all of which are absent from the A. glabripennis genome) were indeed encoded by these beetles' genome. Most of the PCWDE-encoding genes of bacterial origin gained introns after their transfer into the beetle genome. Altogether, we show that draft genome assemblies generated from nanopore long-reads offer meaningful information to the study of the evolution of gene families in insects. We anticipate that our data will support studies aiming to better understand the biology of the Cerambycidae and other beetles in general.

Toxicoepigenetics and Environmental Health: Challenges and Opportunities

The rapidly growing field of toxicoepigenetics seeks to understand how toxicant exposures interact with the epigenome to influence disease risk. Toxicoepigenetics is a promising field of environmental health research, as integrating epigenetics into the field of toxicology will enable a more thorough evaluation of toxicant-induced disease mechanisms as well as the elucidation of the role of the epigenome as a biomarker of exposure and disease and possible mediator of exposure effects. Likewise, toxicoepigenetics will enhance our knowledge of how environmental exposures, lifestyle factors, and diet interact to influence health. Ultimately, an understanding of how the environment impacts the epigenome to cause disease may inform risk assessment, permit noninvasive biomonitoring, and provide potential opportunities for therapeutic intervention. However, the translation of research from this exciting field into benefits for human and animal health presents several challenges and opportunities. Here, we describe four significant areas in which we see opportunity to transform the field and improve human health by reducing the disease burden caused by environmental exposures. These include (1) research into the mechanistic role for epigenetic change in environment-induced disease, (2) understanding key factors influencing vulnerability to the adverse effects of environmental exposures, (3) identifying appropriate biomarkers of environmental exposures and their associated diseases, and (4) determining whether the adverse effects of environment on the epigenome and human health are reversible through pharmacologic, dietary, or behavioral interventions. We then highlight several initiatives currently underway to address these challenges.

Characterisation of the symbionts in the Mediterranean fruit fly gut

Symbioses between bacteria and their insect hosts can range from loose associations through to obligate interdependence. While fundamental evolutionary insights have been gained from the in-depth study of obligate mutualisms, there is increasing interest in the evolutionary potential of flexible symbiotic associations between hosts and their gut microbiomes. Understanding relationships between microbes and hosts also offers the potential for exploitation for insect control. Here, we investigate the gut microbiome of a global agricultural pest, the Mediterranean fruit fly (Ceratitis capitata). We used 16S rRNA profiling to compare the gut microbiomes of laboratory and wild strains raised on different diets and from flies collected from various natural plant hosts. The results showed that medfly guts harbour a simple microbiome that is primarily determined by the larval diet. However, regardless of the laboratory diet or natural plant host on which flies were raised, Klebsiella spp. dominated medfly microbiomes and were resistant to removal by antibiotic treatment. We sequenced the genome of the dominant putative Klebsiella spp. (‘Medkleb’) isolated from the gut of the Toliman wild-type strain. Genome-wide ANI analysis placed Medkleb within the K. oxytoca / michiganensis group. Species level taxonomy for Medkleb was resolved using a mutli-locus phylogenetic approach - and molecular, sequence and phenotypic analyses all supported its identity as K. michiganensis. Medkleb has a genome size (5825435 bp) which is 1.6 standard deviations smaller than the mean genome size of free-living Klebsiella spp. Medkleb also lacks some genes involved in environmental sensing. Moreover, the Medkleb genome contains at least two recently acquired unique genomic islands as well as genes that encode pectinolytic enzymes capable of degrading plant cell walls. This may be advantageous given that the medfly diet includes unripe fruits containing high proportions of pectin. The results suggest that the medfly harbours a commensal gut bacterium that may have developed a mutualistic association with its host and provide nutritional benefits.

The Chinese mitten crab genome provides insights into adaptive plasticity and developmental regulation

The infraorder Brachyura (true or short-tailed crabs) represents a successful group of marine invertebrates yet with limited genomic resources. Here we report a chromosome-anchored reference genome and transcriptomes of the Chinese mitten crab Eriocheir sinensis, a catadromous crab and invasive species with wide environmental tolerance, strong osmoregulatory capacity and high fertility. We show the expansion of specific gene families in the crab, including F-ATPase, which enhances our knowledge on the adaptive plasticity of this successful invasive species. Our analysis of spatio-temporal transcriptomes and the genome of E. sinensis and other decapods shows that brachyurization development is associated with down-regulation of Hox genes at the megalopa stage when tail shortening occurs. A better understanding of the molecular mechanism regulating sexual development is achieved by integrated analysis of multiple omics. These genomic resources significantly expand the gene repertoire of Brachyura, and provide insights into the biology of this group, and Crustacea in general.

A genomic survey of Tc1-mariner transposons in nematodes suggests extensive horizontal transposon transfer events

The number of reports concerning horizontal transposon transfers (HTT) in metazoan species is considerably increased, alongside with the exponential growth of genomic sequence data However, our understanding of the mechanisms of such phenomenon is still at an early stage. Nematodes constitute an animal phylum successfully adapted to almost every ecosystem and for this reason could potentially contribute to spreading the genetic information through horizontal transfer. To date, few studies describe HTT of nematode retrotransposons. This is due to the lack of annotation of transposable elements in the sequenced nematode genomes, especially DNA transposons, which are acknowledged as the best horizontal travelers among mobile sequences. We have therefore started a survey of DNA transposons and their possible involvement in HTT in sequenced nematode genomes. Here, we describe 83 new Tc1/mariner elements distributed in 17 nematode species. Among them, nine families were possibly horizontally transferred between nematodes and the most diverse animal species, including ants as preferred partner of HTT. The results obtained suggest that HTT events involving nematodes Tc1/mariner elements are not uncommon, and that nematodes could have a possible role as transposon reservoir that, in turn, can be redistributed among animal genomes. Overall, this could be relevant to understand how the inter-species genetic flows shape the landscape of genetic variation of organisms inhabiting specific environmental communities.

R2 and Non-Site-Specific R2-Like Retrotransposons of the German Cockroach, Blattella germanica

The structural and functional organization of the ribosomal RNA gene cluster and the full-length R2 non-LTR retrotransposon (integrated into a specific site of 28S ribosomal RNA genes) of the German cockroach, Blattella germanica, is described. A partial sequence of the R2 retrotransposon of the cockroach Rhyparobia maderae is also analyzed. The analysis of previously published next-generation sequencing data from the B. germanica genome reveals a new type of retrotransposon closely related to R2 retrotransposons but with a random distribution in the genome. Phylogenetic analysis reveals that these newly described retrotransposons form a separate clade. It is shown that proteins corresponding to the open reading frames of newly described retrotransposons exhibit unequal structural domains. Within these retrotransposons, a recombination event is described. New mechanism of transposition activity is discussed. The essential structural features of R2 retrotransposons are conserved in cockroaches and are typical of previously described R2 retrotransposons. However, the investigation of the number and frequency of 5′-truncated R2 retrotransposon insertion variants in eight B. germanica populations suggests recent mobile element activity. It is shown that the pattern of 5′-truncated R2 retrotransposon copies can be an informative molecular genetic marker for revealing genetic distances between insect populations.

In Vitro Synthesis, Delivery, and Bioavailability of Exogenous mRNA in Gene Transfer Mediated by PiggyBac Transposition

Nowadays, nonviral gene transfer is currently of great importance for introducing exogenous genes into genomes and for ensuring that transgene expression is suitable for therapeutic and bioproduction purposes. The piggyBac transposon-based system is particularly interesting since it is easy to engineer and has a large cargo capacity, up to 100 kb. In its setup, the system requires only the piggyBac transposase protein and the transgene delineated by the two piggyBac-specific inverted terminal repeats. Usually the source of transposase is carried by a DNA plasmid. However, the principal drawback of this method is the lasting presence of the transposase, due to episomal persistence or possible integration of the transposase gene vector into the cell's genome. This can lead to genotoxic effects such as multiple genomic integration events and remobilization of the transposon vector once it has been integrated. One alternative to improve the safety of the system is to deliver the transposase as in vitro-synthesized messenger RNA in order to define a very narrow expression window during which a one-shot transposition process would occur. Issues that can be encountered when working on mRNA cell transfer are related to the quality of the synthetic mRNA, the system used to introduce mRNA into the cells and the bioavailability of the mRNA molecules. Here we describe a method to produce mRNA, verify its quality, determine which transfecting reagents can be used and how this mRNA is available to promote the transposition process in HeLa cells. Additionally, we illustrate this method in stromal mesenchymal cell lines in order to support hematopoiesis.

Lack of population differentiation patterns of previously identified putatively adaptive transposable element insertions at microgeographic scales

Background

Transposable elements (TEs) play an important role in genome function and evolution. It has been shown that TEs are a considerable source of adaptive changes in the genome of Drosophila melanogaster. Specifically, footprints of selection at the DNA level, the presence of population differentiation patterns across environmental gradients, and detailed mechanistic and fitness analyses of a few candidate adaptive TEs pointed to the role of TEs in environmental adaptation. However, whether the population differentiation patterns observed at large geographic scales can be replicated at a microgeographic scale has never been assessed before.

Results

In this work, we explored the population patterns of putatively adaptive TEs at a micro-spatial scale level. We compared the frequencies of TEs, previously identified as putatively adaptive and putatively neutral, in populations collected in opposite slopes of the Evolution Canyon at Mt. Carmel in Israel separated by 200 m on average. However, the differentiation patterns previously observed across large geographic distances (2000–2200 km) were not replicated at the microscale level of the Evolution Canyon populations.

Conclusion

TE insertions previously associated with D. melanogaster environmental adaptation at a macro scale level do not play such a role at the microscale level of the Evolution Canyon populations. However, these results do not exclude a role of TEs in microgeographic adaptation because the dataset analyzed in this work is restricted to TEs identified in a single North American strain and as such is highly biased and incomplete.

Reviewers

This article was reviewed by Eugene Koonin, Limsoon Wong and Fyodor Kondrashov.

LTR retroelements are intrinsic components of transcriptional networks in frogs

Background

LTR retroelements (LTR REs) constitute a major group of transposable elements widely distributed in eukaryotic genomes. Through their own mechanism of retrotranscription LTR REs enrich the genomic landscape by providing genetic variability, thus contributing to genome structure and organization. Nonetheless, transcriptomic activity of LTR REs still remains an obscure domain within cell, developmental, and organism biology.

Results

Here we present a first comparative analysis of LTR REs for anuran amphibians based on a full depth coverage transcriptome of the European pool frog, Pelophylax lessonae, the genome of the African clawed frog, Silurana tropicalis (release v7.1), and additional transcriptomes of S. tropicalis and Cyclorana alboguttata. We identified over 1000 copies of LTR REs from all four families (Bel/Pao, Ty1/Copia, Ty3/Gypsy, Retroviridae) in the genome of S. tropicalis and discovered transcripts of several of these elements in all RNA-seq datasets analyzed. Elements of the Ty3/Gypsy family were most active, especially Amn-san elements, which accounted for approximately 0.27% of the genome in Silurana. Some elements exhibited tissue specific expression patterns, for example Hydra1.1 and MuERV-like elements in Pelophylax. In S. tropicalis considerable transcription of LTR REs was observed during embryogenesis as soon as the embryonic genome became activated, i.e. at midblastula transition. In the course of embryonic development the spectrum of transcribed LTR REs changed; during gastrulation and neurulation MuERV-like and SnRV like retroviruses were abundantly transcribed while during organogenesis transcripts of the XEN1 retroviruses became much more active.

Conclusions

The differential expression of LTR REs during embryogenesis in concert with their tissue-specificity and the protein domains they encode are evidence for the functional roles these elements play as integrative parts of complex regulatory networks. Our results support the meanwhile widely accepted concept that retroelements are not simple “junk DNA” or “harmful genomic parasites” but essential components of the transcriptomic machinery in vertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-626) contains supplementary material, which is available to authorized users.

Genome-wide comparative analysis of 20 miniature inverted-repeat transposable element families in Brassica rapa and B. oleracea

Miniature inverted-repeat transposable elements (MITEs) are ubiquitous, non-autonomous class II transposable elements. Here, we conducted genome-wide comparative analysis of 20 MITE families in B. rapa, B. oleracea, and Arabidopsis thaliana. A total of 5894 and 6026 MITE members belonging to the 20 families were found in the whole genome pseudo-chromosome sequences of B. rapa and B. oleracea, respectively. Meanwhile, only four of the 20 families, comprising 573 members, were identified in the Arabidopsis genome, indicating that most of the families were activated in the Brassica genus after divergence from Arabidopsis. Copy numbers varied from 4 to 1459 for each MITE family, and there was up to 6-fold variation between B. rapa and B. oleracea. In particular, analysis of intact members showed that whereas eleven families were present in similar copy numbers in B. rapa and B. oleracea, nine families showed copy number variation ranging from 2- to 16-fold. Four of those families (BraSto-3, BraTo-3, 4, 5) were more abundant in B. rapa, and the other five (BraSto-1, BraSto-4, BraTo-1, 7 and BraHAT-1) were more abundant in B. oleracea. Overall, 54% and 51% of the MITEs resided in or within 2 kb of a gene in the B. rapa and B. oleracea genomes, respectively. Notably, 92 MITEs were found within the CDS of annotated genes, suggesting that MITEs might play roles in diversification of genes in the recently triplicated Brassica genome. MITE insertion polymorphism (MIP) analysis of 289 MITE members showed that 52% and 23% were polymorphic at the inter- and intra-species levels, respectively, indicating that there has been recent MITE activity in the Brassica genome. These recently activated MITE families with abundant MIP will provide useful resources for molecular breeding and identification of novel functional genes arising from MITE insertion.

Large-scale mapping of transposable element insertion sites using digital encoding of sample identity

Determining the genomic locations of transposable elements is a common experimental goal. When mapping large collections of transposon insertions, individualized amplification and sequencing is both time consuming and costly. We describe an approach in which large numbers of insertion lines can be simultaneously mapped in a single DNA sequencing reaction by using digital error-correcting codes to encode line identity in a unique set of barcoded pools.

Optimization of the piggyBac transposon using mRNA and insulators: toward a more reliable gene delivery system

Integrating and expressing stably a transgene into the cellular genome remain major challenges for gene-based therapies and for bioproduction purposes. While transposon vectors mediate efficient transgene integration, expression may be limited by epigenetic silencing, and persistent transposase expression may mediate multiple transposition cycles. Here, we evaluated the delivery of the piggyBac transposase messenger RNA combined with genetically insulated transposons to isolate the transgene from neighboring regulatory elements and stabilize expression. A comparison of piggyBac transposase expression from messenger RNA and DNA vectors was carried out in terms of expression levels, transposition efficiency, transgene expression and genotoxic effects, in order to calibrate and secure the transposition-based delivery system. Messenger RNA reduced the persistence of the transposase to a narrow window, thus decreasing side effects such as superfluous genomic DNA cleavage. Both the CTF/NF1 and the D4Z4 insulators were found to mediate more efficient expression from a few transposition events. We conclude that the use of engineered piggyBac transposase mRNA and insulated transposons offer promising ways of improving the quality of the integration process and sustaining the expression of transposon vectors.

Exogenous mRNA delivery and bioavailability in gene transfer mediated by piggyBac transposition

Background

Up to now, the different uptake pathways and the subsequent intracellular trafficking of plasmid DNA have been largely explored. By contrast, the mode of internalization and the intracellular routing of an exogenous mRNA in transfected cells are poorly investigated and remain to be elucidated. The bioavailability of internalized mRNA depends on its intracellular routing and its potential accumulation in dynamic sorting sites for storage: stress granules and processing bodies. This question is of particular significance when a secure transposon-based system able to integrate a therapeutic transgene into the genome is used. Transposon vectors usually require two components: a plasmid DNA, carrying the gene of interest, and a source of transposase allowing the integration of the transgene. The principal drawback is the lasting presence of the transposase, which could remobilize the transgene once it has been inserted. Our study focused on the pharmacokinetics of the transposition process mediated by the piggyBac transposase mRNA transfection. Exogenous mRNA internalization and trafficking were investigated towards a better apprehension and fine control of the piggyBac transposase bioavailability.

Results

The mRNA prototype designed in this study provides a very narrow expression window of transposase, which allows high efficiency transposition with no cytotoxicity. Our data reveal that exogenous transposase mRNA enters cells by clathrin and caveolae-mediated endocytosis, before finishing in late endosomes 3 h after transfection. At this point, the mRNA is dissociated from its carrier and localized in stress granules, but not in cytoplasmic processing bodies. Some weaker signals have been observed in stress granules at 18 h and 48 h without causing prolonged production of the transposase. So, we designed an mRNA that is efficiently translated with a peak of transposase production 18 h post-transfection without additional release of the molecule. This confines the integration of the transgene in a very small time window.

Conclusion

Our results shed light on processes of exogenous mRNA trafficking, which are crucial to estimate the mRNA bioavailability, and increase the biosafety of transgene integration mediated by transposition. This approach provides a new way for limiting the transgene copy in the genome and their remobilization by mRNA engineering and trafficking.

Transposase concentration controls transposition activity: myth or reality?

Deciphering the mechanisms underlying the regulation of DNA transposons might be central to understanding their function and dynamics in genomes. From results obtained under artificial experimental conditions, it has been proposed that some DNA transposons self-regulate their activity via overproduction inhibition (OPI), a mechanism by which transposition activity is down-regulated when the transposase is overconcentrated in cells. However, numerous studies have given contradictory results depending on the experimental conditions. Moreover, we do not know in which cellular compartment this phenomenon takes place, or whether transposases assemble to form dense foci when they are highly expressed in cells. In the present review, we focus on investigating the data available about eukaryotic transposons to explain the mechanisms underlying OPI. Data in the literature indicate that members of the IS630-Tc1-mariner, Hobo-Ac-Tam, and piggyBac superfamilies are able to use OPI to self-regulate their transposition activity in vivo in most eukaryotic cells, and that some of them are able to assemble so as to form higher order soluble oligomers. We also investigated the localization and behavior of GFP-fused transposases belonging to the mariner, Tc1-like, and piggyBac families, investigating their ability to aggregate in cells when they are overexpressed. Transposases are able to form dense foci when they are highly expressed. Moreover, the cellular compartments in which these foci are concentrated depend on the transposase, and on its expression. The data presented here suggest that sequestration in cytoplasmic or nucleoplasmic foci, or within the nucleoli, might protect the genome against the potentially genotoxic effects of the non-specific nuclease activities of eukaryotic transposases.

Endonuclease domain of the Drosophila melanogaster R2 non-LTR retrotransposon and related retroelements: a new model for transposition

The molecular mechanisms of the transposition of non-long terminal repeat (non-LTR) retrotransposons are not well understood; the key questions of how the 3′-ends of cDNA copies integrate and how site-specific integration occurs remain unresolved. Integration depends on properties of the endonuclease (EN) domain of retrotransposons. Using the EN domain of the Drosophila R2 retrotransposon as a model for other, closely related non-LTR retrotransposons, we investigated the EN domain and found that it resembles archaeal Holliday-junction resolvases. We suggest that these non-LTR retrotransposons are co-transcribed with the host transcript. Combined with the proposed resolvase activity of the EN domain, this model yields a novel mechanism for site-specific retrotransposition within this class of retrotransposons, with resolution proceeding via a Holliday junction intermediate.