The transposable element profile of the anolis genome: How a lizard can provide insights into the evolution of vertebrate genome size and structure

The Satellite DNA PcH-Sat, Isolated and Characterized in the Limpet Patella caerulea (Mollusca, Gastropoda), Suggests the Origin from a Nin-SINE Transposable Element

Satellite DNA (sat-DNA) was previously described as junk and selfish DNA in the cellular economy, without a clear functional role. However, during the last two decades, evidence has been accumulated about the roles of sat-DNA in different cellular functions and its probable involvement in tumorigenesis and adaptation to environmental changes. In molluscs, studies on sat-DNAs have been performed mainly on bivalve species, especially those of economic interest. Conversely, in Gastropoda (which includes about 80% of the currently described molluscs species), studies on sat-DNA have been largely neglected. In this study, we isolated and characterized a sat-DNA, here named PcH-sat, in the limpet Patella caerulea using the restriction enzyme method, particularly HaeIII. Monomeric units of PcH-sat are 179 bp long, AT-rich (58.7%), and with an identity among monomers ranging from 91.6 to 99.8%. Southern blot showed that PcH-sat is conserved in P. depressa and P. ulyssiponensis, while a smeared signal of hybridization was present in the other three investigated limpets (P. ferruginea, P. rustica and P. vulgata). Dot blot showed that PcH-sat represents about 10% of the genome of P. caerulea, 5% of that of P. depressa, and 0.3% of that of P. ulyssiponensis. FISH showed that PcH-sat was mainly localized on pericentromeric regions of chromosome pairs 2 and 4-7 of P. caerulea (2n = 18). A database search showed that PcH-sat contains a large segment (of 118 bp) showing high identity with a homologous trait of the Nin-SINE transposable element (TE) of the patellogastropod Lottia gigantea, supporting the hypothesis that TEs are involved in the rising and tandemization processes of sat-DNAs.

Chromosome-Level Genome Assembly of the Blacktail Brush Lizard, Urosaurus nigricaudus, Reveals Dosage Compensation in an Endemic Lizard

Urosaurus nigricaudus is a phrynosomatid lizard endemic to the Baja California Peninsula in Mexico. This work presents a chromosome-level genome assembly and annotation from a male individual. We used PacBio long reads and HiRise scaffolding to generate a high-quality genomic assembly of 1.87 Gb distributed in 327 scaffolds, with an N50 of 279 Mb and an L50 of 3. Approximately 98.4% of the genome is contained in 14 scaffolds, with 6 large scaffolds (334-127 Mb) representing macrochromosomes and 8 small scaffolds (63-22 Mb) representing microchromosomes. Using standard gene modeling and transcriptomic data, we predicted 17,902 protein-coding genes on the genome. The repeat content is characterized by a large proportion of long interspersed nuclear elements that are relatively old. Synteny analysis revealed some microchromosomes with high repeat content are more prone to rearrangements but that both macro- and microchromosomes are well conserved across reptiles. We identified scaffold 14 as the X chromosome. This microchromosome presents perfect dosage compensation where the single X of males has the same expression levels as two X chromosomes in females. Finally, we estimated the effective population size for U. nigricaudus was extremely low, which may reflect a reduction in polymorphism related to it becoming a peninsular endemic.

The State of Squamate Genomics: Past, Present, and Future of Genome Research in the Most Speciose Terrestrial Vertebrate Order

Squamates include more than 11,000 extant species of lizards, snakes, and amphisbaenians, and display a dazzling diversity of phenotypes across their over 200-million-year evolutionary history on Earth. Here, we introduce and define squamates (Order Squamata) and review the history and promise of genomic investigations into the patterns and processes governing squamate evolution, given recent technological advances in DNA sequencing, genome assembly, and evolutionary analysis. We survey the most recently available whole genome assemblies for squamates, including the taxonomic distribution of available squamate genomes, and assess their quality metrics and usefulness for research. We then focus on disagreements in squamate phylogenetic inference, how methods of high-throughput phylogenomics affect these inferences, and demonstrate the promise of whole genomes to settle or sustain persistent phylogenetic arguments for squamates. We review the role transposable elements play in vertebrate evolution, methods of transposable element annotation and analysis, and further demonstrate that through the understanding of the diversity, abundance, and activity of transposable elements in squamate genomes, squamates can be an ideal model for the evolution of genome size and structure in vertebrates. We discuss how squamate genomes can contribute to other areas of biological research such as venom systems, studies of phenotypic evolution, and sex determination. Because they represent more than 30% of the living species of amniote, squamates deserve a genome consortium on par with recent efforts for other amniotes (i.e., mammals and birds) that aim to sequence most of the extant families in a clade.

