An Annotated Chromosome-Level Reference Genome of the Red-Eared Slider Turtle (Trachemys scripta elegans)

slc26a12-A novel member of the slc26 family, is located in tandem with slc26a2 in coelacanths, amphibians, reptiles, and birds

Solute carrier family 26 (Slc26) is a family of anion exchangers with 11 members in mammals (named Slc26a1-a11). Here, we identified a novel member of the slc26 family, slc26a12, located in tandem with slc26a2 in the genomes of several vertebrate lineages. BLAST and synteny analyses of various jawed vertebrate genome databases revealed that slc26a12 is present in coelacanths, amphibians, reptiles, and birds but not in cartilaginous fishes, lungfish, mammals, or ray-finned fishes. In some avian and reptilian lineages such as owls, penguins, egrets, and ducks, and most turtles examined, slc26a12 was lost or pseudogenized. Phylogenetic analysis showed that Slc26a12 formed an independent branch with the other Slc26 members and Slc26a12, Slc26a1 and Slc26a2 formed a single branch, suggesting that these three members formed a subfamily in Slc26. In jawless fish, hagfish have two genes homologous to slc26a2 and slc26a12, whereas lamprey has a single gene homologous to slc26a2. African clawed frogs express slc26a12 in larval gills, skin, and fins. These results show that slc26a12 was present at least before the separation of lobe-finned fish and tetrapods; the name slc26a12 is appropriate because the gene duplication occurred in the distant past.

A chromosome-level genome assembly of the pig-nosed turtle (Carettochelys insculpta)

The pig-nosed turtle (Carettochelys insculpta) represents the only extant species within the Carettochelyidae family, is a unique Trionychia member fully adapted to aquatic life and currently facing endangerment. To enhance our understanding of this species and contribute to its conservation efforts, we employed high-fidelity (HiFi) and Hi-C sequencing technology to generate its genome assembly at the chromosome level. The assembly result spans 2.18 Gb, with a contig N50 of 126 Mb, encompassing 34 chromosomes that account for 99.6% of the genome. The assembly has a BUSCO score above 95% with different databases and strong collinearity with Yangtze giant softshell turtles (Rafetus swinhoei), indicating its completeness and continuity. A total of 19,175 genes and 46.86% repetitive sequences were annotated. The availability of this chromosome-scale genome represents a valuable resource for the pig-nosed turtle, providing insights into its aquatic adaptation and serving as a foundation for future turtle research.

Experiencing short heat waves early in development changes thermal responsiveness of turtle embryos to later heat waves

Although physiological responses to the thermal environment are most frequently investigated using constant temperatures, the incorporation of thermal variability can allow for a more accurate prediction of how thermally sensitive species respond to a rapidly changing climate. In species with temperature-dependent sex determination (TSD), developmental responses to incubation temperature are mediated by several genes involved in gonadal differentiation. Kdm6b and Dmrt1 respond to cool incubation temperatures and are associated with testis development, while FoxL2 and Cyp19A1 respond to warm incubation temperatures and are associated with ovary development. Using fluctuating incubation temperatures, we designed two studies, one investigating how conflicting thermal cues affect the timing of commitment to gonadal development, and another investigating the rapid molecular responses to conflicting thermal cues in the red-eared slider turtle (Trachemys scripta). Using gene expression as a proxy of timing of commitment to gonadal fate, results from the first study show that exposure to high amounts of conflicting thermal cues during development delays commitment to gonadal fate. Results from the second study show that Kdm6b splice variants exhibit differential responses to early heat wave exposure, but rapidly (within 2 days) recover to pre-exposure levels after the heat wave. Despite changes in the expression of Kdm6b splice variants, there was no effect on Dmrt1 expression. Collectively, these findings demonstrate how short exposures to heat early in development can change how embryos respond to heat later in development.

The first high-quality chromosome-level genome of Eretmochelys imbricata using HiFi and Hi-C data

