Discovery of highly divergent repeat landscapes in snake genomes using high-throughput sequencing

Genome Report: First whole genome assembly of Python regius (ball python), a model of extreme physiological and metabolic plasticity

The study of nontraditional model organisms, particularly those exhibiting extreme phenotypes, offers unique insights into adaptive mechanisms of stress response and survival. Snakes, with their remarkable physiological, metabolic, and morphological adaptations, serve as powerful models for investigating these processes. Pythons are a unique model organism that have been studied for their extreme metabolic and physiological plasticity. To date, the Burmese python (Python bivittatus) is the only member of the Pythonidae family to have been sequenced. The low contiguity of this genome and rising challenges in obtaining Burmese pythons for study prompted us to sequence, assemble, and annotate the genome of the closely related ball python (Python regius). Using a hybrid sequencing approach, we generated a 1.45 Gb genome assembly with a contig N50 greater than 18 Mb and a BUSCO score of 98%, representing the highest quality genome to date for a member of the Pythonidae family. This assembly provides a valuable resource for studying snake-specific traits and evolutionary biology. Furthermore, it enables exploration of the molecular mechanisms underlying the remarkable cardiac and muscular adaptations in pythons, such as their ability to rapidly remodel organs following feeding and resist muscular atrophy during prolonged fasting. These insights have potential applications in human health, particularly in the development of therapies targeting cardiac hypertrophy and muscular atrophy.

Chromosomal and molecular perspectives on Potamotrygon motoro (Müller & Henle, 1841) from central Amazon

Cytogenetic studies on Potamotrygon motoro (Müller & Henle, 1841) are limited to classical cytogenetic techniques, but they do reveal great karyotypic variation. The main differences are related to the karyotypic formula and the absence/presence of sex chromosome systems. Thus, this study aimed to expand knowledge of the karyotypic composition of Potamotrygon motoro from different locations of the Central Amazon using Fluorescence in situ Hybridization to investigate the distribution of ribosomal DNAs (rDNA) and microsatellites sequences (SSRs). In addition, we used the mitochondrial DNA cytochrome oxidase subunit I (mtDNA COI) to perform neighbor-joining analysis to investigate the relationships among the individuals sampled. In our study, Potamotrygon motoro presented 2n = 66 chromosomes, with 18m + 12sm + 10st + 26a and heterochromatic blocks on centromeric region of all chromosomes. The 18S rDNA is present in three chromosomal pairs and 5S rDNA is located in the pair 16, which is a feature shared among freshwater stingray species. Regarding the mapping of SSRs, dinucleotide sequences showed a greater number of sites, usually on terminal regions of chromosomal pairs, with an accumulation throughout the long arms of the pair 17. Our molecular analyses did not reveal differences between the sequences used. In general, the karyotypic differences previously reported for Potamotrygon motoro indicate the presence of different cytotypes within the species.

Comparison of the Distribution Patterns of Microsatellites Across the Genomes of Reptiles

Microsatellites or simple sequence repeats (SSRs) are prevalent across various organisms' genomes. However, their distribution patterns and evolutionary dynamics in reptile genomes are rarely studied systematically. We herein conducted a comprehensive analysis of SSRs in the genomes of 36 reptile species. Our findings revealed that the total number of SSRs ranged from 1,840,965 to 7,664,452, accounting for 2.16%-8.19% of the genomes analyzed. The relative density ranged from 21,567.82 to 81,889.41 bp per megabase (Mbp). The abundance of different SSR categories followed the pattern of imperfect SSR (I-SSR) > perfect SSR (P-SSR) > compound SSR (C-SSR). A significant positive correlation was observed between the number of SSRs and genome size (p = 0.0034), whereas SSR frequency (p = 0.013) or density (p = 0.0099) showed a negative correlation with genome size. Furthermore, no correlation was found between SSR length and genome size. Mononucleotide repeats were the most common P-SSRs in crocodilians and turtles, whereas mononucleotides, trinucleotides, or tetranucleotides were the most common P-SSRs in snakes, lizards, and tuatara. P-SSRs of varying motif sizes showed nonrandom distribution across different genic regions, with AT-rich repeats being predominant. The genomic SSR content of the squamate lineage ranked the highest in abundance and variability, whereas crocodilians and turtles showed a slowly evolving and reduced microsatellite landscape. Gene ontology enrichment and Kyoto Encyclopedia of Genes and Genomes pathway analyses indicated that genes harboring P-SSRs in the coding DNA sequence regions were primarily involved in the regulation of transcription and translation processes. The SSR dataset generated in this study provides potential candidates for functional analysis and calls for broader-scale analyses across the evolutionary spectrum.

Differential Conservation and Loss of Chicken Repeat 1 (CR1) Retrotransposons in Squamates Reveal Lineage-Specific Genome Dynamics Across Reptiles

Transposable elements (TEs) are repetitive DNA sequences which create mutations and generate genetic diversity across the tree of life. In amniote vertebrates, TEs have been mainly studied in mammals and birds, whose genomes generally display low TE diversity. Squamates (Order Squamata; including ∼11,000 extant species of lizards and snakes) show as much variation in TE abundance and activity as they do in species and phenotypes. Despite this high TE activity, squamate genomes are remarkably uniform in size. We hypothesize that novel, lineage-specific genome dynamics have evolved over the course of squamate evolution. To understand the interplay between TEs and host genomes, we analyzed the evolutionary history of the chicken repeat 1 (CR1) retrotransposon, a TE family found in most tetrapod genomes which is the dominant TE in most reptiles. We compared 113 squamate genomes to the genomes of turtles, crocodilians, and birds and used ancestral state reconstruction to identify shifts in the rate of CR1 copy number evolution across reptiles. We analyzed the repeat landscapes of CR1 in squamate genomes and determined that shifts in the rate of CR1 copy number evolution are associated with lineage-specific variation in CR1 activity. We then used phylogenetic reconstruction of CR1 subfamilies across amniotes to reveal both recent and ancient CR1 subclades across the squamate tree of life. The patterns of CR1 evolution in squamates contrast other amniotes, suggesting key differences in how TEs interact with different host genomes and at different points across evolutionary history.

Cytogenetic Analysis of Satellitome of Madagascar Leaf-Tailed Geckos

Satellite DNA (satDNA) consists of sequences of DNA that form tandem repetitions across the genome, and it is notorious for its diversity and fast evolutionary rate. Despite its importance, satDNA has been only sporadically studied in reptile lineages. Here, we sequenced genomic DNA and PCR-amplified microdissected W chromosomes on the Illumina platform in order to characterize the monomers of satDNA from the Henkel's leaf-tailed gecko U. henkeli and to compare their topology by in situ hybridization in the karyotypes of the closely related Günther's flat-tail gecko U. guentheri and gold dust day gecko P. laticauda. We identified seventeen different satDNAs; twelve of them seem to accumulate in centromeres, telomeres and/or the W chromosome. Notably, centromeric and telomeric regions seem to share similar types of satDNAs, and we found two that seem to accumulate at both edges of all chromosomes in all three species. We speculate that the long-term stability of all-acrocentric karyotypes in geckos might be explained from the presence of specific satDNAs at the centromeric regions that are strong meiotic drivers, a hypothesis that should be further tested.