Abundance of Transgene Transcript Variants Associated with Somatically Active Transgenic Helitrons from Multiple T-DNA Integration Sites in Maize

Helitrons, a novel type of mysterious DNA transposons discovered computationally prior to bench work confirmation, are components ubiquitous in most sequenced genomes of various eukaryotes, including plants, animals, and fungi. There is a paucity of empirical evidence to elucidate the mechanism of Helitrons transposition in plants. Here, by constructing several artificial defective Helitron (dHel) reporter systems, we aim to identify the autonomous Helitrons (aHel) in maize genetically and to demonstrate the transposition and repair mechanisms of Helitrons upon the dHel-GFP excision in maize. When crossing with various inbred lines, several transgenic lines produced progeny of segregated, purple-blotched kernels, resulting from a leaky expression of the C1 gene driven by the dHel-interrupted promoter. Transcription analysis indicated that the insertion of different dHels into the C1 promoter or exon would lead to multiple distinct mRNA transcripts corresponding to transgenes in the host genome. Simple excision products and circular intermediates of dHel-GFP transposition have been detected from the leaf tissue of the seedlings in F1 hybrids of transgenic lines with corresponding c1 tester, although they failed to be detected in all primary transgenic lines. These results revealed the transposition and repair mechanism of Helitrons in maize. It is strongly suggested that this reporter system can detect the genetic activity of autonomic Helitron at the molecular level. Sequence features of dHel itself, together with the flanking regions, impact the excision activity of dHel and the regulation of the dHel on the transcription level of the host gene.

Do Ty3/Gypsy Transposable Elements Play Preferential Roles in Sex Chromosome Differentiation?

Transposable elements (TEs) comprise a substantial portion of eukaryotic genomes. They have the unique ability to integrate into new locations and serve as the main source of genomic novelties by mediating chromosomal rearrangements and regulating portions of functional genes. Recent studies have revealed that TEs are abundant in sex chromosomes. In this review, we propose evolutionary relationships between specific TEs, such as Ty3/Gypsy, and sex chromosomes in different lineages based on the hypothesis that these elements contributed to sex chromosome differentiation processes. We highlight how TEs can drive the dynamics of sex-determining regions via suppression recombination under a selective force to affect the organization and structural evolution of sex chromosomes. The abundance of TEs in the sex-determining regions originates from TE-poor genomic regions, suggesting a link between TE accumulation and the emergence of the sex-determining regions. TEs are generally considered to be a hallmark of chromosome degeneration. Finally, we outline recent approaches to identify TEs and study their sex-related roles and effects in the differentiation and evolution of sex chromosomes.

A chromosome-level genome assembly for the eastern fence lizard (Sceloporus undulatus), a reptile model for physiological and evolutionary ecology

Background

High-quality genomic resources facilitate investigations into behavioral ecology, morphological and physiological adaptations, and the evolution of genomic architecture. Lizards in the genus Sceloporus have a long history as important ecological, evolutionary, and physiological models, making them a valuable target for the development of genomic resources.

Findings

We present a high-quality chromosome-level reference genome assembly, SceUnd1.0 (using 10X Genomics Chromium, HiC, and Pacific Biosciences data), and tissue/developmental stage transcriptomes for the eastern fence lizard, Sceloporus undulatus. We performed synteny analysis with other snake and lizard assemblies to identify broad patterns of chromosome evolution including the fusion of micro- and macrochromosomes. We also used this new assembly to provide improved reference-based genome assemblies for 34 additional Sceloporus species. Finally, we used RNAseq and whole-genome resequencing data to compare 3 assemblies, each representing an increased level of cost and effort: Supernova Assembly with data from 10X Genomics Chromium, HiRise Assembly that added data from HiC, and PBJelly Assembly that added data from Pacific Biosciences sequencing. We found that the Supernova Assembly contained the full genome and was a suitable reference for RNAseq and single-nucleotide polymorphism calling, but the chromosome-level scaffolds provided by the addition of HiC data allowed synteny and whole-genome association mapping analyses. The subsequent addition of PacBio data doubled the contig N50 but provided negligible gains in scaffold length.

