Mobilization of retrotransposons as a cause of chromosomal diversification and rapid speciation: the case for the Antarctic teleost genus Trematomus

The complete mitochondrial genome of Trematomus hansoni Boulenger, 1902 (Perciformes, Nototheniidae)

The striped notothen Trematomus hansoni is an Antarctic fish species belonging to the family Nototheniidae (cod icefishes) that is distributed throughout the Southern Ocean. In this study, the complete mitochondrial genome of T. hansoni was sequenced using an Illumina MiSeq platform. The circular mitochondrial genome is 19,218 bp long and contains 13 protein-coding genes, 23 tRNA genes, two rRNA genes, and one control region. Notably, there are two trnG-UCC genes and the second gene, located between trnE-UUC and trnI-GAU, has no D-arm structure. The base composition is 56.18% of A + T and 43.82% of G + C. The phylogenetic analysis supports that T. hansoni is grouped into a single clade with T. bernacchii. This study will be a valuable resource for further research on the phylogeny and evolution of the genus Trematomus.

A macroevolutionary role for chromosomal fusion and fission in Erebia butterflies

The impact of large-scale chromosomal rearrangements, such as fusions and fissions, on speciation is a long-standing conundrum. We assessed whether bursts of change in chromosome numbers resulting from chromosomal fusion or fission are related to increased speciation rates in Erebia, one of the most species-rich and karyotypically variable butterfly groups. We established a genome-based phylogeny and used state-dependent birth-death models to infer trajectories of karyotype evolution. We demonstrated that rates of anagenetic chromosomal changes (i.e., along phylogenetic branches) exceed cladogenetic changes (i.e., at speciation events), but, when cladogenetic changes occur, they are mostly associated with chromosomal fissions rather than fusions. We found that the relative importance of fusion and fission differs among Erebia clades of different ages and that especially in younger, more karyotypically diverse clades, speciation is more frequently associated with cladogenetic chromosomal changes. Overall, our results imply that chromosomal fusions and fissions have contrasting macroevolutionary roles and that large-scale chromosomal rearrangements are associated with bursts of species diversification.

Hiker, a new family of DNA transposons encoding transposases with DD35E motifs, displays a distinct phylogenetic relationship with most known DNA transposon families of IS630-Tc1-mariner (ITm)

DNA transposons play a crucial role in determining the size and structure of eukaryotic genomes. In this study, a new family of IS630-Tc1-mariner (ITm) DNA transposons, named Hiker (HK), was identified. HK is characterized by a DD35E catalytic domain and is distinct from all previously known families of the ITm group. Phylogenetic analyses showed that DD35E/Hiker forms a monophyletic clade with DD34E/Gambol, indicating that they may represent a separate superfamily of ITm. A total of 178 Hiker species were identified, with 170 found mainly in Actinopterygii, one in Chondrichthyes, six in Anura and one in Mollusca. Gambol (GM), on the other hand, are found in invertebrates, with 18 in Arthropoda and one in Platyhelminthes. Hiker transposons have a total length ranging from 2.14 to 3.67 kb and contain a single open reading frame that encodes a protein of approximately 370 amino acids (range 311-413 aa). They are flanked by short terminal inverted repeats (TIRs) of 16-30 base pairs and two base pair (TA) target-site duplications. In contrast, most transposons of the Gambol family have a total length of 1.35-5.96 kb, encode a transposase protein of approximately 350 amino acids (range 306-374 aa), and are flanked by TIRs that range from 32 to 1097 bp in length. Both Hiker and Gambol transposases have several conserved motifs, including helix-turn-helix (HTH) motifs and a DDE domain. Our study observed multiple amplification waves and repeated horizontal transfer (HT) events of HK transposons in vertebrate genomes, indicating their role in diversifying and shaping the genomes of Actinopterygii, Chondrichthyes, and Anura. Conversely, GM transposons showed few Horizontal transfer events. According to cell-based transposition assays, most HK transposons are likely inactive due to the truncated DNA binding domains of their transposases. We present an updated classification of the ITm group based on these findings, which will enhance the understanding of both the evolution of ITm transposons and that of their hosts.

The genome of the glasshouse plant noble rhubarb (Rheum nobile) provides a window into alpine adaptation

Glasshouse plants are species that trap warmth via specialized morphology and physiology, mimicking a human glasshouse. In the Himalayan alpine region, the highly specialized glasshouse morphology has independently evolved in distinct lineages to adapt to intensive UV radiation and low temperature. Here we demonstrate that the glasshouse structure - specialized cauline leaves - is highly effective in absorbing UV light but transmitting visible and infrared light, creating an optimal microclimate for the development of reproductive organs. We reveal that this glasshouse syndrome has evolved at least three times independently in the rhubarb genus Rheum. We report the genome sequence of the flagship glasshouse plant Rheum nobile and identify key genetic network modules in association with the morphological transition to specialized glasshouse leaves, including active secondary cell wall biogenesis, upregulated cuticular cutin biosynthesis, and suppression of photosynthesis and terpenoid biosynthesis. The distinct cell wall organization and cuticle development might be important for the specialized optical property of glasshouse leaves. We also find that the expansion of LTRs has likely played an important role in noble rhubarb adaptation to high elevation environments. Our study will enable additional comparative analyses to identify the genetic basis underlying the convergent occurrence of glasshouse syndrome.

Insights into the Genomics of Clownfish Adaptive Radiation: The Genomic Substrate of the Diversification