Differential Conservation and Loss of CR1 Retrotransposons in Squamates Reveals Lineage-Specific Genome Dynamics across Reptiles

Transposable elements (TEs) are repetitive DNA sequences which create mutations and generate genetic diversity across the tree of life. In amniotic vertebrates, TEs have been mainly studied in mammals and birds, whose genomes generally display low TE diversity. Squamates (Order Squamata; ~11,000 extant species of lizards and snakes) show as much variation in TE abundance and activity as they do in species and phenotypes. Despite this high TE activity, squamate genomes are remarkably uniform in size. We hypothesize that novel, lineage-specific dynamics have evolved over the course of squamate evolution to constrain genome size across the order. Thus, squamates may represent a prime model for investigations into TE diversity and evolution. To understand the interplay between TEs and host genomes, we analyzed the evolutionary history of the CR1 retrotransposon, a TE family found in most tetrapod genomes. We compared 113 squamate genomes to the genomes of turtles, crocodilians, and birds, and used ancestral state reconstruction to identify shifts in the rate of CR1 copy number evolution across reptiles. We analyzed the repeat landscapes of CR1 in squamate genomes and determined that shifts in the rate of CR1 copy number evolution are associated with lineage-specific variation in CR1 activity. We then used phylogenetic reconstruction of CR1 subfamilies across amniotes to reveal both recent and ancient CR1 subclades across the squamate tree of life. The patterns of CR1 evolution in squamates contrast other amniotes, suggesting key differences in how TEs interact with different host genomes and at different points across evolutionary history.

Pattern and variation in simple sequence repeat (SSR) at different genomic regions and its implications to maize evolution and breeding

BACKGROUND

Repetitive DNA sequences accounts for over 80% of maize genome. Although simple sequence repeats (SSRs) account for only 0.03% of the genome, they have been widely used in maize genetic research and breeding as highly informative codominant DNA markers. The genome-wide distribution and polymorphism of SSRs are not well studied due to the lack of high-quality genome DNA sequence data.

RESULTS

In this study, using data from high-quality de novo-sequenced maize genomes of five representative maize inbred lines, we revealed that SSRs were more densely present in telomeric region than centromeric region, and were more abundant in genic sequences than intergenic sequences. On genic sequences, tri- and hexanucleotide motifs were more abundant in CDS sequence and some mono- and dinucleotide motifs were more abundant in UTR sequences. Median length and chromosomal density of SSRs were both narrowly range-bound, with median length of 14-18 bp and genome-wide average density of 3355.77 bp/Mbp. LTR-RTs of < 0.4 Mya had higher SSR density (4498-4992 bp/Mbp). The genome-specific and motif-specific SSR polymorphism were studied. Their potential breeding applications were discussed.

CONCLUSIONS

We found that the median length of SSR sequences of different SSR motifs was nearly constant. SSR density in genic regions was much higher than intergenic regions. In addition, SSR density at LTR-RTs of different evolutionary ages varied in a narrow range. The SSRs and their LTR-RT carriers evolved at an equal rate. All these observations indicated that SSR length and density were under control of yet unknown evolutionary forces. The chromosome region-specific and motif-specific SSR polymorphisms we observed supported the notion that SSR polymorphism was invaluable genome resource for developing highly informative genome and gene markers in maize genetic research and molecular breeding.

Genome Evolution and the Future of Phylogenomics of Non-Avian Reptiles

Non-avian reptiles comprise a large proportion of amniote vertebrate diversity, with squamate reptiles-lizards and snakes-recently overtaking birds as the most species-rich tetrapod radiation. Despite displaying an extraordinary diversity of phenotypic and genomic traits, genomic resources in non-avian reptiles have accumulated more slowly than they have in mammals and birds, the remaining amniotes. Here we review the remarkable natural history of non-avian reptiles, with a focus on the physical traits, genomic characteristics, and sequence compositional patterns that comprise key axes of variation across amniotes. We argue that the high evolutionary diversity of non-avian reptiles can fuel a new generation of whole-genome phylogenomic analyses. A survey of phylogenetic investigations in non-avian reptiles shows that sequence capture-based approaches are the most commonly used, with studies of markers known as ultraconserved elements (UCEs) especially well represented. However, many other types of markers exist and are increasingly being mined from genome assemblies in silico, including some with greater information potential than UCEs for certain investigations. We discuss the importance of high-quality genomic resources and methods for bioinformatically extracting a range of marker sets from genome assemblies. Finally, we encourage herpetologists working in genomics, genetics, evolutionary biology, and other fields to work collectively towards building genomic resources for non-avian reptiles, especially squamates, that rival those already in place for mammals and birds. Overall, the development of this cross-amniote phylogenomic tree of life will contribute to illuminate interesting dimensions of biodiversity across non-avian reptiles and broader amniotes.

Snake venom gene expression is coordinated by novel regulatory architecture and the integration of multiple co-opted vertebrate pathways

Understanding how regulatory mechanisms evolve is critical for understanding the processes that give rise to novel phenotypes. Snake venom systems represent a valuable and tractable model for testing hypotheses related to the evolution of novel regulatory networks, yet the regulatory mechanisms underlying venom production remain poorly understood. Here, we use functional genomics approaches to investigate venom regulatory architecture in the prairie rattlesnake and identify cis-regulatory sequences (enhancers and promoters), trans-regulatory transcription factors, and integrated signaling cascades involved in the regulation of snake venom genes. We find evidence that two conserved vertebrate pathways, the extracellular signal-regulated kinase and unfolded protein response pathways, were co-opted to regulate snake venom. In one large venom gene family (snake venom serine proteases), this co-option was likely facilitated by the activity of transposable elements. Patterns of snake venom gene enhancer conservation, in some cases spanning 50 million yr of lineage divergence, highlight early origins and subsequent lineage-specific adaptations that have accompanied the evolution of venom regulatory architecture. We also identify features of chromatin structure involved in venom regulation, including topologically associated domains and CTCF loops that underscore the potential importance of novel chromatin structure to coevolve when duplicated genes evolve new regulatory control. Our findings provide a model for understanding how novel regulatory systems may evolve through a combination of genomic processes, including tandem duplication of genes and regulatory sequences, cis-regulatory sequence seeding by transposable elements, and diverse transcriptional regulatory proteins controlled by a co-opted regulatory cascade.

Matamatas Chelus spp. (Testudines, Chelidae) have a remarkable evolutionary history of sex chromosomes with a long-term stable XY microchromosome system

