Genomic dynamics of transposable elements in the western clawed frog (Silurana tropicalis)

The genomic ecosystem of transposable elements in maize

Transposable elements (TEs) constitute the majority of flowering plant DNA, reflecting their tremendous success in subverting, avoiding, and surviving the defenses of their host genomes to ensure their selfish replication. More than 85% of the sequence of the maize genome can be ascribed to past transposition, providing a major contribution to the structure of the genome. Evidence from individual loci has informed our understanding of how transposition has shaped the genome, and a number of individual TE insertions have been causally linked to dramatic phenotypic changes. Genome-wide analyses in maize and other taxa have frequently represented TEs as a relatively homogeneous class of fragmentary relics of past transposition, obscuring their evolutionary history and interaction with their host genome. Using an updated annotation of structurally intact TEs in the maize reference genome, we investigate the family-level dynamics of TEs in maize. Integrating a variety of data, from descriptors of individual TEs like coding capacity, expression, and methylation, as well as similar features of the sequence they inserted into, we model the relationship between attributes of the genomic environment and the survival of TE copies and families. In contrast to the wholesale relegation of all TEs to a single category of junk DNA, these differences reveal a diversity of survival strategies of TE families. Together these generate a rich ecology of the genome, with each TE family representing the evolution of a distinct ecological niche. We conclude that while the impact of transposition is highly family- and context-dependent, a family-level understanding of the ecology of TEs in the genome can refine our ability to predict the role of TEs in generating genetic and phenotypic diversity.

A Field Guide to Eukaryotic Transposable Elements

Transposable elements (TEs) are mobile DNA sequences that propagate within genomes. Through diverse invasion strategies, TEs have come to occupy a substantial fraction of nearly all eukaryotic genomes, and they represent a major source of genetic variation and novelty. Here we review the defining features of each major group of eukaryotic TEs and explore their evolutionary origins and relationships. We discuss how the unique biology of different TEs influences their propagation and distribution within and across genomes. Environmental and genetic factors acting at the level of the host species further modulate the activity, diversification, and fate of TEs, producing the dramatic variation in TE content observed across eukaryotes. We argue that cataloging TE diversity and dissecting the idiosyncratic behavior of individual elements are crucial to expanding our comprehension of their impact on the biology of genomes and the evolution of species.

Heterogeneity of synonymous substitution rates in the Xenopus frog genome

With the increasing availability of high quality genomic data, there is opportunity to deeply explore the genealogical relationships of different gene loci between closely related species. In this study, we utilized genomes of Xenopus laevis (XLA, a tetraploid species with (L) and (S) sub-genomes) and X. tropicalis (XTR, a diploid species) to investigate whether synonymous substitution rates among orthologous or homoeologous genes displayed any heterogeneity. From over 1500 orthologous/homoeologous genes collected, we calculated proportion of synonymous substitutions between genomes/sub-genomes (k) and found variation within and between chromosomes. Within most chromosomes, we identified higher k with distance from the centromere, likely attributed to higher substitution rates and recombination in these regions. Using maximum likelihood methods, we identified further evidence supporting rate heterogeneity, and estimated species divergence times and ancestral population sizes. Estimated species divergence times (XLA.L-XLA.S: ~25.5 mya; XLA-XTR: ~33.0 mya) were slightly younger compared to a past study, attributed to consideration of population size in our study. Meanwhile, we found very large estimated population size in the ancestral populations of the two species (N_A = 2.55 x 10⁶). Local hybridization and population structure, which have not yet been well elucidated in frogs, may be a contributing factor to these possible large population sizes.

Environmental Change-Dependent Inherited Epigenetic Response

Epigenetic modifications are a mechanism conveying environmental information to subsequent generations via parental germ lines. Research on epigenetic responses to environmental changes in wild mammals has been widely neglected, as well as studies that compare responses to changes in different environmental factors. Here, we focused on the transmission of DNA methylation changes to naive male offspring after paternal exposure to either diet (~40% less protein) or temperature increase (10 °C increased temperature). Because both experiments focused on the liver as the main metabolic and thermoregulation organ, we were able to decipher if epigenetic changes differed in response to different environmental changes. Reduced representation bisulfite sequencing (RRBS) revealed differentially methylated regions (DMRs) in annotated genomic regions in sons sired before (control) and after the fathers’ treatments. We detected both a highly specific epigenetic response dependent on the environmental factor that had changed that was reflected in genes involved in specific metabolic pathways, and a more general response to changes in outer stimuli reflected by epigenetic modifications in a small subset of genes shared between both responses. Our results indicated that fathers prepared their offspring for specific environmental changes by paternally inherited epigenetic modifications, suggesting a strong paternal contribution to adaptive processes.