Conclusions

These new genomic resources provide valuable tools for advanced molecular analysis of an organism that has become a model in physiology and evolutionary ecology.

Existence of Bov-B LINE Retrotransposons in Snake Lineages Reveals Recent Multiple Horizontal Gene Transfers with Copy Number Variation

Transposable elements (TEs) are dynamic elements present in all eukaryotic genomes. They can “jump” and amplify within the genome and promote segmental genome rearrangements on both autosomes and sex chromosomes by disruption of gene structures. The Bovine-B long interspersed nuclear element (Bov-B LINE) is among the most abundant TE-retrotransposon families in vertebrates due to horizontal transfer (HT) among vertebrate lineages. Recent studies have shown multiple HTs or the presence of diverse Bov-B LINE groups in the snake lineage. It is hypothesized that Bov-B LINEs are highly dynamic and that the diversity reflects multiple HTs in snake lineages. Partial sequences of Bov-B LINE from 23 snake species were characterized. Phylogenetic analysis resolved at least two Bov-B LINE groups that might correspond to henophidian and caenophidian snakes; however, the tree topology differed from that based on functional nuclear and mitochondrial gene sequences. Several Bov-B LINEs of snakes showed greater than 80% similarity to sequences obtained from insects, whereas the two Bov-B LINE groups as well as sequences from the same snake species classified in different Bov-B LINE groups showed sequence similarities of less than 80%. Calculation of estimated divergence time and pairwise divergence between all individual Bov-B LINE copies suggest invasion times ranging from 79.19 to 98.8 million years ago in snakes. Accumulation of elements in a lineage-specific fashion ranged from 9 × 10⁻⁶% to 5.63 × 10⁻²% per genome. The genomic proportion of Bov-B LINEs varied among snake species but was not directly associated with genome size or invasion time. No differentiation in Bov-B LINE copy number between males and females was observed in any of the snake species examined. Incongruence in tree topology between Bov-B LINEs and other snake phylogenies may reflect past HT events. Sequence divergence of Bov-B LINEs between copies suggests that recent multiple HTs occurred within the same evolutionary timeframe in the snake lineage. The proportion of Bov-B LINEs varies among species, reflecting species specificity in TE invasion. The rapid speciation of snakes, coinciding with Bov-B LINE invasion in snake genomes, leads us to better understand the effect of Bov-B LINEs on snake genome evolution.

Consequence of Paradigm Shift with Repeat Landscapes in Reptiles: Powerful Facilitators of Chromosomal Rearrangements for Diversity and Evolution

Reptiles are notable for the extensive genomic diversity and species richness among amniote classes, but there is nevertheless a need for detailed genome-scale studies. Although the monophyletic amniotes have recently been a focus of attention through an increasing number of genome sequencing projects, the abundant repetitive portion of the genome, termed the "repeatome", remains poorly understood across different lineages. Consisting predominantly of transposable elements or mobile and satellite sequences, these repeat elements are considered crucial in causing chromosomal rearrangements that lead to genomic diversity and evolution. Here, we propose major repeat landscapes in representative reptilian species, highlighting their evolutionary dynamics and role in mediating chromosomal rearrangements. Distinct karyotype variability, which is typically a conspicuous feature of reptile genomes, is discussed, with a particular focus on rearrangements correlated with evolutionary reorganization of micro- and macrochromosomes and sex chromosomes. The exceptional karyotype variation and extreme genomic diversity of reptiles are used to test several hypotheses concerning genomic structure, function, and evolution.

Recent Secondary Contacts, Linked Selection, and Variable Recombination Rates Shape Genomic Diversity in the Model Species Anolis carolinensis

Gaining a better understanding on how selection and neutral processes affect genomic diversity is essential to gain better insights into the mechanisms driving adaptation and speciation. However, the evolutionary processes affecting variation at a genomic scale have not been investigated in most vertebrate lineages. Here, we present the first population genomics survey using whole genome resequencing in the green anole (Anolis carolinensis). Anoles have been intensively studied to understand mechanisms underlying adaptation and speciation. The green anole in particular is an important model to study genome evolution. We quantified how demography, recombination, and selection have led to the current genetic diversity of the green anole by using whole-genome resequencing of five genetic clusters covering the entire species range. The differentiation of green anole's populations is consistent with a northward expansion from South Florida followed by genetic isolation and subsequent gene flow among adjacent genetic clusters. Dispersal out-of-Florida was accompanied by a drastic population bottleneck followed by a rapid population expansion. This event was accompanied by male-biased dispersal and/or selective sweeps on the X chromosome. We show that the interaction between linked selection and recombination is the main contributor to the genomic landscape of differentiation in the anole genome.