Eretmochelys imbricata, a critically endangered sea turtle inhabiting tropical oceans and protected across the world, had an unknown genome sequence until now. In this study, we used HiFi reads and Hi-C technology to assemble a high-quality, chromosome-level genome of E. imbricata. The genome size was 2,138.26 Mb, with contig N50 length of 123.49 Mb and scaffold N50 of 137.21 Mb. Approximately 97.52% of the genome sequence was anchored to 28 chromosomes. A total of 20,206 protein-coding genes were predicted. We also analyzed the evolutionary relationships, gene family expansions, and positive selection of E. imbricata. Our results revealed that E. imbricata diverged from Chelonia mydas 38 million years ago and had enriched olfactory receptors and aging-related genes. Our genome will be useful for studying E. imbricata and its conservation.

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.

Genome assembly and isoform analysis of a highly heterozygous New Zealand fisheries species, the tarakihi (Nemadactylus macropterus)

Although being some of the most valuable and heavily exploited wild organisms, few fisheries species have been studied at the whole-genome level. This is especially the case in New Zealand, where genomics resources are urgently needed to assist fisheries management. Here, we generated 55 Gb of short Illumina reads (92× coverage) and 73 Gb of long Nanopore reads (122×) to produce the first genome assembly of the marine teleost tarakihi [Nemadactylus macropterus (Forster, 1801)], a highly valuable fisheries species in New Zealand. An additional 300 Mb of Iso-Seq reads were obtained to assist in gene annotation. The final genome assembly was 568 Mb long with an N50 of 3.37 Mb. The genome completeness was high, with 97.8% of complete Actinopterygii Benchmarking Universal Single-Copy Orthologs. Heterozygosity values estimated through k-mer counting (1.00%) and bi-allelic SNPs (0.64%) were high compared with the same values reported for other fishes. Iso-Seq analysis recovered 91,313 unique transcripts from 15,515 genes (mean ratio of 5.89 transcripts per gene), and the most common alternative splicing event was intron retention. This highly contiguous genome assembly and the isoform-resolved transcriptome will provide a useful resource to assist the study of population genomics and comparative eco-evolutionary studies in teleosts and related organisms.

Chromosome-level genome assembly for the Aldabra giant tortoise enables insights into the genetic health of a threatened population

Background

The Aldabra giant tortoise (Aldabrachelys gigantea) is one of only two giant tortoise species left in the world. The species is endemic to Aldabra Atoll in Seychelles and is listed as Vulnerable on the International Union for Conservation of Nature Red List (v2.3) due to its limited distribution and threats posed by climate change. Genomic resources for A. gigantea are lacking, hampering conservation efforts for both wild and ex situpopulations. A high-quality genome would also open avenues to investigate the genetic basis of the species’ exceptionally long life span.

Findings

We produced the first chromosome-level de novo genome assembly of A. gigantea using PacBio High-Fidelity sequencing and high-throughput chromosome conformation capture. We produced a 2.37-Gbp assembly with a scaffold N50 of 148.6 Mbp and a resolution into 26 chromosomes. RNA sequencing–assisted gene model prediction identified 23,953 protein-coding genes and 1.1 Gbp of repetitive sequences. Synteny analyses among turtle genomes revealed high levels of chromosomal collinearity even among distantly related taxa. To assess the utility of the high-quality assembly for species conservation, we performed a low-coverage resequencing of 30 individuals from wild populations and two zoo individuals. Our genome-wide population structure analyses detected genetic population structure in the wild and identified the most likely origin of the zoo-housed individuals. We further identified putatively deleterious mutations to be monitored.

Conclusions

We establish a high-quality chromosome-level reference genome for A. gigantea and one of the most complete turtle genomes available. We show that low-coverage whole-genome resequencing, for which alignment to the reference genome is a necessity, is a powerful tool to assess the population structure of the wild population and reveal the geographic origins of ex situ individuals relevant for genetic diversity management and rewilding efforts.

Transcriptome Analyses Provide Insights into the Auditory Function in Trachemys scripta elegans