The genus Chelus, commonly known as Matamata is one of the most emblematic and remarkable species among the Neotropical chelids. It is an Amazonian species with an extensive distribution throughout Negro/Orinoco and Amazonas River basins. Currently, two species are formally recognized: Chelus orinocensis and Chelus fimbriata and although it is still classified as "Least Concern" in the IUCN, the Matamatas are very appreciated and illegally sold in the international pet trade. Regardless, little is known regarding many aspects of its natural history. Chromosomal features for Chelus, for instance, are meagre and practically restricted to the description of the diploid number (2n = 50) for Chelus fimbriata, and its sex determining strategies are yet to be fully investigated. Here, we examined the karyotype of Chelus fimbriata and the newly described Chelus orinocensis, applying an extensive conventional and molecular cytogenetic approach. This allowed us to identify a genetic sex determining mechanism with a micro XY sex chromosome system in both species, a system that was likely present in their most common recent ancestor Chelus colombiana. Furthermore, the XY system found in Chelus orinocensis and Chelus fimbriata, as seen in other chelid species, recruited several repeat motifs, possibly prior to the split of South America and Australasian lineages, indicating that such system indeed dates back to the earliest lineages of Chelid species.

Active Expression of Genes for Protein Modification Enzymes in Habu Venom Glands

Genes encoding snake venom toxins have been studied extensively. However, genes involved in the modification and functioning of venom proteins are little known. Protobothrops is a genus of pit vipers, which are venomous and inhabit the Nansei (Southwest) islands of Japan, Taiwan China, Vietnam, Thailand, Myanmar, Nepal, Bhutan, and India. Our previous study decoded the genome of Protobothrops flavoviridis, a species endemic to the Nansei Islands, Japan, and revealed unique evolutionary processes of some venom genes. In this study, we analyzed genes that are highly expressed in venom glands to survey genes for candidate enzymes or chaperone proteins involved in toxin folding and modification. We found that, in addition to genes that encode venom proteins and ribosomal proteins, genes that encode protein disulfide isomerase (PDI) family members (orthologs of human P4HB and PDIA3), Selenoprotein M (SELENOM), and Calreticulin (CALR) are highly expressed in venom glands. Since these enzymes or chaperones are involved in protein modification and potentially possess protein folding functions, we propose that P4HB, SELENOM, CALR, and PDIA3 encode candidate enzymes or chaperones to confer toxic functions upon the venom transcriptome.

Snake Venomics: Fundamentals, Recent Updates, and a Look to the Next Decade

Venomic research, powered by techniques adapted from proteomics, transcriptomics, and genomics, seeks to unravel the diversity and complexity of venom through which knowledge can be applied in the treatment of envenoming, biodiscovery, and conservation. Snake venom proteomics is most extensively studied, but the methods varied widely, creating a massive amount of information which complicates data comparison and interpretation. Advancement in mass spectrometry technology, accompanied by growing databases and sophisticated bioinformatic tools, has overcome earlier limitations of protein identification. The progress, however, remains challenged by limited accessibility to samples, non-standardized quantitative methods, and biased interpretation of -omic data. Next-generation sequencing (NGS) technologies enable high-throughput venom-gland transcriptomics and genomics, complementing venom proteomics by providing deeper insights into the structural diversity, differential expression, regulation and functional interaction of the toxin genes. Venomic tissue sampling is, however, difficult due to strict regulations on wildlife use and transfer of biological materials in some countries. Limited resources for techniques and funding are among other pertinent issues that impede the progress of venomics, particularly in less developed regions and for neglected species. Genuine collaboration between international researchers, due recognition of regional experts by global organizations (e.g., WHO), and improved distribution of research support, should be embraced.

Horizontal Transposon Transfer and Its Implications for the Ancestral Ecology of Hydrophiine Snakes

Transposable elements (TEs), also known as jumping genes, are sequences able to move or copy themselves within a genome. As TEs move throughout genomes they often act as a source of genetic novelty, hence understanding TE evolution within lineages may help in understanding environmental adaptation. Studies into the TE content of lineages of mammals such as bats have uncovered horizontal transposon transfer (HTT) into these lineages, with squamates often also containing the same TEs. Despite the repeated finding of HTT into squamates, little comparative research has examined the evolution of TEs within squamates. Here we examine a diverse family of Australo-Melanesian snakes (Hydrophiinae) to examine if the previously identified, order-wide pattern of variable TE content and activity holds true on a smaller scale. Hydrophiinae diverged from Asian elapids ~30 Mya and have since rapidly diversified into six amphibious, ~60 marine and ~100 terrestrial species that fill a broad range of ecological niches. We find TE diversity and expansion differs between hydrophiines and their Asian relatives and identify multiple HTTs into Hydrophiinae, including three likely transferred into the ancestral hydrophiine from fish. These HTT events provide the first tangible evidence that Hydrophiinae reached Australia from Asia via a marine route.

The mitochondrial genome sequence analysis of Ophidascaris baylisi from the Burmese python (Python molurus bivittatus)

Ophidascaris species are parasitic roundworms that inhabit the python gut, resulting in severe granulomatous lesions or even death. However, the classification and nomenclature of these roundworms are still controversial. Our study aims to identify a snake roundworm from the Burmese python (Python molurus bivittatus) and analyze the mitochondrial genome. We identified this roundworm as Ophidascaris baylisi based on the morphology and cytochrome c oxidase subunit I (cox1) sequence. Ophidascaris baylisi complete mitochondrial genome was 14,784 bp in length, consisting of two non-coding regions and 36 mitochondrial genes (12 protein-coding genes, 22 tRNA genes, and two rRNA genes). The protein-coding genes used TTG, ATG, ATT, or TTA as start codons and TAG, TAA, or T as stop codons. All tRNA genes showed a TV-loop structure, except trnS1^AGN and trnS2^UCN revealed a D-loop structure. The mitochondrial large ribosomal subunit 16S (rrnL) and small ribosomal subunit 12S (rrnS) were 956 bp and 700 bp long, respectively. Phylogenetic analysis based on O. baylisi mitochondrial protein-coding genes demonstrated that O. baylisi clustered with the family Ascarididae members and was most closely related to Ophidascaris wangi. These results may enhance the nematode mitochondrial genome database and provide valuable molecular markers for further research on the taxonomy, phylogeny, and genetic relationships of Ophidascaris nematodes.

Impact of Repetitive DNA Elements on Snake Genome Biology and Evolution

The distinctive biology and unique evolutionary features of snakes make them fascinating model systems to elucidate how genomes evolve and how variation at the genomic level is interlinked with phenotypic-level evolution. Similar to other eukaryotic genomes, large proportions of snake genomes contain repetitive DNA, including transposable elements (TEs) and satellite repeats. The importance of repetitive DNA and its structural and functional role in the snake genome, remain unclear. This review highlights the major types of repeats and their proportions in snake genomes, reflecting the high diversity and composition of snake repeats. We present snakes as an emerging and important model system for the study of repetitive DNA under the impact of sex and microchromosome evolution. We assemble evidence to show that certain repetitive elements in snakes are transcriptionally active and demonstrate highly dynamic lineage-specific patterns as repeat sequences. We hypothesize that particular TEs can trigger different genomic mechanisms that might contribute to driving adaptive evolution in snakes. Finally, we review emerging approaches that may be used to study the expression of repetitive elements in complex genomes, such as snakes. The specific aspects presented here will stimulate further discussion on the role of genomic repeats in shaping snake evolution.