Clownfishes are an iconic group of coral reef fishes that evolved a mutualistic interaction with sea anemones, which triggered the rapid diversification of the group. Following the emergence of this mutualism, clownfishes diversified into different ecological niches and developed convergent phenotypes associated with their host use. The genetic basis of the initial acquisition of the mutualism with host anemones has been described, but the genomic architecture underlying clownfish diversification once the mutualism was established and the extent to which clownfish phenotypic convergence originated through shared genetic mechanisms are still unknown. Here, we investigated these questions by performing comparative genomic analyses on the available genomic data of five pairs of closely related but ecologically divergent clownfish species. We found that clownfish diversification was characterized by bursts of transposable elements, an overall accelerated coding evolution, incomplete lineage sorting, and ancestral hybridization events. Additionally, we detected a signature of positive selection in 5.4% of the clownfish genes. Among them, five presented functions associated with social behavior and ecology, and they represent candidate genes involved in the evolution of the size-based hierarchical social structure so particular to clownfishes. Finally, we found genes with patterns of either relaxation or intensification of purifying selection and signals of positive selection linked with clownfish ecological divergence, suggesting some level of parallel evolution during the diversification of the group. Altogether, this work provides the first insights into the genomic substrate of clownfish adaptive radiation and integrates the growing collection of studies investigating the genomic mechanisms governing species diversification.

Chromosome-Level Genome Assembly and Circadian Gene Repertoire of the Patagonia Blennie Eleginops maclovinus-The Closest Ancestral Proxy of Antarctic Cryonotothenioids

The basal South American notothenioid Eleginops maclovinus (Patagonia blennie or róbalo) occupies a uniquely important phylogenetic position in Notothenioidei as the singular closest sister species to the Antarctic cryonotothenioid fishes. Its genome and the traits encoded therein would be the nearest representatives of the temperate ancestor from which the Antarctic clade arose, providing an ancestral reference for deducing polar derived changes. In this study, we generated a gene- and chromosome-complete assembly of the E. maclovinus genome using long read sequencing and HiC scaffolding. We compared its genome architecture with the more basally divergent Cottoperca gobio and the derived genomes of nine cryonotothenioids representing all five Antarctic families. We also reconstructed a notothenioid phylogeny using 2918 proteins of single-copy orthologous genes from these genomes that reaffirmed E. maclovinus' phylogenetic position. We additionally curated E. maclovinus' repertoire of circadian rhythm genes, ascertained their functionality by transcriptome sequencing, and compared its pattern of gene retention with C. gobio and the derived cryonotothenioids. Through reconstructing circadian gene trees, we also assessed the potential role of the retained genes in cryonotothenioids by referencing to the functions of the human orthologs. Our results found E. maclovinus to share greater conservation with the Antarctic clade, solidifying its evolutionary status as the direct sister and best suited ancestral proxy of cryonotothenioids. The high-quality genome of E. maclovinus will facilitate inquiries into cold derived traits in temperate to polar evolution, and conversely on the paths of readaptation to non-freezing habitats in various secondarily temperate cryonotothenioids through comparative genomic analyses.

Genomics of Secondarily Temperate Adaptation in the Only Non-Antarctic Icefish

White-blooded Antarctic icefishes, a family within the adaptive radiation of Antarctic notothenioid fishes, are an example of extreme biological specialization to both the chronic cold of the Southern Ocean and life without hemoglobin. As a result, icefishes display derived physiology that limits them to the cold and highly oxygenated Antarctic waters. Against these constraints, remarkably one species, the pike icefish Champsocephalus esox, successfully colonized temperate South American waters. To study the genetic mechanisms underlying secondarily temperate adaptation in icefishes, we generated chromosome-level genome assemblies of both C. esox and its Antarctic sister species, Champsocephalus gunnari. The C. esox genome is similar in structure and organization to that of its Antarctic congener; however, we observe evidence of chromosomal rearrangements coinciding with regions of elevated genetic divergence in pike icefish populations. We also find several key biological pathways under selection, including genes related to mitochondria and vision, highlighting candidates behind temperate adaptation in C. esox. Substantial antifreeze glycoprotein (AFGP) pseudogenization has occurred in the pike icefish, likely due to relaxed selection following ancestral escape from Antarctica. The canonical AFGP locus organization is conserved in C. esox and C. gunnari, but both show a translocation of two AFGP copies to a separate locus, previously unobserved in cryonotothenioids. Altogether, the study of this secondarily temperate species provides an insight into the mechanisms underlying adaptation to ecologically disparate environments in this otherwise highly specialized group.

Pathways to polar adaptation in fishes revealed by long-read sequencing

Long-read sequencing is driving a new reality for genome science in which highly contiguous assemblies can be produced efficiently with modest resources. Genome assemblies from long-read sequences are particularly exciting for understanding the evolution of complex genomic regions that are often difficult to assemble. In this study, we utilized long-read sequencing data to generate a high-quality genome assembly for an Antarctic eelpout, Ophthalmolycus amberensis, the first for the globally distributed family Zoarcidae. We used this assembly to understand how O. amberensis has adapted to the harsh Southern Ocean and compared it to another group of Antarctic fishes: the notothenioids. We showed that selection has largely acted on different targets in eelpouts relative to notothenioids. However, we did find some overlap; in both groups, genes involved in membrane structure, thermal tolerance and vision have evidence of positive selection. We found evidence for historical shifts of transposable element activity in O. amberensis and other polar fishes, perhaps reflecting a response to environmental change. We were specifically interested in the evolution of two complex genomic loci known to underlie key adaptations to polar seas: haemoglobin and antifreeze proteins (AFPs). We observed unique evolution of the haemoglobin MN cluster in eelpouts and related fishes in the suborder Zoarcoidei relative to other Perciformes. For AFPs, we identified the first species in the suborder with no evidence of afpIII sequences (Cebidichthys violaceus) in the genomic region where they are found in all other Zoarcoidei, potentially reflecting a lineage-specific loss of this cluster. Beyond polar fishes, our results highlight the power of long-read sequencing to understand genome evolution.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