An auditory ability is essential for communication in vertebrates, and considerable attention has been paid to auditory sensitivity in mammals, birds, and frogs. Turtles were thought to be deaf for a long time; however, recent studies have confirmed the presence of an auditory ability in Trachemys scripta elegans as well as sex-related differences in hearing sensitivity. Earlier studies mainly focused on the morphological and physiological functions of the hearing organ in turtles; thus, the gene expression patterns remain unclear. In this study, 36 transcriptomes from six tissues (inner ear, tympanic membrane, brain, eye, lung, and muscle) were sequenced to explore the gene expression patterns of the hearing system in T. scripta elegans. A weighted gene co-expression network analysis revealed that hub genes related to the inner ear and tympanic membrane are involved in development and signal transduction. Moreover, we identified six differently expressed genes (GABRA1, GABRG2, GABBR2, GNAO1, SLC38A1, and SLC12A5) related to the GABAergic synapse pathway as candidate genes to explain the differences in sexually dimorphic hearing sensitivity. Collectively, this study provides a critical foundation for genetic research on auditory functions in turtles.

Effects of Constant and Fluctuating Temperatures on Gene Expression During Gonadal Development

There is ample research demonstrating that temperature can have complex effects on biological processes, including the timing of when organisms respond to temperature; some responses occur rapidly while others require an extended exposure time. However, most of what we know about temperature effects comes from studies using constant temperature conditions, which are not reflective of natural, fluctuating temperatures. Species with temperature-dependent sex determination (TSD) present an ideal system to study the temporal aspects of the temperature response because prior research has established a number of temperature-responsive genes involved in TSD, albeit under constant temperatures. To investigate potential differences in timing of sexual development between constant and fluctuating incubation temperatures, we exposed Trachemys scripta embryos to two conditions that produce males (constant 26 °C and 26 ± 3 °C) and two that produce females (constant 31 °C and 31 ± 3 °C) and sampled embryonic gonads for gene expression analysis via qPCR. We analyzed three genes involved in testis differentiation (Kdm6b, Dmrt1, and Sox9) and two genes involved in ovary differentiation (Foxl2 and Cyp19A1). Results show that Kdm6b expression was significantly lower under fluctuating temperatures compared to constant temperatures. Foxl2 and Cyp19A1 expression were also lower under fluctuating temperatures, but not at all stages of development. These results suggest that constant temperatures caused increases in both Foxl2 and Cyp19A1 expression earlier (developmental stage 20) than fluctuating temperatures (stages 22-23). Dmrt1 and Sox9 expression did not differ between constant and fluctuating temperatures. These results highlight that not all genes in a temperature-dependent process respond to temperature in the same manner. Whether there are functional consequences of this variation remains to be determined.

A chromosome-level genome assembly and annotation of the desert horned lizard, Phrynosoma platyrhinos, provides insight into chromosomal rearrangements among reptiles

BACKGROUND

The increasing number of chromosome-level genome assemblies has advanced our knowledge and understanding of macroevolutionary processes. Here, we introduce the genome of the desert horned lizard, Phrynosoma platyrhinos, an iguanid lizard occupying extreme desert conditions of the American southwest. We conduct analysis of the chromosomal structure and composition of this species and compare these features across genomes of 12 other reptiles (5 species of lizards, 3 snakes, 3 turtles, and 1 bird).

FINDINGS

The desert horned lizard genome was sequenced using Illumina paired-end reads and assembled and scaffolded using Dovetail Genomics Hi-C and Chicago long-range contact data. The resulting genome assembly has a total length of 1,901.85 Mb, scaffold N50 length of 273.213 Mb, and includes 5,294 scaffolds. The chromosome-level assembly is composed of 6 macrochromosomes and 11 microchromosomes. A total of 20,764 genes were annotated in the assembly. GC content and gene density are higher for microchromosomes than macrochromosomes, while repeat element distributions show the opposite trend. Pathway analyses provide preliminary evidence that microchromosome and macrochromosome gene content are functionally distinct. Synteny analysis indicates that large microchromosome blocks are conserved among closely related species, whereas macrochromosomes show evidence of frequent fusion and fission events among reptiles, even between closely related species.