Tracking the recruitment and evolution of snake toxins using the evolutionary context provided by the Bothrops jararaca genome

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option.

Gigantic Genomes Provide Empirical Tests of Transposable Element Dynamics Models

Transposable elements (TEs) are a major determinant of eukaryotic genome size. The collective properties of a genomic TE community reveal the history of TE/host evolutionary dynamics and impact present-day host structure and function, from genome to organism levels. In rare cases, TE community/genome size has greatly expanded in animals, associated with increased cell size and changes to anatomy and physiology. Here, we characterize the TE landscape of the genome and transcriptome in an amphibian with a giant genome — the caecilianIchthyophis bannanicus, which we show has a genome size of 12.2 Gb. Amphibians are an important model system because the clade includes independent cases of genomic gigantism. The I. bannanicus genome differs compositionally from other giant amphibian genomes, but shares a low rate of ectopic recombination-mediated deletion. We examine TE activity using expression and divergence plots; TEs account for 15% of somatic transcription, and most superfamilies appear active. We quantify TE diversity in the caecilian, as well as other vertebrates with a range of genome sizes, using diversity indices commonly applied in community ecology. We synthesize previous models that integrate TE abundance, diversity, and activity, and test whether the caecilian meets model predictions for genomes with high TE abundance. We propose thorough, consistent characterization of TEs to strengthen future comparative analyses. Such analyses will ultimately be required to reveal whether the divergent TE assemblages found across convergent gigantic genomes reflect fundamental shared features of TE/host genome evolutionary dynamics.

Transposable Elements and Stress in Vertebrates: An Overview

Since their identification as genomic regulatory elements, Transposable Elements (TEs) were considered, at first, molecular parasites and later as an important source of genetic diversity and regulatory innovations. In vertebrates in particular, TEs have been recognized as playing an important role in major evolutionary transitions and biodiversity. Moreover, in the last decade, a significant number of papers has been published highlighting a correlation between TE activity and exposition to environmental stresses and dietary factors. In this review we present an overview of the impact of TEs in vertebrate genomes, report the silencing mechanisms adopted by host genomes to regulate TE activity, and finally we explore the effects of environmental and dietary factor exposures on TE activity in mammals, which is the most studied group among vertebrates. The studies here reported evidence that several factors can induce changes in the epigenetic status of TEs and silencing mechanisms leading to their activation with consequent effects on the host genome. The study of TE can represent a future challenge for research for developing effective markers able to detect precocious epigenetic changes and prevent human diseases.

Snake W Sex Chromosome: The Shadow of Ancestral Amniote Super-Sex Chromosome

: Heteromorphic sex chromosomes, particularly the ZZ/ZW sex chromosome system of birds and some reptiles, undergo evolutionary dynamics distinct from those of autosomes. The W sex chromosome is a unique karyological member of this heteromorphic pair, which has been extensively studied in snakes to explore the origin, evolution, and genetic diversity of amniote sex chromosomes. The snake W sex chromosome offers a fascinating model system to elucidate ancestral trajectories that have resulted in genetic divergence of amniote sex chromosomes. Although the principal mechanism driving evolution of the amniote sex chromosome remains obscure, an emerging hypothesis, supported by studies of W sex chromosomes of squamate reptiles and snakes, suggests that sex chromosomes share varied genomic blocks across several amniote lineages. This implies the possible split of an ancestral super-sex chromosome via chromosomal rearrangements. We review the major findings pertaining to sex chromosomal profiles in amniotes and discuss the evolution of an ancestral super-sex chromosome by collating recent evidence sourced mainly from the snake W sex chromosome analysis. We highlight the role of repeat-mediated sex chromosome conformation and present a genomic landscape of snake Z and W chromosomes, which reveals the relative abundance of major repeats, and identifies the expansion of certain transposable elements. The latest revolution in chromosomics, i.e., complete telomere-to-telomere assembly, offers mechanistic insights into the evolutionary origin of sex chromosomes.

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.

Genome mapping of a LYST mutation in corn snakes indicates that vertebrate chromatophore vesicles are lysosome-related organelles

Reptiles exhibit a spectacular diversity of skin colors and patterns brought about by the interactions among three chromatophore types: black melanophores with melanin-packed melanosomes, red and yellow xanthophores with pteridine- and/or carotenoid-containing vesicles, and iridophores filled with light-reflecting platelets generating structural colors. Whereas the melanosome, the only color-producing endosome in mammals and birds, has been documented as a lysosome-related organelle, the maturation paths of xanthosomes and iridosomes are unknown. Here, we first use 10x Genomics linked-reads and optical mapping to assemble and annotate a nearly chromosome-quality genome of the corn snake Pantherophis guttatus The assembly is 1.71 Gb long, with an N50 of 16.8 Mb and L50 of 24. Second, we perform mapping-by-sequencing analyses and identify a 3.9-Mb genomic interval where the lavender variant resides. The lavender color morph in corn snakes is characterized by gray, rather than red, blotches on a pink, instead of orange, background. Third, our sequencing analyses reveal a single nucleotide polymorphism introducing a premature stop codon in the lysosomal trafficking regulator gene (LYST) that shortens the corresponding protein by 603 amino acids and removes evolutionary-conserved domains. Fourth, we use light and transmission electron microscopy comparative analyses of wild type versus lavender corn snakes and show that the color-producing endosomes of all chromatophores are substantially affected in the LYST mutant. Our work provides evidence characterizing xanthosomes in xanthophores and iridosomes in iridophores as lysosome-related organelles.

The Amazonian Red Side-Necked Turtle Rhinemys rufipes (Spix, 1824) (Testudines, Chelidae) Has a GSD Sex-Determining Mechanism with an Ancient XY Sex Microchromosome System

The Amazonian red side-necked turtle Rhynemis rufipes is an endemic Amazonian Chelidae species that occurs in small streams throughout Colombia and Brazil river basins. Little is known about various biological aspects of this species, including its sex determination strategies. Among chelids, the greatest karyotype diversity is found in the Neotropical species, with several 2n configurations, including cases of triploidy. Here, we investigate the karyotype of Rhinemys rufipes by applying combined conventional and molecular cytogenetic procedures. This allowed us to discover a genetic sex-determining mechanism that shares an ancestral micro XY sex chromosome system. This ancient micro XY system recruited distinct repeat motifs before it diverged from several South America and Australasian species. We propose that such a system dates back to the earliest lineages of the chelid species before the split of South America and Australasian lineages.

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.

First Record on Sex Chromosomes in a Species of the Family Cynodontidae: Cynodon gibbus (Agassiz, 1829)