Accumulation of transposable elements in Hox gene clusters during adaptive radiation of Anolis lizards

Transposable elements (TEs) are DNA sequences that can insert elsewhere in the genome and modify genome structure and gene regulation. The role of TEs in evolution is contentious. One hypothesis posits that TE activity generates genomic incompatibilities that can cause reproductive isolation between incipient species. This predicts that TEs will accumulate during speciation events. Here, I tested the prediction that extant lineages with a relatively high rate of speciation have a high number of TEs in their genomes. I sequenced and analysed the TE content of a marker genomic region (Hox clusters) in Anolis lizards, a classic case of an adaptive radiation. Unlike other vertebrates, including closely related lizards, Anolis lizards have high numbers of TEs in their Hox clusters, genomic regions that regulate development of the morphological adaptations that characterize habitat specialists in these lizards. Following a burst of TE activity in the lineage leading to extant Anolis, TEs have continued to accumulate during or after speciation events, resulting in a positive relationship between TE density and lineage speciation rate. These results are consistent with the prediction that TE activity contributes to adaptive radiation by promoting speciation. Although there was no evidence that TE density per se is associated with ecological morphology, the activity of TEs in Hox clusters could have been a rich source for phenotypic variation that may have facilitated the rapid parallel morphological adaptation to microhabitats seen in extant Anolis lizards.

Multiple lineages of ancient CR1 retroposons shaped the early genome evolution of amniotes

Chicken repeat 1 (CR1) retroposons are long interspersed elements (LINEs) that are ubiquitous within amniote genomes and constitute the most abundant family of transposed elements in birds, crocodilians, turtles, and snakes. They are also present in mammalian genomes, where they reside as numerous relics of ancient retroposition events. Yet, despite their relevance for understanding amniote genome evolution, the diversity and evolution of CR1 elements has never been studied on an amniote-wide level. We reconstruct the temporal and quantitative activity of CR1 subfamilies via presence/absence analyses across crocodilian phylogeny and comparative analyses of 12 crocodilian genomes, revealing relative genomic stasis of retroposition during genome evolution of extant Crocodylia. Our large-scale phylogenetic analysis of amniote CR1 subfamilies suggests the presence of at least seven ancient CR1 lineages in the amniote ancestor; and amniote-wide analyses of CR1 successions and quantities reveal differential retention (presence of ancient relics or recent activity) of these CR1 lineages across amniote genome evolution. Interestingly, birds and lepidosaurs retained the fewest ancient CR1 lineages among amniotes and also exhibit smaller genome sizes. Our study is the first to analyze CR1 evolution in a genome-wide and amniote-wide context and the data strongly suggest that the ancestral amniote genome contained myriad CR1 elements from multiple ancient lineages, and remnants of these are still detectable in the relatively stable genomes of crocodilians and turtles. Early mammalian genome evolution was thus characterized by a drastic shift from CR1 prevalence to dominance and hyperactivity of L2 LINEs in monotremes and L1 LINEs in therians.

Early mesozoic coexistence of amniotes and hepadnaviridae