Tail Wags Dog’s SINE: Retropositional Mechanisms of Can SINE Depend on Its A-Tail Structure

Simple Summary

The genomes of higher organisms including humans are invaded by millions of repetitive elements (transposons), which can sometimes be deleterious or beneficial for hosts. Many aspects of the mechanisms underlying the expansion of transposons in the genomes remain unclear. Short retrotransposons (SINEs) are one of the most abundant classes of genomic repeats. Their amplification relies on two major processes: transcription and reverse transcription. Here, short retrotransposons of dogs and other canids called Can SINE were analyzed. Their amplification was extraordinarily active in the wolf and, particularly, dog breeds relative to other canids. We also studied a variation of their transcription mechanism involving the polyadenylation of transcripts. An analysis of specific signals involved in this process allowed us to conclude that Can SINEs could alternate amplification with and without polyadenylation in their evolution. Understanding the mechanisms of transposon replication can shed light on the mechanisms of genome function.

Abstract

SINEs, non-autonomous short retrotransposons, are widespread in mammalian genomes. Their transcripts are generated by RNA polymerase III (pol III). Transcripts of certain SINEs can be polyadenylated, which requires polyadenylation and pol III termination signals in their sequences. Our sequence analysis divided Can SINEs in canids into four subfamilies, older a1 and a2 and younger b1 and b2. Can_b2 and to a lesser extent Can_b1 remained retrotranspositionally active, while the amplification of Can_a1 and Can_a2 ceased long ago. An extraordinarily high Can amplification was revealed in different dog breeds. Functional polyadenylation signals were analyzed in Can subfamilies, particularly in fractions of recently amplified, i.e., active copies. The transcription of various Can constructs transfected into HeLa cells proposed AATAAA and (TC)_n as functional polyadenylation signals. Our analysis indicates that older Can subfamilies (a1, a2, and b1) with an active transcription terminator were amplified by the T⁺ mechanism (with polyadenylation of pol III transcripts). In the currently active Can_b2 subfamily, the amplification mechanisms with (T⁺) and without the polyadenylation of pol III transcripts (T⁻) irregularly alternate. The active transcription terminator tends to shorten, which renders it nonfunctional and favors a switch to the T⁻ retrotransposition. The activity of a truncated terminator is occasionally restored by its elongation, which rehabilitates the T⁺ retrotransposition for a particular SINE copy.

Comparative Analysis of the pIgR Gene from the Antarctic Teleost Trematomus bernacchii Reveals Distinctive Features of Cold-Adapted Notothenioidei

The IgM and IgT classes were previously identified and characterized in the Antarctic teleost Trematomus bernacchii, a species belonging to the Perciform suborder Notothenoidei. Herein, we characterized the gene encoding the polymeric immunoglobulin receptor (pIgR) in the same species and compared it to the pIgR of multiple teleost species belonging to five perciform suborders, including 11 Antarctic and 1 non-Antarctic (Cottoperca gobio) notothenioid species, the latter living in the less-cold peri-Antarctic sea. Antarctic pIgR genes displayed particularly long introns marked by sites of transposable elements and transcription factors. Furthermore, analysis of T. bernacchii pIgR cDNA unveiled multiple amino acid substitutions unique to the Antarctic species, all introducing adaptive features, including N-glycosylation sequons. Interestingly, C. gobio shared most features with the other perciforms rather than with the cold-adapted relatives. T. bernacchii pIgR transcripts were predominantly expressed in mucosal tissues, as indicated by q-PCR and in situ hybridization analysis. These results suggest that in cold-adapted species, pIgR preserved its fundamental role in mucosal immune defense, although remarkable gene structure modifications occurred.

Chromosome-level genome assembly of the aquatic plant Nymphoides indica reveals transposable element bursts and NBS-LRR gene family expansion shedding light on its invasiveness

Nymphoides indica, an aquatic plant, is an invasive species that causes both ecological and economic damage in North America and elsewhere. However, the lack of genomic data of N. indica limits the in-depth analysis of this invasive species. Here, we report a chromosome-level genome assembly of nine pseudochromosomes of N. indica with a total size of ∼ 520 Mb. More than half of the N. indica genome consists of transposable elements (TEs), and a higher density of TEs around genes may play a significant role in response to an ever-changing environment by regulating the nearby gene. Additionally, our analysis revealed that N. indica only experienced a gamma (γ) whole-genome triplication event. Functional enrichment of the N. indica-specific and expanded gene families highlighted genes involved in the responses to hypoxia and plant–pathogen interactions, which may strengthen the ability to adapt to external challenges and improve ecological fitness. Furthermore, we identified 160 members of the nucleotide-binding site and leucine-rich repeat gene family, which may be linked to the defence response. Collectively, the high-quality N. indica genome reported here opens a novel avenue to understand the evolution and rapid invasion of Nymphoides spp.

Rapid karyotypic evolution with high diploid number variation in a rare genus of bromeligenous frogs