The fish family Cynodontidae belongs to the superfamily Curimatoidea, together with the Hemiodontidae, Serrasalmidae, Parodontidae, Prochilodontidae, Chilodontidae, Curimatidae, and Anostomidae. The majority of the species of this superfamily that have been analyzed to date have a diploid chromosome number of 2n = 54. Differentiated sex chromosomes (with female heterogamety) have been observed only in the Prochilodontidae, Parodontidae, and Anostomidae. The present study provides the first description of differentiated sex chromosomes in the cynodontid species Cynodon gibbus, which has a ZZ/ZW system, and shows that repetitive DNA has played a fundamental role in the differentiation of these sex chromosomes.

Draft Genome of a Blister Beetle Mylabris aulica

Mylabris aulica is a widely distributed blister beetle of the Meloidae family. It has the ability to synthesize a potent defensive secretion that includes cantharidin, a toxic compound used to treat many major illnesses. However, owing to the lack of genetic studies on cantharidin biosynthesis in M. aulica, the commercial use of this species is less extensive than that of other blister beetle species in China. This study reports a draft assembly and possible genes and pathways related to cantharidin biosynthesis for the M. aulica blister beetle using nanopore sequencing data. The draft genome assembly size was 288.5 Mb with a 467.8 Kb N50, and a repeat content of 50.62%. An integrated gene finding pipeline performed for assembly obtained 16,500 protein coding genes. Benchmarking universal single-copy orthologs assessment showed that this gene set included 94.4% complete Insecta universal single-copy orthologs. Over 99% of these genes were assigned functional annotations in the gene ontology, Kyoto Encyclopedia of Genes and Genomes, or Genbank non-redundant databases. Comparative genomic analysis showed that the completeness and continuity of our assembly was better than those of Hycleus cichorii and Hycleus phaleratus blister beetle genomes. The analysis of homologous orthologous genes and inference from evolutionary history imply that the Mylabris and Hycleus genera are genetically close, have a similar genetic background, and have differentiated within one million years. This M. aulica genome assembly provides a valuable resource for future blister beetle studies and will contribute to cantharidin biosynthesis.

Omics Technologies for Profiling Toxin Diversity and Evolution in Snake Venom: Impacts on the Discovery of Therapeutic and Diagnostic Agents

Snake venoms are primarily composed of proteins and peptides, and these toxins have developed high selectivity to their biological targets. This makes venoms interesting for exploration into protein evolution and structure-function relationships. A single venom protein superfamily can exhibit a variety of pharmacological effects; these variations in activity originate from differences in functional sites, domains, posttranslational modifications, and the formations of toxin complexes. In this review, we discuss examples of how the major venom protein superfamilies have diversified, as well as how newer technologies in the omics fields, such as genomics, transcriptomics, and proteomics, can be used to characterize both known and unknown toxins.Because toxins are bioactive molecules with a rich diversity of activities, they can be useful as therapeutic and diagnostic agents, and successful examples of toxin applications in these areas are also reviewed. With the current rapid pace of technology, snake venom research and its applications will only continue to expand.

Alternative mRNA Splicing in Three Venom Families Underlying a Possible Production of Divergent Venom Proteins of the Habu Snake, Protobothrops flavoviridis

Snake venoms are complex mixtures of toxic proteins encoded by various gene families that function synergistically to incapacitate prey. A huge repertoire of snake venom genes and proteins have been reported, and alternative splicing is suggested to be involved in the production of divergent gene transcripts. However, a genome-wide survey of the transcript repertoire and the extent of alternative splicing still remains to be determined. In this study, the comprehensive analysis of transcriptomes in the venom gland was achieved by using PacBio sequencing. Extensive alternative splicing was observed in three venom protein gene families, metalloproteinase (MP), serine protease (SP), and vascular endothelial growth factors (VEGF). Eleven MP and SP genes and a VEGF gene are expressed as a total of 81, 61, and 8 transcript variants, respectively. In the MP gene family, individual genes are transcribed into different classes of MPs by alternative splicing. We also observed trans-splicing among the clustered SP genes. No other venom genes as well as non-venom counterpart genes exhibited alternative splicing. Our results thus indicate a potential contribution of mRNA alternative and trans-splicing in the production of highly variable transcripts of venom genes in the habu snake.

Genome-wide mining of microsatellites in king cobra (Ophiophagus hannah) and cross-species development of tetranucleotide SSR markers in Chinese cobra (Naja atra)

The complete genome sequence provides the opportunity for genome-wide and coding region analysis of SSRs in the king cobra and for cross-species identification of microsatellite markers in the Chinese cobra. In the Ophiophagus hannah genome, tetranucleotide repeats (38.03%) were the most abundant category, followed by dinucleotides (23.03%), pentanucleotides (13.07%), mononucleotides (11.78%), trinucleotides (11.49%) and hexanucleotides (2.6%). Twenty predominant motifs in the O. hannah genome were (A)n (C)n, (AC)n, (AG)n, (AT)n, (AGG)n, (AAT)n, (AAG)n, (AAC)n, (ATG)n, (ATAG)n, (AAGG)n, (ATCT)n, (CCTT)n, (ATTT)n, (AAAT)n, (AATAG)n, (ATTCT)n, (ATATGT)n, (AGATAT)n. In total, 4344 SSRs were found in coding sequences (CDSs). Tetranucleotides (52.79%) were the most abundant microsatellite type in CDS, followed by trinucleotides (28.50%), dinucleotides (11.02%), pentanucleotides (4.42%), mononucleotides (1.77%), and hexanucleotides (1.50%). A total of 984 CDSs containing microsatellites were assigned 11152 Gene Ontology (GO) functional terms. Gene Ontology (GO) analysis demonstrated that cellular process, cell and binding were the most frequent GO terms in biological process, cellular component and molecular function, respectively. Thirty-two novel highly polymorphic (PIC > 0.5) SSR markers for Naja atra were developed from cross-species amplification based on the tetranucleotide microsatellite sequences in the king cobra genome. The number of alleles (N_A) per locus had between 3 and 11 alleles with an average of 6.5, the polymorphism information content (PIC) value ranged from 0.521 to 0.858 (average = 0.707), the observed heterozygosity (Ho) of 32 microsatellite loci ranged from 0.292 to 0.875 (mean = 0.678), the expected heterozygosity (H_E) ranged from 0.561 to 0.889 (average = 0.761), and 3 microsatellite loci exhibited statistically significant departure from Hardy-Weinberg equilibrium (HWE) after Bonferroni correction (p < 0.003).

Molecular Adaptations for Sensing and Securing Prey and Insight into Amniote Genome Diversity from the Garter Snake Genome