Hepadnaviridae are double-stranded DNA viruses that infect some species of birds and mammals. This includes humans, where hepatitis B viruses (HBVs) are prevalent pathogens in considerable parts of the global population. Recently, endogenized sequences of HBVs (eHBVs) have been discovered in bird genomes where they constitute direct evidence for the coexistence of these viruses and their hosts from the late Mesozoic until present. Nevertheless, virtually nothing is known about the ancient host range of this virus family in other animals. Here we report the first eHBVs from crocodilian, snake, and turtle genomes, including a turtle eHBV that endogenized >207 million years ago. This genomic “fossil” is >125 million years older than the oldest avian eHBV and provides the first direct evidence that Hepadnaviridae already existed during the Early Mesozoic. This implies that the Mesozoic fossil record of HBV infection spans three of the five major groups of land vertebrates, namely birds, crocodilians, and turtles. We show that the deep phylogenetic relationships of HBVs are largely congruent with the deep phylogeny of their amniote hosts, which suggests an ancient amniote–HBV coexistence and codivergence, at least since the Early Mesozoic. Notably, the organization of overlapping genes as well as the structure of elements involved in viral replication has remained highly conserved among HBVs along that time span, except for the presence of the X gene. We provide multiple lines of evidence that the tumor-promoting X protein of mammalian HBVs lacks a homolog in all other hepadnaviruses and propose a novel scenario for the emergence of X via segmental duplication and overprinting of pre-existing reading frames in the ancestor of mammalian HBVs. Our study reveals an unforeseen host range of prehistoric HBVs and provides novel insights into the genome evolution of hepadnaviruses throughout their long-lasting association with amniote hosts.

Viruses are not known to leave physical fossil traces, which makes our understanding of their evolutionary prehistory crucially dependent on the detection of endogenous viruses. Ancient endogenous viruses, also known as paleoviruses, are relics of viral genomes or fragments thereof that once infiltrated their host's germline and then remained as molecular “fossils” within the host genome. The massive genome sequencing of recent years has unearthed vast numbers of paleoviruses from various animal genomes, including the first endogenous hepatitis B viruses (eHBVs) in bird genomes. We screened genomes of land vertebrates (amniotes) for the presence of paleoviruses and identified ancient eHBVs in the recently sequenced genomes of crocodilians, snakes, and turtles. We report an eHBV that is >207 million years old, making it the oldest endogenous virus currently known. Furthermore, our results provide direct evidence that the Hepadnaviridae virus family infected birds, crocodilians and turtles during the Mesozoic Era, and suggest a long-lasting coexistence of these viruses and their amniote hosts at least since the Early Mesozoic. We challenge previous views on the origin of the oncogenic X gene and provide an evolutionary explanation as to why only mammalian hepatitis B infection leads to hepatocellular carcinoma.

Transposons, environmental changes, and heritable induced phenotypic variability

The mechanisms of biological evolution have always been, and still are, the subject of intense debate and modeling. One of the main problems is how the genetic variability is produced and maintained in order to make the organisms adaptable to environmental changes and therefore capable of evolving. In recent years, it has been reported that, in flies and plants, mutations in Hsp90 gene are capable to induce, with a low frequency, many different developmental abnormalities depending on the genetic backgrounds. This has suggested that the reduction of Hsp90 amount makes different development pathways more sensitive to hidden genetic variability. This suggestion revitalized a classical debate around the original Waddington hypothesis of canalization and genetic assimilation making Hsp90 the prototype of morphological capacitor. Other data have also suggested a different mechanism that revitalizes another classic debate about the response of genome to physiological and environmental stress put forward by Barbara McClintock. That data demonstrated that Hsp90 is involved in repression of transposon activity by playing a significant role in piwi-interacting RNA (piRNAs)-dependent RNA interference (RNAi) silencing. The important implication is that the fixed phenotypic abnormalities observed in Hsp90 mutants are probably related to de novo induced mutations by transposon activation. In this case, Hsp90 could be considered as a mutator. In the present theoretical paper, we discuss several possible implications about environmental stress, transposon, and evolution offering also a support to the concept of evolvability.

Lizards and LINEs: selection and demography affect the fate of L1 retrotransposons in the genome of the green anole (Anolis carolinensis)

Autonomous retrotransposons lacking long terminal repeats (LTR) account for much of the variation in genome size and structure among vertebrates. Mammalian genomes contain hundreds of thousands of non-LTR retrotransposon copies, mostly resulting from the amplification of a single clade known as L1. The genomes of teleost fish and squamate reptiles contain a much more diverse array of non-LTR retrotransposon families, whereas copy number is relatively low. The majority of non-LTR retrotransposon insertions in nonmammalian vertebrates also appear to be very recent, suggesting strong purifying selection limits the accumulation of non-LTR retrotransposon copies. It is however unclear whether this turnover model, originally proposed in Drosophila, applies to nonmammalian vertebrates. Here, we studied the population dynamics of L1 in the green anole lizard (Anolis carolinensis). We found that although most L1 elements are recent in this genome, truncated insertions accumulate readily, and many are fixed at both the population and species level. In contrast, full-length L1 insertions are found at lower population frequencies, suggesting that the turnover model only applies to longer L1 elements in Anolis. We also found that full-length L1 inserts are more likely to be fixed in populations of small effective size, suggesting that the strength of purifying selection against deleterious alleles is highly dependent on host demographic history. Similar mechanisms seem to be controlling the fate of non-LTR retrotransposons in both Anolis and teleostean fish, which suggests that mammals have considerably diverged from the ancestral vertebrate in terms of how they interact with their intragenomic parasites.