Divergent Evolutionary Trajectories of Two Young, Homomorphic, and Closely Related Sex Chromosome Systems

There exists extraordinary variation among species in the degree and nature of sex chromosome divergence. However, much of our knowledge about sex chromosomes is based on comparisons between deeply diverged species with different ancestral sex chromosomes, making it difficult to establish how fast and why sex chromosomes acquire variable levels of divergence. To address this problem, we studied sex chromosome evolution in two species of African clawed frog (Xenopus), both of whom acquired novel systems for sex determination from a recent common ancestor, and both of whom have female (ZW/ZZ) heterogamy. Derived sex chromosomes of one species, X. laevis, have a small region of suppressed recombination that surrounds the sex determining locus, and have remained this way for millions of years. In the other species, X. borealis, a younger sex chromosome system exists on a different pair of chromosomes, but the region of suppressed recombination surrounding an unidentified sex determining gene is vast, spanning almost half of the sex chromosomes. Differences between these sex chromosome systems are also apparent in the extent of nucleotide divergence between the sex chromosomes carried by females. Our analyses also indicate that in autosomes of both of these species, recombination during oogenesis occurs more frequently and in different genomic locations than during spermatogenesis. These results demonstrate that new sex chromosomes can assume radically different evolutionary trajectories, with far-reaching genomic consequences. They also suggest that in some instances the origin of new triggers for sex determination may be coupled with rapid evolution sex chromosomes, including recombination suppression of large genomic regions.

LINE Insertion Polymorphisms are Abundant but at Low Frequencies across Populations of Anolis carolinensis

Vertebrate genomes differ considerably in size and structure. Among the features that show the most variation is the abundance of Long Interspersed Nuclear Elements (LINEs). Mammalian genomes contain 100,000s LINEs that belong to a single clade, L1, and in most species a single family is usually active at a time. In contrast, non-mammalian vertebrates (fish, amphibians and reptiles) contain multiple active families, belonging to several clades, but each of them is represented by a small number of recently inserted copies. It is unclear why vertebrate genomes harbor such drastic differences in LINE composition. To address this issue, we conducted whole genome resequencing to investigate the population genomics of LINEs across 13 genomes of the lizard Anolis carolinensis sampled from two geographically and genetically distinct populations in the Eastern Florida and the Gulf Atlantic regions of the United States. We used the Mobile Element Locator Tool to identify and genotype polymorphic insertions from five major clades of LINEs (CR1, L1, L2, RTE and R4) and the 41 subfamilies that constitute them. Across these groups we found large variation in the frequency of polymorphic insertions and the observed length distributions of these insertions, suggesting these groups vary in their activity and how frequently they successfully generate full-length, potentially active copies. Though we found an abundance of polymorphic insertions (over 45,000) most of these were observed at low frequencies and typically appeared as singletons. Site frequency spectra for most LINEs showed a significant shift toward low frequency alleles compared to the spectra observed for total genomic single nucleotide polymorphisms. Using Tajima's D, FST and the mean number of pairwise differences in LINE insertion polymorphisms, we found evidence that negative selection is acting on LINE families in a length-dependent manner, its effects being stronger in the larger Eastern Florida population. Our results suggest that a large effective population size and negative selection limit the expansion of polymorphic LINE insertions across these populations and that the probability of LINE polymorphisms reaching fixation is extremely low.

Heterochromatic histone modifications at transposons in Xenopus tropicalis embryos