CONCLUSIONS

Our results demonstrate dynamic karyotypic evolution across Reptilia, with frequent inferred splits, fusions, and rearrangements that have resulted in shuffling of chromosomal blocks between macrochromosomes and microchromosomes. Our analyses also provide new evidence for distinct gene content and chromosomal structure between microchromosomes and macrochromosomes within reptiles.

Chromosome-level genome assembly of the Chinese three-keeled pond turtle (Mauremys reevesii) provides insights into freshwater adaptation

Mauremys reevesii is an endangered freshwater turtle that symbolizes longevity in Chinese culture. Despite its importance, genetic studies of this species remain limited, with no genomic sequence reported to date. Here, we report a high-quality, chromosome-level genomic sequence of M. reevesii obtained using a combination of Nanopore and Hi-C sequencing technologies. The 2.37 Gb M. reevesii genome was assembled from a total of ~226.80 Gb of Nanopore sequencing data. The M. reevesii genome contig N50 is 34.73 Mb, the highest value in published turtle genomes. In total, 18,238 genes were functionally annotated. The contigs were clustered and ordered onto 27 pseudochromosomes covering ~96.55% of the genome assembled with Hi-C data. To explore genome evolution, synteny analysis was performed between M. reevesii (freshwater turtle) and Gopherus evgoodei (terrestrial turtle) genomes. In general, each chromosome of M. reevesii corresponded to one chromosome of Gopherus evgoodei, but some interchromosomal rearrangements occurred between the two species based on the assembled genomes. These interchromosomal rearrangements were further confirmed by mapping of the long-read nanopore data to the assembly. The reconstructed demographic history showed varied effective population size among freshwater, marine and terrestrial turtles. We also discovered expansion of genes related to the innate immune system in M. reevesii that may provide defence against freshwater pathogens. The high-quality genomic sequence provides a valuable genetic resource for further studies of genetics and genome evolution in turtles.

ClineHelpR: an R package for genomic cline outlier detection and visualization

Background

Patterns of multi-locus differentiation (i.e., genomic clines) often extend broadly across hybrid zones and their quantification can help diagnose how species boundaries are shaped by adaptive processes, both intrinsic and extrinsic. In this sense, the transitioning of loci across admixed individuals can be contrasted as a function of the genome-wide trend, in turn allowing an expansion of clinal theory across a much wider array of biodiversity. However, computational tools that serve to interpret and consequently visualize ‘genomic clines’ are limited, and users must often write custom, relatively complex code to do so.

Results

Here, we introduce the ClineHelpR R-package for visualizing genomic clines and detecting outlier loci using output generated by two popular software packages, bgc and Introgress. ClineHelpR bundles both input generation (i.e., filtering datasets and creating specialized file formats) and output processing (e.g., MCMC thinning and burn-in) with functions that directly facilitate interpretation and hypothesis testing. Tools are also provided for post-hoc analyses that interface with external packages such as ENMeval and RIdeogram.

Conclusions

Our package increases the reproducibility and accessibility of genomic cline methods, thus allowing an expanded user base and promoting these methods as mechanisms to address diverse evolutionary questions in both model and non-model organisms. Furthermore, the ClineHelpR extended functionality can evaluate genomic clines in the context of spatial and environmental features, allowing users to explore underlying processes potentially contributing to the observed patterns and helping facilitate effective conservation management strategies.

Sex chromosome evolution among amniotes: is the origin of sex chromosomes non-random?

Sex chromosomes are a great example of a convergent evolution at the genomic level, having evolved dozens of times just within amniotes. An intriguing question is whether this repeated evolution was random, or whether some ancestral syntenic blocks have significantly higher chance to be co-opted for the role of sex chromosomes owing to their gene content related to gonad development. Here, we summarize current knowledge on the evolutionary history of sex determination and sex chromosomes in amniotes and evaluate the hypothesis of non-random emergence of sex chromosomes. The current data on the origin of sex chromosomes in amniotes suggest that their evolution is indeed non-random. However, this non-random pattern is not very strong, and many syntenic blocks representing putatively independently evolved sex chromosomes are unique. Still, repeatedly co-opted chromosomes are an excellent model system, as independent co-option of the same genomic region for the role of sex chromosome offers a great opportunity for testing evolutionary scenarios on the sex chromosome evolution under the explicit control for the genomic background and gene identity. Future studies should use these systems more to explore the convergent/divergent evolution of sex chromosomes. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