Genetic and epigenetic changes involving (retro)transposons in animal hybrids and polyploids

Transposable elements (TEs) are discrete genetic units that have the ability to change their location within chromosomal DNA, and constitute a major and rapidly evolving component of eukaryotic genomes. They can be subdivided into 2 distinct types: retrotransposons, which use an RNA intermediate for transposition, and DNA transposons, which move only as DNA. Rapid advances in genome sequencing significantly improved our understanding of TE roles in genome shaping and restructuring, and studies of transcriptomes and epigenomes shed light on the previously unknown molecular mechanisms underlying genetic and epigenetic TE controls. Knowledge of these control systems may be important for better understanding of reticulate evolution and speciation in the context of bringing different genomes together by hybridization and perturbing the established regulatory balance by ploidy changes. See also sister article focusing on plants by Bento et al. in this themed issue.

Genome reannotation of the lizard Anolis carolinensis based on 14 adult and embryonic deep transcriptomes

BACKGROUND

The green anole lizard, Anolis carolinensis, is a key species for both laboratory and field-based studies of evolutionary genetics, development, neurobiology, physiology, behavior, and ecology. As the first non-avian reptilian genome sequenced, A. carolinesis is also a prime reptilian model for comparison with other vertebrate genomes. The public databases of Ensembl and NCBI have provided a first generation gene annotation of the anole genome that relies primarily on sequence conservation with related species. A second generation annotation based on tissue-specific transcriptomes would provide a valuable resource for molecular studies.

RESULTS

Here we provide an annotation of the A. carolinensis genome based on de novo assembly of deep transcriptomes of 14 adult and embryonic tissues. This revised annotation describes 59,373 transcripts, compared to 16,533 and 18,939 currently for Ensembl and NCBI, and 22,962 predicted protein-coding genes. A key improvement in this revised annotation is coverage of untranslated region (UTR) sequences, with 79% and 59% of transcripts containing 5' and 3' UTRs, respectively. Gaps in genome sequence from the current A. carolinensis build (Anocar2.0) are highlighted by our identification of 16,542 unmapped transcripts, representing 6,695 orthologues, with less than 70% genomic coverage.

CONCLUSIONS

Incorporation of tissue-specific transcriptome sequence into the A. carolinensis genome annotation has markedly improved its utility for comparative and functional studies. Increased UTR coverage allows for more accurate predicted protein sequence and regulatory analysis. This revised annotation also provides an atlas of gene expression specific to adult and embryonic tissues.

Unexpected diversity and expression of avian endogenous retroviruses

Endogenous retroviruses (ERVs) were identified and characterized in three avian genomes to gain insight into early retroviral evolution. Using the computer program RetroTector to detect relatively intact ERVs, we identified 500 ERVs in the chicken genome, 150 in the turkey genome, and 1,200 in the zebra finch genome. Previous studies suggested that endogenous alpharetroviruses were present in chicken genomes. In this analysis, a small number of alpharetroviruses were seen in the chicken and turkey genomes; however, these were greatly outnumbered by beta-like, gamma-like, and alphabeta proviruses. While the avian ERVs belonged to the same major groups as mammalian ERVs, they were more heterogeneous. In particular, the beta-like viruses revealed an evolutionary continuum with the gradual acquisition and loss of betaretroviral markers and a transition from beta to alphabeta and then to alpharetroviruses. Thus, it appears that birds may resemble a melting pot for early ERV evolution. Many of the ERVs were integrated in clusters on chromosomes, often near centromeres. About 25% of the chicken ERVs were in or near cellular transcription units; this is nearly random. The majority of these integrations were in the sense orientation in introns. A higher-than-random number of integrations were >100 kb from the nearest gene. Deep-sequencing studies of chicken embryo fibroblasts revealed that about 20% of the 500 ERVs were transcribed and translated. A subset of these were also transcribed in vivo in chickens, showing tissue-specific patterns of expression. IMPORTANCE Studies of avian endogenous retroviruses (ERVs) have given us a glimpse of an earlier retroviral world. Three different classes of ERVs were observed with many features of mammalian retroviruses, as well as some important differences. Many avian ERVs were transcribed and translated.