Bromeligenous Crossodactylodes is a leptodactylid genus closely related to Paratelmatobius and Scythrophrys. The diploid number in all karyotyped species of these two latter genera is 24, which diverges from the modal diploid number (2n = 22) in the family. Here, we analyzed three species of Crossodactylodes and found karyotypes with 2n = 30, 2n = 32, and 2n = 36, diploid numbers that have not been reported in any other diploid leptodactylid species to date. Reconstruction of the ancestral chromosome number indicated that the diploid number changed from 22 to 24 in the common ancestor of Crossodactylodes, Paratelmatobius, and Scythrophrys, and that progressive increases in diploid number have occurred in Crossodactylodes. The large number of telocentric/subtelocentric chromosomes in karyotypes with higher diploid numbers raises the possibility that centric fissions may have occurred during the evolution of Paratelmatobiinae. Three metacentric chromosomes, probably involved in fission events, were inferred to be present in the common ancestor of all species of Crossodactylodes, but in C. bokermanni. Chromosome mapping of the satellite DNA PcP190 suggests homology between one arm of metacentric chromosome 1 of Crossodactylodes sp. 3 and telocentric chromosome 2 of C. itambe, supporting one of the presumed centric fission events.

OUP accepted manuscript

Animals rely on their sensory systems to inform them of ecologically relevant environmental variation. In the Southern Ocean, the thermal environment has remained between −1.9 and 5 °C for 15 Myr, yet we have no knowledge of how an Antarctic marine organism might sense their thermal habitat as we have yet to discover a thermosensitive ion channel that gates (opens/closes) below 10 °C. Here, we investigate the evolutionary dynamics of transient receptor potential (TRP) channels, which are the primary thermosensors in animals, within cryonotothenioid fishes—the dominant fish fauna of the Southern Ocean. We found cryonotothenioids have a similar complement of TRP channels as other teleosts (∼28 genes). Previous work has shown that thermosensitive gating in a given channel is species specific, and multiple channels act together to sense the thermal environment. Therefore, we combined evidence of changes in selective pressure, gene gain/loss dynamics, and the first sensory ganglion transcriptome in this clade to identify the best candidate TRP channels that might have a functional dynamic range relevant for frigid Antarctic temperatures. We concluded that TRPV1a, TRPA1b, and TRPM4 are the likeliest putative thermosensors, and found evidence of diversifying selection at sites across these proteins. We also put forward hypotheses for molecular mechanisms of other cryonotothenioid adaptations, such as reduced skeletal calcium deposition, sensing oxidative stress, and unusual magnesium homeostasis. By completing a comprehensive and unbiased survey of these genes, we lay the groundwork for functional characterization and answering long-standing thermodynamic questions of thermosensitive gating and protein adaptation to low temperatures.

Comparative analysis of transposable elements provides insights into genome evolution in the genus Camelus

BACKGROUND

Transposable elements (TEs) are common features in eukaryotic genomes that are known to affect genome evolution critically and to play roles in gene regulation. Vertebrate genomes are dominated by TEs, which can reach copy numbers in the hundreds of thousands. To date, details regarding the presence and characteristics of TEs in camelid genomes have not been made available.

RESULTS