The Easter Egg Weevil (Pachyrhynchus) genome reveals syntenic patterns in Coleoptera across 200 million years of evolution

Patterns of genomic architecture across insects remain largely undocumented or decoupled from a broader phylogenetic context. For instance, it is unknown whether translocation rates differ between insect orders. We address broad scale patterns of genome architecture across Insecta by examining synteny in a phylogenetic framework from open-source insect genomes. To accomplish this, we add a chromosome level genome to a crucial lineage, Coleoptera. Our assembly of the Pachyrhynchus sulphureomaculatus genome is the first chromosome scale genome for the hyperdiverse Phytophaga lineage and currently the largest insect genome assembled to this scale. The genome is significantly larger than those of other weevils, and this increase in size is caused by repetitive elements. Our results also indicate that, among beetles, there are instances of long-lasting (>200 Ma) localization of genes to a particular chromosome with few translocation events. While some chromosomes have a paucity of translocations, intra-chromosomal synteny was almost absent, with gene order thoroughly shuffled along a chromosome. This large amount of reshuffling within chromosomes with few inter-chromosomal events contrasts with patterns seen in mammals in which the chromosomes tend to exchange larger blocks of material more readily. To place our findings in an evolutionary context, we compared syntenic patterns across Insecta in a phylogenetic framework. For the first time, we find that synteny decays at an exponential rate relative to phylogenetic distance. Additionally, there are significant differences in decay rates between insect orders, this pattern was not driven by Lepidoptera alone which has a substantially different rate.

Brief exposure to warm temperatures reduces intron retention in Kdm6b in a species with temperature-dependent sex determination

Animals with temperature-dependent sex determination (TSD) respond to thermal cues during early embryonic development to trigger gonadal differentiation. TSD has primarily been studied using constant temperature incubations, where embryos are exposed to constant male- or female-producing temperatures, and these studies have identified genes that display sex-specific expression in response to incubation temperature. Kdm6b, a histone demethylase gene, has received specific attention as it is among the initial genes to respond to incubation temperature and is necessary for testis development. Interestingly, Kdm6b retains an intron when eggs are incubated at a constant male-producing temperature, but the role of thermal variability in this developmental process is relatively understudied. Species with TSD regularly experience thermal cues that fluctuate between male- and female-producing temperatures throughout development but it is unclear how Kdm6b responds to such variable temperatures. In this study, we investigate temperature-sensitive splicing in Kdm6b by exposing embryos to male- and female-producing thermal conditions. We show a rapid decrease in levels of the intron retaining transcript of Kdm6b upon exposure to female-producing conditions. These results demonstrate that, under ecologically relevant conditions, temperature-sensitive splicing can differentially regulate genes critical to TSD.

Predicting the Evolution of Syntenies—An Algorithmic Review

Syntenies are genomic segments of consecutive genes identified by a certain conservation in gene content and order. The notion of conservation may vary from one definition to another, the more constrained requiring identical gene contents and gene orders, while more relaxed definitions just require a certain similarity in gene content, and not necessarily in the same order. Regardless of the way they are identified, the goal is to characterize homologous genomic regions, i.e., regions deriving from a common ancestral region, reflecting a certain gene co-evolution that can enlighten important functional properties. In addition of being able to identify them, it is also necessary to infer the evolutionary history that has led from the ancestral segment to the extant ones. In this field, most algorithmic studies address the problem of inferring rearrangement scenarios explaining the disruption in gene order between segments with the same gene content, some of them extending the evolutionary model to gene insertion and deletion. However, syntenies also evolve through other events modifying their content in genes, such as duplications, losses or horizontal gene transfers, i.e., the movement of genes from one species to another. Although the reconciliation approach between a gene tree and a species tree addresses the problem of inferring such events for single-gene families, little effort has been dedicated to the generalization to segmental events and to syntenies. This paper reviews some of the main algorithmic methods for inferring ancestral syntenies and focus on those integrating both gene orders and gene trees.