New insights on the sialidase protein family revealed by a phylogenetic analysis in metazoa

Sialidases are glycohydrolytic enzymes present from virus to mammals that remove sialic acid from oligosaccharide chains. Four different sialidase forms are known in vertebrates: the lysosomal NEU1, the cytosolic NEU2 and the membrane-associated NEU3 and NEU4. These enzymes modulate the cell sialic acid content and are involved in several cellular processes and pathological conditions. Molecular defects in NEU1 are responsible for sialidosis, an inherited disease characterized by lysosomal storage disorder and neurodegeneration. The studies on the biology of sialic acids and sialyltransferases, the anabolic counterparts of sialidases, have revealed a complex picture with more than 50 sialic acid variants selectively present in the different branches of the tree of life. The gain/loss of specific sialoconjugates have been proposed as key events in the evolution of deuterostomes and Homo sapiens, as well as in the host-pathogen interactions. To date, less attention has been paid to the evolution of sialidases. Thus we have conducted a survey on the state of the sialidase family in metazoan. Using an in silico approach, we identified and characterized sialidase orthologs from 21 different organisms distributed among the evolutionary tree: Metazoa relative (Monosiga brevicollis), early Deuterostomia, precursor of Chordata and Vertebrata (teleost fishes, amphibians, reptiles, avians and early and recent mammals). We were able to reconstruct the evolution of the sialidase protein family from the ancestral sialidase NEU1 and identify a new form of the enzyme, NEU5, representing an intermediate step in the evolution leading to the modern NEU3, NEU4 and NEU2. Our study provides new insights on the mechanisms that shaped the substrate specificity and other peculiar properties of the modern mammalian sialidases. Moreover, we further confirm findings on the catalytic residues and identified enzyme loop portions that behave as rapidly diverging regions and may be involved in the evolution of specific properties of sialidases.

Transposon-derived and satellite-derived repetitive sequences play distinct functional roles in Mammalian intron size expansion

Background

Repetitive sequences (RSs) are redundant, complex at times, and often lineage-specific, representing significant “building” materials for genes and genomes. According to their origins, sequence characteristics, and ways of propagation, repetitive sequences are divided into transposable elements (TEs) and satellite sequences (SSs) as well as related subfamilies and subgroups hierarchically. The combined changes attributable to the repetitive sequences alter gene and genome architectures, such as the expansion of exonic, intronic, and intergenic sequences, and most of them propagate in a seemingly random fashion and contribute very significantly to the entire mutation spectrum of mammalian genomes.

Principal findings

Our analysis is focused on evolutional features of TEs and SSs in the intronic sequence of twelve selected mammalian genomes. We divided them into four groups—primates, large mammals, rodents, and primary mammals—and used four non-mammalian vertebrate species as the out-group. After classifying intron size variation in an intron-centric way based on RS-dominance (TE-dominant or SS-dominant intron expansions), we observed several distinct profiles in intron length and positioning in different vertebrate lineages, such as retrotransposon-dominance in mammals and DNA transposon-dominance in the lower vertebrates, amphibians and fishes. The RS patterns of mouse and rat genes are most striking, which are not only distinct from those of other mammals but also different from that of the third rodent species analyzed in this study—guinea pig. Looking into the biological functions of relevant genes, we observed a two-dimensional divergence; in particular, genes that possess SS-dominant and/or RS-free introns are enriched in tissue-specific development and transcription regulation in all mammalian lineages. In addition, we found that the tendency of transposons in increasing intron size is much stronger than that of satellites, and the combined effect of both RSs is greater than either one of them alone in a simple arithmetic sum among the mammals and the opposite is found among the four non-mammalian vertebrates.

Conclusions

TE- and SS-derived RSs represent major mutational forces shaping the size and composition of vertebrate genes and genomes, and through natural selection they either fine-tune or facilitate changes in size expansion, position variation, and duplication, and thus in functions and evolutionary paths for better survival and fitness. When analyzed globally, not only are such changes significantly diversified but also comprehensible in lineages and biological implications.