We conducted a genome-wide comparative analysis of camelid TEs, focusing on the identification of TEs and elucidation of transposition histories in four species: Camelus dromedarius, C. bactrianus, C. ferus, and Vicugna pacos. Our TE library was created using both de novo structure-based and homology-based searching strategies ( https://github.com/kacst-bioinfo-lab/TE_ideintification_pipeline ). Annotation results indicated a similar proportion of each genomes comprising TEs (35-36%). Class I LTR retrotransposons comprised 16-20% of genomes, and mostly consisted of the endogenous retroviruses (ERVs) groups ERVL, ERVL-MaLR, ERV_classI, and ERV_classII. Non-LTR elements comprised about 12% of genomes and consisted of SINEs (MIRs) and the LINE superfamilies LINE1, LINE2, L3/CR1, and RTE clades. Least represented were the Class II DNA transposons (2%), consisting of hAT-Charlie, TcMar-Tigger, and Helitron elements and comprising about 1-2% of each genome.

CONCLUSIONS

The findings of the present study revealed that the distribution of transposable elements across camelid genomes is approximately similar. This investigation presents a characterization of TE content in four camelid to contribute to developing a better understanding of camelid genome architecture and evolution.

Adaptations and Diversity of Antarctic Fishes: A Genomic Perspective

Antarctic notothenioid fishes are the classic example of vertebrate adaptive radiation in a marine environment. Notothenioids diversified from a single common ancestor ∼25 Mya to more than 140 species today, and they represent ∼90% of fish biomass on the continental shelf of Antarctica. As they diversified in the cold Southern Ocean, notothenioids evolved numerous traits, including osteopenia, anemia, cardiomegaly, dyslipidemia, and aglomerular kidneys, that are beneficial or tolerated in their environment but are pathological in humans. Thus, notothenioids are models for understanding adaptive radiations, physiological and biochemical adaptations to extreme environments, and genetic mechanisms of human disease. Since 2014, 16 notothenioid genomes have been published, which enable a first-pass holistic analysis of the notothenioid radiation and the genetic underpinnings of novel notothenioid traits. Here, we review the notothenioid radiation from a genomic perspective and integrate our insights with recent observations from other fish radiations. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Prokaryotic and Eukaryotic Horizontal Transfer of Sailor (DD82E), a New Superfamily of IS630-Tc1-Mariner DNA Transposons

Simple Summary

Transposable elements, including DNA transposons, play a significant role in genetic material exchanges between prokaryotes and eukaryotes. Comparative profiling of the evolution pattern of DNA transposons between prokaryotes and eukaryotes may identify potential genetic material exchanges between them and provide insights into the evolutionary history of prokaryotic and eukaryotic genomes. The members of the IS630-Tc1-mariner (ITm) group may represent the most diverse and widely distributed DNA transposons in nature, and the discovery of new members of this group is highly expected based on the increasing availability of genome sequencing data. We discovered a new superfamily (termed Sailor) belonging to the ITm hyperfamily, which differed from the known superfamilies of Tc1/mariner, DDxD/pogo and DD34E/Gambol, regarding phylogenetic position and catalytic domain. Our data revealed that Sailor was distributed in both prokaryotes and eukaryotes and suggested that horizontal transfer (HT) events of Sailor may occur from prokaryotic to eukaryotic genomes. Finally, internal transmissions of Sailor in prokaryotes and eukaryotes were also detected.

Abstract

Here, a new superfamily of IS630-Tc1-mariner (ITm) DNA transposons, termed Sailor, is identified, that is characterized by a DD82E catalytic domain and is distinct from all previously known superfamilies of the ITm group. Phylogenetic analyses revealed that Sailor forms a monophyletic clade with a more intimate link to the clades of Tc1/mariner and DD34E/Gambol. Sailor was detected in both prokaryotes and eukaryotes and invaded a total of 256 species across six kingdoms. Sailor is present in nine species of bacteria, two species of plantae, four species of protozoa, 23 species of Chromista, 12 species of Fungi and 206 species of animals. Moreover, Sailor is extensively distributed in invertebrates (a total of 206 species from six phyla) but is absent in vertebrates. Sailor transposons are 1.38–6.98 kb in total length and encoded transposases of ~676 aa flanked by TIRs with lengths between 18, 1362 and 4 bp (TATA) target-site duplications. Furthermore, our analysis provided strong evidence of Sailor transmissions from prokaryotes to eukaryotes and internal transmissions in both. These data update the classification of the ITm group and will contribute to the understanding of the evolution of ITm transposons and that of their hosts.

Facilitating population genomics of non-model organisms through optimized experimental design for reduced representation sequencing

BACKGROUND

Genome-wide data are invaluable to characterize differentiation and adaptation of natural populations. Reduced representation sequencing (RRS) subsamples a genome repeatedly across many individuals. However, RRS requires careful optimization and fine-tuning to deliver high marker density while being cost-efficient. The number of genomic fragments created through restriction enzyme digestion and the sequencing library setup must match to achieve sufficient sequencing coverage per locus. Here, we present a workflow based on published information and computational and experimental procedures to investigate and streamline the applicability of RRS.

RESULTS

In an iterative process genome size estimates, restriction enzymes and size selection windows were tested and scaled in six classes of Antarctic animals (Ostracoda, Malacostraca, Bivalvia, Asteroidea, Actinopterygii, Aves). Achieving high marker density would be expensive in amphipods, the malacostracan target taxon, due to the large genome size. We propose alternative approaches such as mitogenome or target capture sequencing for this group. Pilot libraries were sequenced for all other target taxa. Ostracods, bivalves, sea stars, and fish showed overall good coverage and marker numbers for downstream population genomic analyses. In contrast, the bird test library produced low coverage and few polymorphic loci, likely due to degraded DNA.

CONCLUSIONS

Prior testing and optimization are important to identify which groups are amenable for RRS and where alternative methods may currently offer better cost-benefit ratios. The steps outlined here are easy to follow for other non-model taxa with little genomic resources, thus stimulating efficient resource use for the many pressing research questions in molecular ecology.

Extreme Y chromosome polymorphism corresponds to five male reproductive morphs of a freshwater fish

Loss of recombination between sex chromosomes often depletes Y chromosomes of functional content and genetic variation, which might limit their potential to generate adaptive diversity. Males of the freshwater fish Poecilia parae occur as one of five discrete morphs, all of which shoal together in natural populations where morph frequency has been stable for over 50 years. Each morph uses a different complex reproductive strategy and morphs differ dramatically in colour, body size and mating behaviour. Morph phenotype is passed perfectly from father to son, indicating there are five Y haplotypes segregating in the species, which encode the complex male morph characteristics. Here, we examine Y diversity in natural populations of P. parae. Using linked-read sequencing on multiple P. parae females and males of all five morphs, we find that the genetic architecture of the male morphs evolved on the Y chromosome after recombination suppression had occurred with the X. Comparing Y chromosomes between each of the morphs, we show that, although the Ys of the three minor morphs that differ in colour are highly similar, there are substantial amounts of unique genetic material and divergence between the Ys of the three major morphs that differ in reproductive strategy, body size and mating behaviour. Altogether, our results suggest that the Y chromosome is able to overcome the constraints of recombination loss to generate extreme diversity, resulting in five discrete Y chromosomes that control complex reproductive strategies.

Leveraging Short-Read Sequencing to Explore the Genomics of Sepiolid Squid

Due to their large size (∼3-5 Gb) and high repetitive content, the study of cephalopod genomes has historically been problematic. However, with the recent sequencing of several cephalopod genomes, including the Hawaiian bobtail squid (Euprymna scolopes), whole-genome studies of these molluscs are now possible. Of particular interest are the sepiolid or bobtail squids, many of which develop photophores in which bioluminescent bacterial symbionts reside. The variable presence of the symbiosis throughout the family allows us to determine regions of the genome that are under selection in symbiotic lineages, potentially providing a mechanism for identifying genes instrumental in the evolution of these mutualistic associations. To this end, we have used high-throughput sequencing to generate sequence from five bobtail squid genomes, four of which maintain symbioses with luminescent bacteria (E. hyllebergi, E. albatrossae, E. scolopes and Rondeletiola minor), and one of which does not (Sepietta neglecta). When we performed K-mer based heterozygosity and genome size estimations, we found that the Euprymna genus has a higher predicted genome size than other bobtail squid (∼ 5 Gb as compared to ∼ 4 Gb) and lower genomic heterozygosity. When we analyzed the repetitive content of the genomes, we found that genomes in the genus Euprymna appear to have recently acquired a significant quantity of LINE elements that are not found in its sister genus Rondeletiola or the closely related Sepietta. Using Abyss-2.0 and then Chromosomer with the published E. scolopes genome as a reference, we generated E. hyllebergi and E. albatrossae genomes of 1.54-1.57 Gb in size, but containing over 78-81% of eukaryotic single-copy othologs. The data we have generated will enable future whole-genome comparisons between these species to determine gene and regulatory content that differs between symbiotic and non-symbiotic lineages, as well as genes associated with symbiosis that are under selection.

The IS630/Tc1/mariner transposons in three ctenophore genomes

Transposable elements (TEs) exert a significant effect on the structure and functioning of the genomes and also serve as a source of the new genes. The study of the TE diversity and evolution in different taxa is indispensable for the fundamental understanding of their roles in the genomes. IS630/Tc1/mariner (ITm) transposable elements represent the most prevalent and diverse group of DNA transposons. In this work, we studied the diversity, evolutionary dynamics and the phylogenetic relationships of the ITm transposons found in three ctenophore species: Mnemiopsis leidyi, Pleurobrachia bachei, Beroe ovata. We identified 29 ITm transposons, seven of which possess the terminal inverted repeats (TIRs) and an intact transposase, and, thus, are, presumably, active. Four other ITm transposons have the features of domesticated TEs. According to the results of the phylogenetic analysis, the ITm transposons of the ctenophores represent five groups - MLE/DD34D, TLE/DD34-38E, mosquito/DD37E, Visiror/DD41D and pogo/DDxD. Pogo/DDxD superfamily turnes out to be the most diverse and prevalent, since it accounts for more than 40% of the TEs identified. The data obtained in this research will fill the gap of knowledge of the diversity and evolution of the ITm transposons in the multicellular genomes and will lay the ground for the study of the TE effects on the evolution of the ctenophores.

Movement of transposable elements contributes to cichlid diversity

African cichlid fishes are a prime model for studying speciation mechanisms. Despite the development of extensive genomic resources, it has been difficult to determine which sources of genetic variation are responsible for cichlid phenotypic variation. One of their most variable phenotypes is visual sensitivity, with some of the largest spectral shifts among vertebrates. These shifts arise primarily from differential expression of seven cone opsin genes. By mapping expression quantitative trait loci (eQTL) in intergeneric crosses of Lake Malawi cichlids, we previously identified four causative genetic variants that correspond to indels in the promoters of either key transcription factors or an opsin gene. In this comprehensive study, we show that these indels are the result of the movement of transposable elements (TEs) that correlate with opsin expression variation across the Malawi flock. In tracking the evolutionary history of these particular indels, we found they are endemic to Lake Malawi, suggesting that these TEs are recently active and are segregating within the Malawi cichlid lineage. However, an independent indel has arisen at a similar genomic location in one locus outside of the Malawi flock. The convergence in TE movement suggests these loci are primed for TE insertion and subsequent deletions. Increased TE mobility may be associated with interspecific hybridization, which disrupts mechanisms of TE suppression. This might provide a link between cichlid hybridization and accelerated regulatory variation. Overall, our study suggests that TEs may be an important driver of key regulatory changes, facilitating rapid phenotypic change and possibly speciation in African cichlids.

Mobile Elements in Ray-Finned Fish Genomes

Ray-finned fishes (Actinopterygii) are a very diverse group of vertebrates, encompassing species adapted to live in freshwater and marine environments, from the deep sea to high mountain streams. Genome sequencing offers a genetic resource for investigating the molecular bases of this phenotypic diversity and these adaptations to various habitats. The wide range of genome sizes observed in fishes is due to the role of transposable elements (TEs), which are powerful drivers of species diversity. Analyses performed to date provide evidence that class II DNA transposons are the most abundant component in most fish genomes and that compared to other vertebrate genomes, many TE superfamilies are present in actinopterygians. Moreover, specific TEs have been reported in ray-finned fishes as a possible result of an intricate relationship between TE evolution and the environment. The data summarized here underline the biological interest in Actinopterygii as a model group to investigate the mechanisms responsible for the high biodiversity observed in this taxon.

Adaptation to Extreme Antarctic Environments Revealed by the Genome of a Sea Ice Green Alga

The unicellular green alga Chlamydomonas sp. ICE-L thrives in polar sea ice, where it tolerates extreme low temperatures, high salinity, and broad seasonal fluctuations in light conditions. Despite the high interest in biotechnological uses of this species, little is known about the adaptations that allow it to thrive in this harsh and complex environment. Here, we assembled a high-quality genome sequence of ∼542 Mb and found that retrotransposon proliferation contributed to the relatively large genome size of ICE-L when compared to other chlorophytes. Genomic features that may support the extremophilic lifestyle of this sea ice alga include massively expanded gene families involved in unsaturated fatty acid biosynthesis, DNA repair, photoprotection, ionic homeostasis, osmotic homeostasis, and reactive oxygen species detoxification. The acquisition of multiple ice binding proteins through putative horizontal gene transfer likely contributed to the origin of the psychrophilic lifestyle in ICE-L. Additional innovations include the significant upregulation under abiotic stress of several expanded ICE-L gene families, likely reflecting adaptive changes among diverse metabolic processes. Our analyses of the genome, transcriptome, and functional assays advance general understanding of the Antarctic green algae and offer potential explanations for how green plants adapt to extreme environments.

Multiple independent chromosomal fusions accompanied the radiation of the Antarctic teleost genus Trematomus (Notothenioidei:Nototheniidae)

Background

Chromosomal rearrangements are thought to be an important driving force underlying lineage diversification, but their link to speciation continues to be debated. Antarctic teleost fish of the family Nototheniidae (Notothenioidei) diversified in a changing environmental context, which led to ecological, morphological, and genetic differentiation among populations. In addition, extensive chromosomal repatterning accompanied species divergence in several clades. The most striking karyotypic changes involved the recent species radiation (about 10 My) of the genus Trematomus, with chromosomal pair numbers ranging between 29 and 12. These dramatic reductions in chromosome number resulted mostly from large-scale chromosome fusions. Multiple centric and/or tandem fusions have been hypothesized in at least seven of the twelve recognized Trematomus species. To reconstruct their evolutionary history, we employed comparative cytogenomics (BAC-FISH and chromosome painting) to reveal patterns of interspecific chromosomal orthologies across several notothenioid clades.

Results

We defined orthologous chromosomal segments of reference, termed Structural Units (SUs). SUs were identified in a total of 18 notothenioid species. We demonstrated for the first time that SUs were strongly conserved across every specimen examined, with chromosomal syntenies highlighting a paucity of intrachromosomal macro-rearrangements. Multiple independent fusions of these SUs were inferred in the Trematomus species, in contrast to the shared SU fusions in species of the sister lineage Notothenia.

Conclusions

The SU segments were defined units of chromosomal rearrangement in the entire family Nototheiidae, which diverged from the other notothenioid families 20 My ago. Some of the identified chromosomal syntenies within the SUs were even conserved in their closest relatives, the family Eleginopsidae. Comparing the timing of acquisition of the fusions in the closely related genera Notothenia and Trematomus of the nototheniid species family, we conclude that they exhibit distinct chromosomal evolutionary histories, which may be relevant to different speciation scenarios.

The genomic basis for colonizing the freezing Southern Ocean revealed by Antarctic toothfish and Patagonian robalo genomes

Background

The Southern Ocean is the coldest ocean on Earth but a hot spot of evolution. The bottom-dwelling Eocene ancestor of Antarctic notothenioid fishes survived polar marine glaciation and underwent adaptive radiation, forming >120 species that fill all water column niches today. Genome-wide changes enabling physiological adaptations and the rapid expansion of the Antarctic notothenioids remain poorly understood.

Results

We sequenced and compared 2 notothenioid genomes—the cold-adapted and neutrally buoyant Antarctic toothfish Dissostichus mawsoni and the basal Patagonian robalo Eleginops maclovinus, representing the temperate ancestor. We detected >200 protein gene families that had expanded and thousands of genes that had evolved faster in the toothfish, with diverse cold-relevant functions including stress response, lipid metabolism, protein homeostasis, and freeze resistance. Besides antifreeze glycoprotein, an eggshell protein had functionally diversified to aid in cellular freezing resistance. Genomic and transcriptomic comparisons revealed proliferation of selcys–transfer RNA genes and broad transcriptional upregulation across anti-oxidative selenoproteins, signifying their prominent role in mitigating oxidative stress in the oxygen-rich Southern Ocean. We found expansion of transposable elements, temporally correlated to Antarctic notothenioid diversification. Additionally, the toothfish exhibited remarkable shifts in genetic programs towards enhanced fat cell differentiation and lipid storage, and promotion of chondrogenesis while inhibiting osteogenesis in bone development, collectively contributing to the achievement of neutral buoyancy and pelagicism.

Conclusions

Our study revealed a comprehensive landscape of evolutionary changes essential for Antarctic notothenioid cold adaptation and ecological expansion. The 2 genomes are valuable resources for further exploration of mechanisms underlying the spectacular notothenioid radiation in the coldest marine environment.

Evolution and diversity of transposable elements in fish genomes

Transposable elements (TEs) are genomic sequences that can move, multiply, and often form sizable fractions of vertebrate genomes. Fish belong to a unique group of vertebrates, since their karyotypes and genome sizes are more diverse and complex, with probably higher diversity and evolution specificity of TE. To investigate the characteristics of fish TEs, we compared the mobilomes of 39 species, and observed significant variation of TE content in fish (from 5% in pufferfish to 56% in zebrafish), along with a positive correlation between fish genome size and TE content. In different classification hierarchies, retrotransposons (class), long terminal repeat (order), as well as Helitron, Maverick, Kolobok, CMC, DIRS, P, I, L1, L2, and 5S (superfamily) were all positively correlated with fish genome size. Consistent with previous studies, our data suggested fish genomes to not always be dominated by DNA transposons; long interspersed nuclear elements are also prominent in many species. This study suggests CR1 distribution in fish genomes to be obviously regular, and provides new clues concerning important events in vertebrate evolution. Altogether, our results highlight the importance of TEs in the structure and evolution of fish genomes and suggest fish species diversity to parallel transposon content diversification.

The Chromosome-Level Genome Assembly of European Grayling Reveals Aspects of a Unique Genome Evolution Process Within Salmonids

Salmonids represent an intriguing taxonomical group for investigating genome evolution in vertebrates due to their relatively recent last common whole genome duplication event, which occurred between 80 and 100 million years ago. Here, we report on the chromosome-level genome assembly of European grayling (Thymallus thymallus), which represents one of the earliest diverged salmonid subfamilies. To achieve this, we first generated relatively long genomic scaffolds by using a previously published draft genome assembly along with long-read sequencing data and a linkage map. We then merged those scaffolds by applying synteny evidence from the Atlantic salmon (Salmo salar) genome. Comparisons of the European grayling genome assembly to the genomes of Atlantic salmon and Northern pike (Esox lucius), the latter used as a nonduplicated outgroup, detailed aspects of the characteristic chromosome evolution process that has taken place in European grayling. While Atlantic salmon and other salmonid genomes are portrayed by the typical occurrence of numerous chromosomal fusions, European grayling chromosomes were confirmed to be fusion-free and were characterized by a relatively large proportion of paracentric and pericentric inversions. We further reported on transposable elements specific to either the European grayling or Atlantic salmon genome, on the male-specific sdY gene in the European grayling chromosome 11A, and on regions under residual tetrasomy in the homeologous European grayling chromosome pairs 9A-9B and 25A-25B. The same chromosome pairs have been observed under residual tetrasomy in Atlantic salmon and in other salmonids, suggesting that this feature has been conserved since the subfamily split.

Genome of Crucihimalaya himalaica, a close relative of Arabidopsis, shows ecological adaptation to high altitude

Crucihimalaya himalaica, a close relative of Arabidopsis and Capsella, grows on the Qinghai-Tibet Plateau (QTP) about 4,000 m above sea level and represents an attractive model system for studying speciation and ecological adaptation in extreme environments. We assembled a draft genome sequence of 234.72 Mb encoding 27,019 genes and investigated its origin and adaptive evolutionary mechanisms. Phylogenomic analyses based on 4,586 single-copy genes revealed that C. himalaica is most closely related to Capsella (estimated divergence 8.8 to 12.2 Mya), whereas both species form a sister clade to Arabidopsis thaliana and Arabidopsis lyrata, from which they diverged between 12.7 and 17.2 Mya. LTR retrotransposons in C. himalaica proliferated shortly after the dramatic uplift and climatic change of the Himalayas from the Late Pliocene to Pleistocene. Compared with closely related species, C. himalaica showed significant contraction and pseudogenization in gene families associated with disease resistance and also significant expansion in gene families associated with ubiquitin-mediated proteolysis and DNA repair. We identified hundreds of genes involved in DNA repair, ubiquitin-mediated proteolysis, and reproductive processes with signs of positive selection. Gene families showing dramatic changes in size and genes showing signs of positive selection are likely candidates for C. himalaica's adaptation to intense radiation, low temperature, and pathogen-depauperate environments in the QTP. Loss of function at the S-locus, the reason for the transition to self-fertilization of C. himalaica, might have enabled its QTP occupation. Overall, the genome sequence of C. himalaica provides insights into the mechanisms of plant adaptation to extreme environments.

Antarctic blackfin icefish genome reveals adaptations to extreme environments

Icefishes (suborder Notothenioidei; family Channichthyidae) are the only vertebrates that lack functional haemoglobin genes and red blood cells. Here, we report a high-quality genome assembly and linkage map for the Antarctic blackfin icefish Chaenocephalus aceratus, highlighting evolved genomic features for its unique physiology. Phylogenomic analysis revealed that Antarctic fish of the teleost suborder Notothenioidei, including icefishes, diverged from the stickleback lineage about 77 million years ago and subsequently evolved cold-adapted phenotypes as the Southern Ocean cooled to sub-zero temperatures. Our results show that genes involved in protection from ice damage, including genes encoding antifreeze glycoprotein and zona pellucida proteins, are highly expanded in the icefish genome. Furthermore, genes that encode enzymes that help to control cellular redox state, including members of the sod3 and nqo1 gene families, are expanded, probably as evolutionary adaptations to the relatively high concentration of oxygen dissolved in cold Antarctic waters. In contrast, some crucial regulators of circadian homeostasis (cry and per genes) are absent from the icefish genome, suggesting compromised control of biological rhythms in the polar light environment. The availability of the icefish genome sequence will accelerate our understanding of adaptation to extreme Antarctic environments.

Insertion Hot Spots of DIRS1 Retrotransposon and Chromosomal Diversifications among the Antarctic Teleosts Nototheniidae

By their faculty to transpose, transposable elements are known to play a key role in eukaryote genomes, impacting both their structuration and remodeling. Their integration in targeted sites may lead to recombination mechanisms involved in chromosomal rearrangements. The Antarctic fish family Nototheniidae went through several waves of species radiations. It is a suitable model to study transposable element (TE)-mediated mechanisms associated to genome and chromosomal diversifications. After the characterization of Gypsy (GyNoto), Copia (CoNoto), and DIRS1 (YNoto) retrotransposons in the genomes of Nototheniidae (diversity, distribution, conservation), we focused on their chromosome location with an emphasis on the three identified nototheniid radiations (the Trematomus, the plunderfishes, and the icefishes). The strong intrafamily TE conservation and wide distribution across species of the whole family suggest an ancestral acquisition with potential secondary losses in some lineages. GyNoto and CoNoto (including Hydra and GalEa clades) mostly produced interspersed signals along chromosomal arms. On the contrary, insertion hot spots accumulating in localized regions (mainly next to centromeric and pericentromeric regions) highlighted the potential role of YNoto in chromosomal diversifications as facilitator of the fusions which occurred in many nototheniid lineages, but not of the fissions.