Colubridae represents the most phenotypically diverse and speciose family of snakes, yet no well-assembled and annotated genome exists for this lineage. Here, we report and analyze the genome of the garter snake, Thamnophis sirtalis, a colubrid snake that is an important model species for research in evolutionary biology, physiology, genomics, behavior, and the evolution of toxin resistance. Using the garter snake genome, we show how snakes have evolved numerous adaptations for sensing and securing prey, and identify features of snake genome structure that provide insight into the evolution of amniote genomes. Analyses of the garter snake and other squamate reptile genomes highlight shifts in repeat element abundance and expansion within snakes, uncover evidence of genes under positive selection, and provide revised neutral substitution rate estimates for squamates. Our identification of Z and W sex chromosome-specific scaffolds provides evidence for multiple origins of sex chromosome systems in snakes and demonstrates the value of this genome for studying sex chromosome evolution. Analysis of gene duplication and loss in visual and olfactory gene families supports a dim-light ancestral condition in snakes and indicates that olfactory receptor repertoires underwent an expansion early in snake evolution. Additionally, we provide some of the first links between secreted venom proteins, the genes that encode them, and their evolutionary origins in a rear-fanged colubrid snake, together with new genomic insight into the coevolutionary arms race between garter snakes and highly toxic newt prey that led to toxin resistance in garter snakes.

The habu genome reveals accelerated evolution of venom protein genes

Evolution of novel traits is a challenging subject in biological research. Several snake lineages developed elaborate venom systems to deliver complex protein mixtures for prey capture. To understand mechanisms involved in snake venom evolution, we decoded here the ~1.4-Gb genome of a habu, Protobothrops flavoviridis. We identified 60 snake venom protein genes (SV) and 224 non-venom paralogs (NV), belonging to 18 gene families. Molecular phylogeny reveals early divergence of SV and NV genes, suggesting that one of the four copies generated through two rounds of whole-genome duplication was modified for use as a toxin. Among them, both SV and NV genes in four major components were extensively duplicated after their diversification, but accelerated evolution is evident exclusively in the SV genes. Both venom-related SV and NV genes are significantly enriched in microchromosomes. The present study thus provides a genetic background for evolution of snake venom composition.

Squamate reptiles challenge paradigms of genomic repeat element evolution set by birds and mammals

Broad paradigms of vertebrate genomic repeat element evolution have been largely shaped by analyses of mammalian and avian genomes. Here, based on analyses of genomes sequenced from over 60 squamate reptiles (lizards and snakes), we show that patterns of genomic repeat landscape evolution in squamates challenge such paradigms. Despite low variance in genome size, squamate genomes exhibit surprisingly high variation among species in abundance (ca. 25–73% of the genome) and composition of identifiable repeat elements. We also demonstrate that snake genomes have experienced microsatellite seeding by transposable elements at a scale unparalleled among eukaryotes, leading to some snake genomes containing the highest microsatellite content of any known eukaryote. Our analyses of transposable element evolution across squamates also suggest that lineage-specific variation in mechanisms of transposable element activity and silencing, rather than variation in species-specific demography, may play a dominant role in driving variation in repeat element landscapes across squamate phylogeny.

Large-scale patterns of genomic repeat element evolution have been studied mainly in birds and mammals. Here, the authors analyze the genomes of over 60 squamate reptiles and show high variation in repeat elements compared to mammals and birds, and particularly high microsatellite seeding in snakes.

The annotation of repetitive elements in the genome of channel catfish (Ictalurus punctatus)

Channel catfish (Ictalurus punctatus) is a highly adaptive species and has been used as a research model for comparative immunology, physiology, and toxicology among ectothermic vertebrates. It is also economically important for aquaculture. As such, its reference genome was generated and annotated with protein coding genes. However, the repetitive elements in the catfish genome are less well understood. In this study, over 417.8 Megabase (MB) of repetitive elements were identified and characterized in the channel catfish genome. Among them, the DNA/TcMar-Tc1 transposons are the most abundant type, making up ~20% of the total repetitive elements, followed by the microsatellites (14%). The prevalence of repetitive elements, especially the mobile elements, may have provided a driving force for the evolution of the catfish genome. A number of catfish-specific repetitive elements were identified including the previously reported Xba elements whose divergence rate was relatively low, slower than that in untranslated regions of genes but faster than the protein coding sequences, suggesting its evolutionary restrictions.

The Evolution of LINE-1 in Vertebrates

The abundance and diversity of the LINE-1 (L1) retrotransposon differ greatly among vertebrates. Mammalian genomes contain hundreds of thousands L1s that have accumulated since the origin of mammals. A single group of very similar elements is active at a time in mammals, thus a single lineage of active families has evolved in this group. In contrast, non-mammalian genomes (fish, amphibians, reptiles) harbor a large diversity of concurrently transposing families, which are all represented by very small number of recently inserted copies. Why the pattern of diversity and abundance of L1 is so different among vertebrates remains unknown. To address this issue, we performed a detailed analysis of the evolution of active L1 in 14 mammals and in 3 non-mammalian vertebrate model species. We examined the evolution of base composition and codon bias, the general structure, and the evolution of the different domains of L1 (5′UTR, ORF1, ORF2, 3′UTR). L1s differ substantially in length, base composition, and structure among vertebrates. The most variation is found in the 5′UTR, which is longer in amniotes, and in the ORF1, which tend to evolve faster in mammals. The highly divergent L1 families of lizard, frog, and fish share species-specific features suggesting that they are subjected to the same functional constraints imposed by their host. The relative conservation of the 5′UTR and ORF1 in non-mammalian vertebrates suggests that the repression of transposition by the host does not act in a sequence-specific manner and did not result in an arms race, as is observed in mammals.

Evolution and Diversity of Transposable Elements in Vertebrate Genomes

Transposable elements (TEs) are selfish genetic elements that mobilize in genomes via transposition or retrotransposition and often make up large fractions of vertebrate genomes. Here, we review the current understanding of vertebrate TE diversity and evolution in the context of recent advances in genome sequencing and assembly techniques. TEs make up 4-60% of assembled vertebrate genomes, and deeply branching lineages such as ray-finned fishes and amphibians generally exhibit a higher TE diversity than the more recent radiations of birds and mammals. Furthermore, the list of taxa with exceptional TE landscapes is growing. We emphasize that the current bottleneck in genome analyses lies in the proper annotation of TEs and provide examples where superficial analyses led to misleading conclusions about genome evolution. Finally, recent advances in long-read sequencing will soon permit access to TE-rich genomic regions that previously resisted assembly including the gigantic, TE-rich genomes of salamanders and lungfishes.

Snake Genome Sequencing: Results and Future Prospects

Snake genome sequencing is in its infancy—very much behind the progress made in sequencing the genomes of humans, model organisms and pathogens relevant to biomedical research, and agricultural species. We provide here an overview of some of the snake genome projects in progress, and discuss the biological findings, with special emphasis on toxinology, from the small number of draft snake genomes already published. We discuss the future of snake genomics, pointing out that new sequencing technologies will help overcome the problem of repetitive sequences in assembling snake genomes. Genome sequences are also likely to be valuable in examining the clustering of toxin genes on the chromosomes, in designing recombinant antivenoms and in studying the epigenetic regulation of toxin gene expression.

Estimation of long terminal repeat element content in the Helicoverpa zea genome from high-throughput sequencing of bacterial artificial chromosome pools