New Ther1-derived SINE Squam3 in scaled reptiles

BACKGROUND

SINEs comprise a significant part of animal genomes and are used to study the evolution of diverse taxa. Despite significant advances in SINE studies in vertebrates and higher eukaryotes in general, their own evolution is poorly understood.

RESULTS

We have discovered and described in detail a new Squam3 SINE specific for scaled reptiles (Squamata). The subfamilies of this SINE demonstrate different distribution in the genomes of squamates, which together with the data on similar SINEs in the tuatara allowed us to propose a scenario of their evolution in the context of reptilian evolution.

CONCLUSIONS

Ancestral SINEs preserved in small numbers in most genomes can give rise to taxa-specific SINE families. Analysis of this aspect of SINEs can shed light on the history and mechanisms of SINE variation in reptilian genomes.

Cytogenetic Analysis of the Asian Box Turtles of the Genus Cuora (Testudines, Geoemydidae)

The Asian box turtle genus Cuora currently comprises 13 species with a wide distribution in Southeast Asia, including China and the islands of Indonesia and Philippines. The populations of these species are rapidly declining due to human pressure, including pollution, habitat loss, and harvesting for food consumption. Notably, the IUCN Red List identifies almost all species of the genus Cuora as Endangered (EN) or Critically Endangered (CR). In this study, we explore the karyotypes of 10 Cuora species with conventional (Giemsa staining, C-banding, karyogram reconstruction) and molecular cytogenetic methods (in situ hybridization with probes for rDNA loci and telomeric repeats). Our study reveals a diploid chromosome number of 2n = 52 chromosomes in all studied species, with karyotypes of similar chromosomal morphology. In all examined species, rDNA loci are detected at a single medium-sized chromosome pair and the telomeric repeats are restricted to the expected terminal position across all chromosomes. In contrast to a previous report, sex chromosomes are neither detected in Cuoragalbinifrons nor in any other species. Therefore, we assume that these turtles have either environmental sex determination or genotypic sex determination with poorly differentiated sex chromosomes. The conservation of genome organization could explain the numerous observed cases of interspecific hybridization both within the genus Cuora and across geoemydid turtles.

Karyotypic Evolution of Sauropsid Vertebrates Illuminated by Optical and Physical Mapping of the Painted Turtle and Slider Turtle Genomes

Recent sequencing and software enhancements have advanced our understanding of the evolution of genomic structure and function, especially addressing novel evolutionary biology questions. Yet fragmentary turtle genome assemblies remain a challenge to fully decipher the genetic architecture of adaptive evolution. Here, we use optical mapping to improve the contiguity of the painted turtle (Chrysemys picta) genome assembly and use de novo fluorescent in situ hybridization (FISH) of bacterial artificial chromosome (BAC) clones, BAC-FISH, to physically map the genomes of the painted and slider turtles (Trachemys scripta elegans). Optical mapping increased C. picta's N50 by ~242% compared to the previous assembly. Physical mapping permitted anchoring ~45% of the genome assembly, spanning 5544 genes (including 20 genes related to the sex determination network of turtles and vertebrates). BAC-FISH data revealed assembly errors in C. picta and T. s. elegans assemblies, highlighting the importance of molecular cytogenetic data to complement bioinformatic approaches. We also compared C. picta's anchored scaffolds to the genomes of other chelonians, chicken, lizards, and snake. Results revealed a mostly one-to-one correspondence between chromosomes of painted and slider turtles, and high homology among large syntenic blocks shared with other turtles and sauropsids. Yet, numerous chromosomal rearrangements were also evident across chelonians, between turtles and squamates, and between avian and non-avian reptiles.