Transposable elements are parasitic genomic elements that can be deleterious for host gene function and genome integrity. Heterochromatic histone modifications are involved in the repression of transposons. However, it remains unknown how these histone modifications mark different types of transposons during embryonic development. Here we document the variety of heterochromatic epigenetic signatures at parasitic elements during development in Xenopus tropicalis, using genome-wide ChIP-sequencing data and ChIP-qPCR analysis. We show that specific subsets of transposons in various families and subfamilies are marked by different combinations of the heterochromatic histone modifications H4K20me3, H3K9me2/3 and H3K27me3. Many DNA transposons are marked at the blastula stage already, whereas at retrotransposons the histone modifications generally accumulate at the gastrula stage or later. Furthermore, transposons marked by H3K9me3 and H4K20me3 are more prominent in gene deserts. Using intra-subfamily divergence as a proxy for age, we show that relatively young DNA transposons are preferentially marked by early embryonic H4K20me3 and H3K27me3. In contrast, relatively young retrotransposons are marked by increasing H3K9me3 and H4K20me3 during development, and are also linked to piRNA-sized small non-coding RNAs. Our results implicate distinct repression mechanisms that operate in a transposon-selective and developmental stage-specific fashion.

Real-time transposable element activity in individual live cells

The excision and reintegration of transposable elements (TEs) restructure their host genomes, generating cellular diversity involved in evolution, development, and the etiology of human diseases. Our current knowledge of TE behavior primarily results from bulk techniques that generate time and cell ensemble averages, but cannot capture cell-to-cell variation or local environmental and temporal variability. We have developed an experimental system based on the bacterial TE IS608 that uses fluorescent reporters to directly observe single TE excision events in individual cells in real time. We find that TE activity depends upon the TE's orientation in the genome and the amount of transposase protein in the cell. We also find that TE activity is highly variable throughout the lifetime of the cell. Upon entering stationary phase, TE activity increases in cells hereditarily predisposed to TE activity. These direct observations demonstrate that real-time live-cell imaging of evolution at the molecular and individual event level is a powerful tool for the exploration of genome plasticity in stressed cells.

Transposons, Genome Size, and Evolutionary Insights in Animals

The relationship between genome size and the percentage of transposons in 161 animal species evidenced that variations in genome size are linked to the amplification or the contraction of transposable elements. The activity of transposable elements could represent a response to environmental stressors. Indeed, although with different trends in protostomes and deuterostomes, comprehensive changes in genome size were recorded in concomitance with particular periods of evolutionary history or adaptations to specific environments. During evolution, genome size and the presence of transposable elements have influenced structural and functional parameters of genomes and cells. Changes of these parameters have had an impact on morphological and functional characteristics of the organism on which natural selection directly acts. Therefore, the current situation represents a balance between insertion and amplification of transposons and the mechanisms responsible for their deletion or for decreasing their activity. Among the latter, methylation and the silencing action of small RNAs likely represent the most frequent mechanisms.

Distribution of the T2-MITE Family Transposons in the Xenopus (Silurana) tropicalis Genome

The T2 family of miniature inverted-repeat transposable elements (T2-MITE) is a prevalent MITE family found in both Xenopus(Silurana) tropicalis and X. laevis. Some subfamilies, particularly T2-A1 and T2-C, may have originated prior to the diversification of the 2 Xenopus lineages and currently include active members in X. tropicalis, whereas another subfamily, T2-E, may have lost its transposition activity even earlier. The distribution of each T2-MITE subfamily in X. tropicalis was investigated and compared to evaluate the evolutionary dynamics of the T2-MITE subfamilies. The subfamilies showed differences in chromosomal distribution, uniformity of insertion density on scaffolds, ratios of upstream to downstream insertions with respect to genes, and their distance from genes. Among these, the T2-C subfamily was interesting because it was frequently inserted upstream and close to genes and because genes with close insertions of this subfamily showed high correlations in spatial expression patterns. This unique distribution and long-lived transposition activity may reflect a mutual relationship evolved between this subfamily and the host.

Insights on genome size evolution from a miniature inverted repeat transposon driving a satellite DNA

The genome size in eukaryotes does not correlate well with the number of genes they contain. We can observe this so-called C-value paradox in amphibian species. By analyzing an amphibian genome we asked how repetitive DNA can impact genome size and architecture. We describe here our discovery of a Tc1/mariner miniature inverted-repeat transposon family present in Xenopus frogs. These transposons named miDNA4 are unique since they contain a satellite DNA motif. We found that miDNA4 measured 331 bp, contained 25 bp long inverted terminal repeat sequences and a sequence motif of 119 bp present as a unique copy or as an array of 2-47 copies. We characterized the structure, dynamics, impact and evolution of the miDNA4 family and its satellite DNA in Xenopus frog genomes. This led us to propose a model for the evolution of these two repeated sequences and how they can synergize to increase genome size.