The lepidopteran pest insect Helicoverpa zea feeds on cultivated corn and cotton across the Americas where control remains challenging owing to the evolution of resistance to chemical and transgenic insecticidal toxins, yet genomic resources remain scarce for this species. A bacterial artificial chromosome (BAC) library having a mean genomic insert size of 145 ± 20 kbp was created from a laboratory strain of H. zea, which provides ∼12.9-fold coverage of a 362.8 ± 8.8 Mbp (0.37 ± 0.09 pg) flow cytometry estimated haploid genome size. Assembly of Illumina HiSeq 2000 reads generated from 14 pools that encompassed all BAC clones resulted in 165 485 genomic contigs (N₅₀ = 3262 bp; 324.6 Mbp total). Long terminal repeat (LTR) protein coding regions annotated from 181 contigs included 30 Ty1/copia, 78 Ty3/gypsy, and 73 BEL/Pao elements, of which 60 (33.1%) encoded all five functional polyprotein (pol) domains. Approximately 14% of LTR elements are distributed non-randomly across pools of BAC clones.

[STRUCTURE AND DISTRIBUTION OF THE RETROTRANSPOSON BOV-B LINE]

The classification of mobile elements was discussed. Special attention was devoted to the retroelement of the LINE group: retrotransposon Bov-B LINE. The history of its origin and distribution in the nature was considered. The results of the phenomenon of horizontal transition of the retrotransposon Bov-B LINE between evolutionally distant classes were discussed.

Tyrosine Recombinase Retrotransposons and Transposons

Retrotransposons carrying tyrosine recombinases (YR) are widespread in eukaryotes. The first described tyrosine recombinase mobile element, DIRS1, is a retroelement from the slime mold Dictyostelium discoideum. The YR elements are bordered by terminal repeats related to their replication via free circular dsDNA intermediates. Site-specific recombination is believed to integrate the circle without creating duplications of the target sites. Recently a large number of YR retrotransposons have been described, including elements from fungi (mucorales and basidiomycetes), plants (green algae) and a wide range of animals including nematodes, insects, sea urchins, fish, amphibia and reptiles. YR retrotransposons can be divided into three major groups: the DIRS elements, PAT-like and the Ngaro elements. The three groups form distinct clades on phylogenetic trees based on alignments of reverse transcriptase/ribonuclease H (RT/RH) and YR sequences, and also having some structural distinctions. A group of eukaryote DNA transposons, cryptons, also carry tyrosine recombinases. These DNA transposons do not encode a reverse transcriptase. They have been detected in several pathogenic fungi and oomycetes. Sequence comparisons suggest that the crypton YRs are related to those of the YR retrotransposons. We suggest that the YR retrotransposons arose from the combination of a crypton-like YR DNA transposon and the RT/RH encoding sequence of a retrotransposon. This acquisition must have occurred at a very early point in the evolution of eukaryotes.

Microsatellite landscape evolutionary dynamics across 450 million years of vertebrate genome evolution

The evolutionary dynamics of simple sequence repeats (SSRs or microsatellites) across the vertebrate tree of life remain largely undocumented and poorly understood. In this study, we analyzed patterns of genomic microsatellite abundance and evolution across 71 vertebrate genomes. The highest abundances of microsatellites exist in the genomes of ray-finned fishes, squamate reptiles, and mammals, while crocodilian, turtle, and avian genomes exhibit reduced microsatellite landscapes. We used comparative methods to infer evolutionary rates of change in microsatellite abundance across vertebrates and to highlight particular lineages that have experienced unusually high or low rates of change in genomic microsatellite abundance. Overall, most variation in microsatellite content, abundance, and evolutionary rate is observed among major lineages of reptiles, yet we found that several deeply divergent clades (i.e., squamate reptiles and mammals) contained relatively similar genomic microsatellite compositions. Archosauromorph reptiles (turtles, crocodilians, and birds) exhibit reduced genomic microsatellite content and the slowest rates of microsatellite evolution, in contrast to squamate reptile genomes that have among the highest rates of microsatellite evolution. Substantial branch-specific shifts in SSR content in primates, monotremes, rodents, snakes, and fish are also evident. Collectively, our results support multiple major shifts in microsatellite genomic landscapes among vertebrates.

The contribution of transposable elements to size variations between four teleost genomes

Background

Teleosts are unique among vertebrates, with a wide range of haploid genome sizes in very close lineages, varying from less than 400 mega base pairs (Mb) for pufferfish to over 3000 Mb for salmon. The cause of the difference in genome size remains largely unexplained.

Results

In this study, we reveal that the differential success of transposable elements (TEs) correlates with the variation of genome size across four representative teleost species (zebrafish, medaka, stickleback, and tetraodon). The larger genomes represent a higher diversity within each clade (superfamily) and family and a greater abundance of TEs compared with the smaller genomes; zebrafish, representing the largest genome, shows the highest diversity and abundance of TEs in its genome, followed by medaka and stickleback; while the tetraodon, representing the most compact genome, displays the lowest diversity and density of TEs in its genome. Both of Class I (retrotransposons) and Class II TEs (DNA transposons) contribute to the difference of TE accumulation of teleost genomes, however, Class II TEs are the major component of the larger teleost genomes analyzed and the most important contributors to genome size variation across teleost lineages. The hAT and Tc1/Mariner superfamilies are the major DNA transposons of all four investigated teleosts. Divergence distribution revealed contrasting proliferation dynamics both between clades of retrotransposons and between species. The TEs within the larger genomes of the zebrafish and medaka represent relatively stronger activity with an extended time period during the evolution history, in contrast with the very young activity in the smaller stickleback genome, or the very low level of activity in the tetraodon genome.

Conclusion

Overall, our data shows that teleosts represent contrasting profiles of mobilomes with a differential density, diversity and activity of TEs. The differences in TE accumulation, dominated by DNA transposons, explain the main size variations of genomes across the investigated teleost species, and the species differences in both diversity and activity of TEs contributed to the variations of TE accumulations across the four teleost species. TEs play major roles in teleost genome evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s13100-016-0059-7) contains supplementary material, which is available to authorized users.

No Accumulation of Transposable Elements in Asexual Arthropods

Transposable elements (TEs) and other repetitive DNA can accumulate in the absence of recombination, a process contributing to the degeneration of Y-chromosomes and other nonrecombining genome portions. A similar accumulation of repetitive DNA is expected for asexually reproducing species, given their entire genome is effectively nonrecombining. We tested this expectation by comparing the whole-genome TE loads of five asexual arthropod lineages and their sexual relatives, including asexual and sexual lineages of crustaceans (Daphnia water fleas), insects (Leptopilina wasps), and mites (Oribatida). Surprisingly, there was no evidence for increased TE load in genomes of asexual as compared to sexual lineages, neither for all classes of repetitive elements combined nor for specific TE families. Our study therefore suggests that nonrecombining genomes do not accumulate TEs like nonrecombining genomic regions of sexual lineages. Even if a slight but undetected increase of TEs were caused by asexual reproduction, it appears to be negligible compared to variance between species caused by processes unrelated to reproductive mode. It remains to be determined if molecular mechanisms underlying genome regulation in asexuals hamper TE activity. Alternatively, the differences in TE dynamics between nonrecombining genomes in asexual lineages versus nonrecombining genome portions in sexual species might stem from selection for benign TEs in asexual lineages because of the lack of genetic conflict between TEs and their hosts and/or because asexual lineages may only arise from sexual ancestors with particularly low TE loads.

Endogenous hepadnaviruses, bornaviruses and circoviruses in snakes

We report the discovery of endogenous viral elements (EVEs) from Hepadnaviridae, Bornaviridae and Circoviridae in the speckled rattlesnake, Crotalus mitchellii, the first viperid snake for which a draft whole genome sequence assembly is available. Analysis of the draft assembly reveals genome fragments from the three virus families were inserted into the genome of this snake over the past 50 Myr. Cross-species PCR screening of orthologous loci and computational scanning of the python and king cobra genomes reveals that circoviruses integrated most recently (within the last approx. 10 Myr), whereas bornaviruses and hepadnaviruses integrated at least approximately 13 and approximately 50 Ma, respectively. This is, to our knowledge, the first report of circo-, borna- and hepadnaviruses in snakes and the first characterization of non-retroviral EVEs in non-avian reptiles. Our study provides a window into the historical dynamics of viruses in these host lineages and shows that their evolution involved multiple host-switches between mammals and reptiles.

A genome survey sequencing of the Java mouse deer (Tragulus javanicus) adds new aspects to the evolution of lineage specific retrotransposons in Ruminantia (Cetartiodactyla)

Ruminantia, the ruminating, hoofed mammals (cow, deer, giraffe and allies) are an unranked artiodactylan clade. Around 50-60 million years ago the BovB retrotransposon entered the ancestral ruminantian genome through horizontal gene transfer. A survey genome screen using 454-pyrosequencing of the Java mouse deer (Tragulus javanicus) and the lesser kudu (Tragelaphus imberbis) was done to investigate and to compare the landscape of transposable elements within Ruminantia. The family Tragulidae (mouse deer) is the only representative of Tragulina and phylogenetically important, because it represents the earliest divergence in Ruminantia. The data analyses show that, relative to other ruminantian species, the lesser kudu genome has seen an expansion of BovB Long INterspersed Elements (LINEs) and BovB related Short INterspersed Elements (SINEs) like BOVA2. In comparison the genome of Java mouse deer has fewer BovB elements than other ruminants, especially Bovinae, and has in addition a novel CHR-3 SINE most likely propagated by LINE-1. By contrast the other ruminants have low amounts of CHR SINEs but high numbers of actively propagating BovB-derived and BovB-propagated SINEs. The survey sequencing data suggest that the transposable element landscape in mouse deer (Tragulina) is unique among Ruminantia, suggesting a lineage specific evolutionary trajectory that does not involve BovB mediated retrotransposition. This shows that the genomic landscape of mobile genetic elements can rapidly change in any lineage.

Comparative analysis of transposable elements highlights mobilome diversity and evolution in vertebrates

Transposable elements (TEs) are major components of vertebrate genomes, with major roles in genome architecture and evolution. In order to characterize both common patterns and lineage-specific differences in TE content and TE evolution, we have compared the mobilomes of 23 vertebrate genomes, including 10 actinopterygian fish, 11 sarcopterygians, and 2 nonbony vertebrates. We found important variations in TE content (from 6% in the pufferfish tetraodon to 55% in zebrafish), with a more important relative contribution of TEs to genome size in fish than in mammals. Some TE superfamilies were found to be widespread in vertebrates, but most elements showed a more patchy distribution, indicative of multiple events of loss or gain. Interestingly, loss of major TE families was observed during the evolution of the sarcopterygian lineage, with a particularly strong reduction in TE diversity in birds and mammals. Phylogenetic trends in TE composition and activity were detected: Teleost fish genomes are dominated by DNA transposons and contain few ancient TE copies, while mammalian genomes have been predominantly shaped by nonlong terminal repeat retrotransposons, along with the persistence of older sequences. Differences were also found within lineages: The medaka fish genome underwent more recent TE amplification than the related platyfish, as observed for LINE retrotransposons in the mouse compared with the human genome. This study allows the identification of putative cases of horizontal transfer of TEs, and to tentatively infer the composition of the ancestral vertebrate mobilome. Taken together, the results obtained highlight the importance of TEs in the structure and evolution of vertebrate genomes, and demonstrate their major impact on genome diversity both between and within lineages.

The finding of a group IIE phospholipase A2 gene in a specified segment of Protobothrops flavoviridis genome and its possible evolutionary relationship to group IIA phospholipase A2 genes

The genes encoding group IIE phospholipase A₂, abbreviated as IIE PLA₂, and its 5' and 3' flanking regions of Crotalinae snakes such as Protobothrops flavoviridis, P. tokarensis, P. elegans, and Ovophis okinavensis, were found and sequenced. The genes consisted of four exons and three introns and coded for 22 or 24 amino acid residues of the signal peptides and 134 amino acid residues of the mature proteins. These IIE PLA₂s show high similarity to those from mammals and Colubridae snakes. The high expression level of IIE PLA₂s in Crotalinae venom glands suggests that they should work as venomous proteins. The blast analysis indicated that the gene encoding OTUD3, which is ovarian tumor domain-containing protein 3, is located in the 3' downstream of IIE PLA₂ gene. Moreover, a group IIA PLA₂ gene was found in the 5' upstream of IIE PLA₂ gene linked to the OTUD3 gene (OTUD3) in the P. flavoviridis genome. It became evident that the specified arrangement of IIA PLA₂ gene, IIE PLA₂ gene, and OTUD3 in this order is common in the genomes of humans to snakes. The present finding that the genes encoding various secretory PLA₂s form a cluster in the genomes of humans to birds is closely related to the previous finding that six venom PLA₂ isozyme genes are densely clustered in the so-called NIS-1 fragment of the P. flavoviridis genome. It is also suggested that venom IIA PLA₂ genes may be evolutionarily derived from the IIE PLA₂ gene.

Two low coverage bird genomes and a comparison of reference-guided versus de novo genome assemblies

As a greater number and diversity of high-quality vertebrate reference genomes become available, it is increasingly feasible to use these references to guide new draft assemblies for related species. Reference-guided assembly approaches may substantially increase the contiguity and completeness of a new genome using only low levels of genome coverage that might otherwise be insufficient for de novo genome assembly. We used low-coverage (∼3.5-5.5x) Illumina paired-end sequencing to assemble draft genomes of two bird species (the Gunnison Sage-Grouse, Centrocercus minimus, and the Clark's Nutcracker, Nucifraga columbiana). We used these data to estimate de novo genome assemblies and reference-guided assemblies, and compared the information content and completeness of these assemblies by comparing CEGMA gene set representation, repeat element content, simple sequence repeat content, and GC isochore structure among assemblies. Our results demonstrate that even lower-coverage genome sequencing projects are capable of producing informative and useful genomic resources, particularly through the use of reference-guided assemblies.