Using low-coverage whole genome sequencing (genome skimming) to delineate three introgressed species of buffalofish (Ictiobus).

Consumption of buffalofish has been sporadically associated with Haff disease-like illnesses involving sudden onset muscle pain and weakness due to skeletal muscle rhabdomyolysis, but determination of precisely which species are associated with these illnesses has been impeded by a lack of species-specific DNA-based markers. Here, three closely related species of buffalofish native to the Mississippi River Basin (Ictiobus bubalus, Ictiobus cyprinellus and Ictiobus niger) that have previously proven genetically indistinguishable using both mitochondrial and nuclear single-locus sequencing were reliably discriminated using low-coverage whole genome sequencing ('genome skimming'). Using 44 specimens representing the three species collected from the mid/upper (Missouri) and lower (Louisiana) regions of the species' native ranges, the SISRS (Site Identification from Short Read Sequences) bioinformatics pipeline was adapted to (1) identify over 620Mbp of putatively homologous nuclear sequence data and (2) isolate over 140,000 single-nucleotide polymorphisms (SNPs) that supported accurate species delimitation, all without the use of a reference genome or annotation data. These sites were used to classify Ictiobus spp. samples with genome-skim data, along with a larger set (n = 67) where ultraconserved elements (UCEs) were sequenced. Analyses of whole mitochondrial data revealed more limited signal. Nearly all samples matched their purported species based on morphologic identification, but two Missouri samples morphologically identified as I. niger grouped with samples of I. bubalus, albeit with significant enrichment of I. niger SNPs. To our knowledge this is the first report of a DNA-based tool to reliably discriminate these three morphologically distinct species.

A Genomic Perspective on the Evolutionary Diversification of Turtles.

To examine phylogenetic heterogeneity in turtle evolution, we collected thousands of high-confidence single-copy orthologs from 19 genome assemblies representative of extant turtle diversity and estimated a phylogeny with multispecies coalescent and concatenated partitioned methods. We also collected next-generation sequences from 26 turtle species and assembled millions of biallelic markers to reconstruct phylogenies based on annotated regions from the western painted turtle (Chrysemys picta bellii) genome (coding regions, introns, untranslated regions, intergenic, and others). We then measured gene tree-species tree discordance, as well as gene and site heterogeneity at each node in the inferred trees, and tested for temporal patterns in phylogenomic conflict across turtle evolution. We found strong and consistent support for all bifurcations in the inferred turtle species phylogenies. However, a number of genes, sites, and genomic features supported alternate relationships between turtle taxa. Our results suggest that gene tree-species tree discordance in these data sets is likely driven by population-level processes such as incomplete lineage sorting. We found very little effect of substitutional saturation on species tree topologies, and no clear phylogenetic patterns in codon usage bias and compositional heterogeneity. There was no correlation between gene and site concordance, node age, and DNA substitution rate across most annotated genomic regions. Our study demonstrates that heterogeneity is to be expected even in well-resolved clades such as turtles, and that future phylogenomic studies should aim to sample as much of the genome as possible in order to obtain accurate phylogenies for assessing conservation priorities in turtles. [Discordance; genomes; phylogeny; turtles.].

Genome-Scale Profiling Reveals Noncoding Loci Carry Higher Proportions of Concordant Data.

Many evolutionary relationships remain controversial despite whole-genome sequencing data. These controversies arise, in part, due to challenges associated with accurately modeling the complex phylogenetic signal coming from genomic regions experiencing distinct evolutionary forces. Here, we examine how different regions of the genome support or contradict well-established relationships among three mammal groups using millions of orthologous parsimony-informative biallelic sites (PIBS) distributed across primate, rodent, and Pecora genomes. We compared PIBS concordance percentages among locus types (e.g. coding sequences (CDS), introns, intergenic regions), and contrasted PIBS utility over evolutionary timescales. Sites derived from noncoding sequences provided more data and proportionally more concordant sites compared with those from CDS in all clades. CDS PIBS were also predominant drivers of tree incongruence in two cases of topological conflict. PIBS derived from most locus types provided surprisingly consistent support for splitting events spread across the timescales we examined, although we find evidence that CDS and intronic PIBS may, respectively and to a limited degree, inform disproportionately about older and younger splits. In this era of accessible wholegenome sequence data, these results:1) suggest benefits to more intentionally focusing on noncoding loci as robust data for tree inference and 2) reinforce the importance of accurate modeling, especially when using CDS data.

Putative Independent Evolutionary Reversals from Genotypic to Temperature-Dependent Sex Determination are Associated with Accelerated Evolution of Sex-Determining Genes in Turtles.

The evolutionary lability of sex-determining mechanisms across the tree of life is well recognized, yet the extent of molecular changes that accompany these repeated transitions remain obscure. Most turtles retain the ancestral temperature-dependent sex determination (TSD) from which multiple transitions to genotypic sex determination (GSD) occurred independently, and two contrasting hypotheses posit the existence or absence of reversals back to TSD. Here we examined the molecular evolution of the coding regions of a set of gene regulators involved in gonadal development in turtles and several other vertebrates. We found slower molecular evolution in turtles and crocodilians compared to other vertebrates, but an acceleration in Trionychia turtles and at some phylogenetic branches demarcating major taxonomic diversification events. Of all gene classes examined, hormone signaling genes, and Srd5a1 in particular, evolve faster in many lineages and especially in turtles. Our data show that sex-linked genes do not follow a ubiquitous nor uniform pattern of molecular evolution. We then evaluated turtle nucleotide and protein evolution under two evolutionary hypotheses with or without GSD-to-TSD reversals, and found that when GSD-to-TSD reversals are considered, all transitional branches irrespective of direction, exhibit accelerated molecular evolution of nucleotide sequences, while GSD-to-TSD transitional branches also show acceleration in protein evolution. Significant changes in predicted secondary structure that may affect protein function were identified in three genes that exhibited hastened evolution in turtles compared to other vertebrates or in transitional versus non-transitional branches within turtles, rendering them candidates for a key role during SDM evolution in turtles.

Seeing red to being red: conserved genetic mechanism for red cone oil droplets and co-option for red coloration in birds and turtles.

Avian ketocarotenoid pigments occur in both the red retinal oil droplets that contribute to colour vision and bright red coloration used in signalling. Turtles are the only other tetrapods with red retinal oil droplets, and some also display red carotenoid-based coloration. Recently, the CYP2J19 gene was strongly implicated in ketocarotenoid synthesis in birds. Here, we investigate CYP2J19 evolution in relation to colour vision and red coloration in reptiles using genomic and expression data. We show that turtles, but not crocodiles or lepidosaurs, possess a CYP2J19 orthologue, which arose via gene duplication before turtles and archosaurs split, and which is strongly and specifically expressed in the ketocarotenoid-containing retina and red integument. We infer that CYP2J19 initially functioned in colour vision in archelosaurs and conclude that red ketocarotenoid-based coloration evolved independently in birds and turtles via gene regulatory changes of CYP2J19 Our results suggest that red oil droplets contributed to colour vision in dinosaurs and pterosaurs.

Cytogenetic Insights into the Evolution of Chromosomes and Sex Determination Reveal Striking Homology of Turtle Sex Chromosomes to Amphibian Autosomes.

Turtle karyotypes are highly conserved compared to other vertebrates; yet, variation in diploid number (2n = 26-68) reflects profound genomic reorganization, which correlates with evolutionary turnovers in sex determination. We evaluate the published literature and newly collected comparative cytogenetic data (G- and C-banding, 18S-NOR, and telomere-FISH mapping) from 13 species spanning 2n = 28-68 to revisit turtle genome evolution and sex determination. Interstitial telomeric sites were detected in multiple lineages that underwent diploid number and sex determination turnovers, suggesting chromosomal rearrangements. C-banding revealed potential interspecific variation in centromere composition and interstitial heterochromatin at secondary constrictions. 18S-NORs were detected in secondary constrictions in a single chromosomal pair per species, refuting previous reports of multiple NORs in turtles. 18S-NORs are linked to ZW chromosomes in Apalone and Pelodiscus and to X (not Y) in Staurotypus. Notably, comparative genomics across amniotes revealed that the sex chromosomes of several turtles, as well as mammals and some lizards, are homologous to components of Xenopus tropicalis XTR1 (carrying Dmrt1). Other turtle sex chromosomes are homologous to XTR4 (carrying Wt1). Interestingly, all known turtle sex chromosomes, except in Trionychidae, evolved via inversions around Dmrt1 or Wt1. Thus, XTR1 appears to represent an amniote proto-sex chromosome (perhaps linked ancestrally to XTR4) that gave rise to turtle and other amniote sex chromosomes.

Molecular cytogenetic search for cryptic sex chromosomes in painted turtles Chrysemys picta.

Sex determination is triggered by factors ranging from genotypic (GSD) to environmental (ESD), or both GSD + EE (GSD susceptible to environmental effects), and its evolution remains enigmatic. The presence/absence of sex chromosomes purportedly separates species at the ESD end of the continuum from the rest (GSD and GSD + EE) because the evolutionary dynamics of sex chromosomes and autosomes differ. However, studies suggest that turtles with temperature-dependent sex determination (TSD) are cryptically GSD and possess sex chromosomes. Here, we test this hypothesis in painted turtles Chrysemys picta (TSD), using comparative-genome-hybridization (CGH), a technique known to detect morphologically indistinguishable sex chromosomes in other turtles and reptiles. Our results show no evidence for the existence of sex chromosomes in painted turtles. While it remains plausible that cryptic sex chromosomes may exist in TSD turtles that are characterized by minor genetic differences that cannot be detected at the resolution of CGH, previous attempts have failed to identify sex-specific markers. Genomic sequencing should prove useful in providing conclusive evidence in this regard. If such efforts uncover sex chromosomes in TSD turtles, it may reveal the existence of a fundamental constraint for the evolution of a full spectrum of sex determination (from pure GSD to pure TSD) that is predicted theoretically. Finding sex chromosomes in ESD organisms would question whether pure ESD mechanisms exist at all in nature, or whether those systems currently considered pure ESD simply await the characterization of an underlying GSD architecture.

Molecular evolution of Dmrt1 accompanies change of sex-determining mechanisms in reptilia.

In reptiles, sex-determining mechanisms have evolved repeatedly and reversibly between genotypic and temperature-dependent sex determination. The gene Dmrt1 directs male determination in chicken (and presumably other birds), and regulates sex differentiation in animals as distantly related as fruit flies, nematodes and humans. Here, we show a consistent molecular difference in Dmrt1 between reptiles with genotypic and temperature-dependent sex determination. Among 34 non-avian reptiles, a convergently evolved pair of amino acids encoded by sequence within exon 2 near the DM-binding domain of Dmrt1 distinguishes species with either type of sex determination. We suggest that this amino acid shift accompanied the evolution of genotypic sex determination from an ancestral condition of temperature-dependent sex determination at least three times among reptiles, as evident in turtles, birds and squamates. This novel hypothesis describes the evolution of sex-determining mechanisms as turnover events accompanied by one or two small mutations.

Transcriptional evolution underlying vertebrate sexual development.

Explaining the diversity of vertebrate sex-determining mechanisms ranging from genotypic (GSD) to temperature-dependent (TSD) remains a developmental and evolutionary conundrum. Using a phylogenetic framework, we explore the transcriptional evolution during gonadogenesis of several genes involved in sexual development, combining novel data from Chrysemys picta turtles (TSD) and published data from other TSD and GSD vertebrates. Our novel C. picta dataset underscores Sf1 and Wt1 as potential activators of the thermosensitive period and uncovered the first evidence of Dax1 involvement in male development in a TSD vertebrate. Contrasting transcriptional profiles revealed male-biased Wt1 expression in fish while monomorphic expression is found in tetrapods but absent in turtles. Sf1 expression appears highly labile with transitions among testicular, ovarian, and non-sex-specific gonadal formation patterns among and within lineages. Dax1's dual role in ovarian and testicular formation is found in fish and mammals but is dosage-sensitive exclusively in eutherian mammals due to its X-linkage in this group. Contrastingly, Sox9 male-biased and Aromatase female-biased expression appear ancestral and virtually conserved throughout vertebrates despite significant heterochronic changes in expression as other elements likely replaced their function in early gonadogenesis. Finally, research avenues are highlighted to further study the evolution of the regulatory network of sexual development.

Reference-free discovery of nuclear SNPs permits accurate, sensitive identification of Carya (hickory) species and hybrids.

PREMISE: DNA-based species identification is critical when morphological identification is restricted, but DNA-based identification pipelines typically rely on the ability to compare homologous sequence data across species. Because many clades lack robust genomic resources, we present here a bioinformatics pipeline capable of generating genome-wide single-nucleotide polymorphism (SNP) data while circumventing the need for any reference genome or annotation data. METHODS: Using the SISRS bioinformatics pipeline, we generated de novo ortholog data for the genus Carya, isolating sites where genetic variation was restricted to a single Carya species (i.e., species-informative SNPs). We leveraged these SNPs to identify both full-species and hybrid Carya specimens, even at very low sequencing depths. RESULTS: We identified between 46,000 and 476,000 species-identifying SNPs for each of eight diploid Carya species, and all species identifications were concordant with the species of record. For all putative F1 hybrid specimens, both parental species were correctly identified in all cases, and more punctate patterns of introgression were detectable in more cryptic crosses. DISCUSSION: Bioinformatics pipelines that use only short-read sequencing data provide vital new tools enabling rapid expansion of DNA identification assays for model and non-model clades alike.

Haff disease associated with consumption of buffalofish (Ictiobus spp.) in the United States, 2010-2020, with confirmation of the causative species.

BACKGROUND: In the United States, buffalofish (Ictiobus spp.) are sporadically associated with sudden onset muscle pain and weakness due to rhabdomyolysis within 24 h of fish consumption (Haff disease). Previous genetic analyses of case-associated samples were unable to distinguish the three species of buffalofish that occur in the US, Ictiobus cyprinellus (bigmouth buffalo), Ictiobus bubalus (smallmouth buffalo), and Ictiobus niger (black buffalo). METHODS: Ten events were investigated between 2010 and 2020 and demographic and clinical information was collected for 24 individuals. Meal remnants were collected from 5 of 10 events with additional associated samples (n = 24) collected from another five of 10 events. Low-coverage whole-genome sequencing (genome skimming) was used to identify meal remnants. RESULTS: Patients (26-75 years of age) ranged from 1-4 per event, with 90% involving ≥2 individuals. Reported symptoms included muscle tenderness and weakness, nausea/vomiting, and brown/tea-colored urine. Median incubation period was 8 h. Ninety-six percent of cases were hospitalized with a median duration of four days. The most commonly reported laboratory finding was elevated creatine phosphokinase and liver transaminases. Treatment was supportive including intravenous fluids to prevent renal failure. Events occurred in California (1), Illinois (2), Louisiana (1), New York (1), Mississippi (1), Missouri (2), New Jersey (1), and Texas (1) with location of harvest, when known, being Illinois, Louisiana, Mississippi, Missouri, Texas, and Wisconsin. Meal remnants were identified as I. bubalus (n = 4) and I. niger (n = 1). Associated samples were identified as I. bubalus (n = 16), I. cyprinellus (n = 5), and I. niger (n = 3). DISCUSSION: Time course, presentation of illness, and clinical findings were all consistent with previous domestic cases of buffalofish-associated Haff disease. In contrast to previous reports that I. cyprinellus is the causative species in US cases, data indicate that all three buffalofish species are harvested but I. bubalus is most often associated with illness.

Screening the PRISM Library against Staphylococcus aureus Reveals a Sesquiterpene Lactone from Liriodendron tulipifera with Inhibitory Activity.

Infections caused by the bacterium Staphylococcus aureus continue to pose threats to human health and put a financial burden on the healthcare system. The overuse of antibiotics has contributed to mutations leading to the emergence of methicillin-resistant S. aureus, and there is a critical need for the discovery and development of new antibiotics to evade drug-resistant bacteria. Medicinal plants have shown promise as sources of new small-molecule therapeutics with potential uses against pathogenic infections. The principal Rhode Island secondary metabolite (PRISM) library is a botanical extract library generated from specimens in the URI Youngken Medicinal Garden by upper-division undergraduate students. PRISM extracts were screened for activity against strains of methicillin-susceptible S. aureus (MSSA). An extract generated from the tulip tree (Liriodendron tulipifera) demonstrated growth inhibition against MSSA, and a bioassay-guided approach identified a sesquiterpene lactone, laurenobiolide, as the active constituent. Intriguingly, its isomers, tulipinolide and epi-tulipinolide, lacked potent activity against MSSA. Laurenobiolide also proved to be more potent against MSSA than the structurally similar sesquiterpene lactones, costunolide and dehydrocostus lactone. Laurenobiolide was the most abundant in the twig bark of the tulip tree, supporting the twig bark's historical and cultural usage in poultices and teas.

Thermosensitive sex chromosome dosage compensation in ZZ/ZW softshell turtles, Apalone spinifera.

Sex chromosome dosage compensation (SCDC) overcomes gene-dose imbalances that disturb transcriptional networks, as when ZW females or XY males are hemizygous for Z/X genes. Mounting data from non-model organisms reveal diverse SCDC mechanisms, yet their evolution remains obscure, because most informative lineages with variable sex chromosomes are unstudied. Here, we discovered SCDC in turtles and an unprecedented thermosensitive SCDC in eukaryotes. We contrasted RNA-seq expression of Z-genes, their autosomal orthologues, and control autosomal genes in Apalone spinifera (ZZ/ZW) and Chrysemys picta turtles with temperature-dependent sex determination (TSD) (proxy for ancestral expression). This approach disentangled chromosomal context effects on Z-linked and autosomal expression, from lineage effects owing to selection or drift. Embryonic Apalone SCDC is tissue- and age-dependent, regulated gene-by-gene, complete in females via Z-upregulation in both sexes (Type IV) but partial and environmentally plastic via Z-downregulation in males (accentuated at colder temperature), present in female hatchlings and a weakly suggestive in adult liver (Type I). Results indicate that embryonic SCDC evolved with/after sex chromosomes in Apalone's family Tryonichidae, while co-opting Z-gene upregulation present in the TSD ancestor. Notably, Apalone's SCDC resembles pygmy snake's, and differs from the full-SCDC of Anolis lizards who share homologous sex chromosomes (XY), advancing our understanding of how XX/XY and ZZ/ZW systems compensate gene-dose imbalance. 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)'.

Utilizing Big Data to Identify Tiny Toxic Components: Digitalis.

The botanical genus Digitalis is equal parts colorful, toxic, and medicinal, and its bioactive compounds have a long history of therapeutic use. However, with an extremely narrow therapeutic range, even trace amounts of Digitalis can cause adverse effects. Using chemical methods, the United States Food and Drug Administration traced a 1997 case of Digitalis toxicity to a shipment of Plantago (a common ingredient in dietary supplements marketed to improve digestion) contaminated with Digitalis lanata. With increased accessibility to next generation sequencing technology, here we ask whether this case could have been cracked rapidly using shallow genome sequencing strategies (e.g., genome skims). Using a modified implementation of the Site Identification from Short Read Sequences (SISRS) bioinformatics pipeline with whole-genome sequence data, we generated over 2 M genus-level single nucleotide polymorphisms in addition to species-informative single nucleotide polymorphisms. We simulated dietary supplement contamination by spiking low quantities (0-10%) of Digitalis whole-genome sequence data into a background of commonly used ingredients in products marketed for "digestive cleansing" and reliably detected Digitalis at the genus level while also discriminating between Digitalis species. This work serves as a roadmap for the development of novel DNA-based assays to quickly and reliably detect the presence of toxic species such as Digitalis in food products or dietary supplements using genomic methods and highlights the power of harnessing the entire genome to identify botanical species.

Extensive genome-wide duplications in the eastern oyster (Crassostrea virginica).

Genomic structural variation is an important source of genetic and phenotypic diversity, playing a critical role in evolution. The recent availability of a high-quality reference genome for the eastern oyster, Crassostrea virginica, and whole-genome sequence data of samples from across the species range in the USA, provides an opportunity to explore structural variation across the genome of this species. Our analysis shows significantly greater individual-level duplications of regions across the genome than that of most model vertebrate species. Duplications are widespread across all ten chromosomes with variation in frequency per chromosome. The eastern oyster shows a large interindividual variation in duplications as well as particular chromosomal regions with a higher density of duplications. A high percentage of duplications seen in C. virginica lie completely within genes and exons, suggesting the potential for impacts on gene function. These results support the hypothesis that structural changes may play a significant role in standing genetic variation in C. virginica, and potentially have a role in their adaptive and evolutionary success. Altogether, these results suggest that copy number variation plays an important role in the genomic variation of C. virginica. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Extreme thermal fluctuations from climate change unexpectedly accelerate demographic collapse of vertebrates with temperature-dependent sex determination.

Global climate is warming rapidly, threatening vertebrates with temperature-dependent sex determination (TSD) by disrupting sex ratios and other traits. Less understood are the effects of increased thermal fluctuations predicted to accompany climate change. Greater fluctuations could accelerate feminization of species that produce females under warmer conditions (further endangering TSD animals), or counter it (reducing extinction risk). Here we use novel experiments exposing eggs of Painted Turtles (Chrysemys picta) to replicated profiles recorded in field nests plus mathematically-modified profiles of similar shape but wider oscillations, and develop a new mathematical model for analysis. We show that broadening fluctuations around naturally male-producing (cooler) profiles feminizes developing embryos, whereas embryos from warmer profiles remain female or die. This occurs presumably because wider oscillations around cooler profiles expose embryos to very low temperatures that inhibit development, and to feminizing temperatures where most embryogenesis accrues. Likewise, embryos incubated under broader fluctuations around warmer profiles experience mostly feminizing temperatures, some dangerously high (which increase mortality), and fewer colder values that are insufficient to induce male development. Therefore, as thermal fluctuations escalate with global warming, the feminization of TSD turtle populations could accelerate, facilitating extinction by demographic collapse. Aggressive global CO2 mitigation scenarios (RCP2.6) could prevent these risks, while intermediate actions (RCP4.5 and RCP6.0 scenarios) yield moderate feminization, highlighting the peril that insufficient reductions of greenhouse gas emissions pose for TSD taxa. If our findings are generalizable, TSD squamates, tuatara, and crocodilians that produce males at warmer temperatures could suffer accelerated masculinization, underscoring the broad taxonomic threats of climate change.

Thermal Response of Epigenetic Genes Informs Turtle Sex Determination with and without Sex Chromosomes.

Vertebrate sexual fate can be established by environmental cues (e.g., temperature-dependent sex determination, TSD) or by genetic content (genotypic sex determination, GSD). While methylation is implicated in TSD, the influence of broader epigenetic processes in sexual development remains obscure. Here, we investigated for the first time the embryonic gonadal expression of the genome-wide epigenetic machinery in turtles, including genes and noncoding RNAs (ncRNAs) involved in DNA/histone acetylation, methylation, ubiquitination, phosphorylation, and RNAi. This machinery was active and differentially thermosensitive in TSD versus GSD (ZZ/ZW) turtles. Methylation and histone acetylation genes responded the strongest. The results suggest these working hypotheses: (i) TSD might be mediated by epigenetically controlled hormonal pathways (via acetylation, methylation, and ncRNAs), or by (ii) hormonally controlled epigenetic processes, and (iii) key epigenetic events prior to the canonical thermosensitive period may explain differences between TSD and GSD. Novel epigenetic candidate regulators other than methylation were identified, including previously unknown ncRNAs that could potentially mediate gonadogenesis. These findings illuminate the molecular ecology of reptilian sex determination and permitted hypothesis building to help guide future functional studies on the epigenetic transduction of external cues in TSD versus GSD systems.

MeDIP-seq and nCpG analyses illuminate sexually dimorphic methylation of gonadal development genes with high historic methylation in turtle hatchlings with temperature-dependent sex determination.

BACKGROUND: DNA methylation alters gene expression but not DNA sequence and mediates some cases of phenotypic plasticity. Temperature-dependent sex determination (TSD) epitomizes phenotypic plasticity where environmental temperature drives embryonic sexual fate, as occurs commonly in turtles. Importantly, the temperature-specific transcription of two genes underlying gonadal differentiation is known to be induced by differential methylation in TSD fish, turtle and alligator. Yet, how extensive is the link between DNA methylation and TSD remains unclear. Here we test for broad differences in genome-wide DNA methylation between male and female hatchling gonads of the TSD painted turtle Chrysemys picta using methyl DNA immunoprecipitation sequencing, to identify differentially methylated candidates for future study. We also examine the genome-wide nCpG distribution (which affects DNA methylation) in painted turtles and test for historic methylation in genes regulating vertebrate gonadogenesis. RESULTS: Turtle global methylation was consistent with other vertebrates (57% of the genome, 78% of all CpG dinucleotides). Numerous genes predicted to regulate turtle gonadogenesis exhibited sex-specific methylation and were proximal to methylated repeats. nCpG distribution predicted actual turtle DNA methylation and was bimodal in gene promoters (as other vertebrates) and introns (unlike other vertebrates). Differentially methylated genes, including regulators of sexual development, had lower nCpG content indicative of higher historic methylation. CONCLUSIONS: Ours is the first evidence suggesting that sexually dimorphic DNA methylation is pervasive in turtle gonads (perhaps mediated by repeat methylation) and that it targets numerous regulators of gonadal development, consistent with the hypothesis that it may regulate thermosensitive transcription in TSD vertebrates. However, further research during embryogenesis will help test this hypothesis and the alternative that instead, most differential methylation observed in hatchlings is the by-product of sexual differentiation and not its cause.

Transcriptomic responses to environmental temperature by turtles with temperature-dependent and genotypic sex determination assessed by RNAseq inform the genetic architecture of embryonic gonadal development.

Vertebrate sexual fate is decided primarily by the individual's genotype (GSD), by the environmental temperature during development (TSD), or both. Turtles exhibit TSD and GSD, making them ideal to study the evolution of sex determination. Here we analyze temperature-specific gonadal transcriptomes (RNA-sequencing validated by qPCR) of painted turtles (Chrysemys picta TSD) before and during the thermosensitive period, and at equivalent stages in soft-shell turtles (Apalone spinifera-GSD), to test whether TSD's and GSD's transcriptional circuitry is identical but deployed differently between mechanisms. Our data show that most elements of the mammalian urogenital network are active during turtle gonadogenesis, but their transcription is generally more thermoresponsive in TSD than GSD, and concordant with their sex-specific function in mammals [e.g., upregulation of Amh, Ar, Esr1, Fog2, Gata4, Igf1r, Insr, and Lhx9 at male-producing temperature, and of ß-catenin, Foxl2, Aromatase (Cyp19a1), Fst, Nf-kb, Crabp2 at female-producing temperature in Chrysemys]. Notably, antagonistic elements in gonadogenesis (e.g., ß-catenin and Insr) were thermosensitive only in TSD early-embryos. Cirbp showed warm-temperature upregulation in both turtles disputing its purported key TSD role. Genes that may convert thermal inputs into sex-specific development (e.g., signaling and hormonal pathways, RNA-binding and heat-shock) were differentially regulated. Jak-Stat, Nf-κB, retinoic-acid, Wnt, and Mapk-signaling (not Akt and Ras-signaling) potentially mediate TSD thermosensitivity. Numerous species-specific ncRNAs (including Xist) were differentially-expressed, mostly upregulated at colder temperatures, as were unannotated loci that constitute novel TSD candidates. Cirbp showed warm-temperature upregulation in both turtles. Consistent transcription between turtles and alligator revealed putatively-critical reptilian TSD elements for male (Sf1, Amh, Amhr2) and female (Crabp2 and Hspb1) gonadogenesis. In conclusion, while preliminary, our data helps illuminate the regulation and evolution of vertebrate sex determination, and contribute genomic resources to guide further research into this fundamental biological process.

Physical Mapping and Refinement of the Painted Turtle Genome (Chrysemys picta) Inform Amniote Genome Evolution and Challenge Turtle-Bird Chromosomal Conservation.

Comparative genomics continues illuminating amniote genome evolution, but for many lineages our understanding remains incomplete. Here, we refine the assembly (CPI 3.0.3 NCBI AHGY00000000.2) and develop a cytogenetic map of the painted turtle (Chrysemys picta-CPI) genome, the first in turtles and in vertebrates with temperature-dependent sex determination. A comparison of turtle genomes with those of chicken, selected nonavian reptiles, and human revealed shared and novel genomic features, such as numerous chromosomal rearrangements. The largest conserved syntenic blocks between birds and turtles exist in four macrochromosomes, whereas rearrangements were evident in these and other chromosomes, disproving that turtles and birds retain fully conserved macrochromosomes for greater than 300 Myr. C-banding revealed large heterochromatic blocks in the centromeric region of only few chromosomes. The nucleolar-organizing region (NOR) mapped to a single CPI microchromosome, whereas in some turtles and lizards the NOR maps to nonhomologous sex-chromosomes, thus revealing independent translocations of the NOR in various reptilian lineages. There was no evidence for recent chromosomal fusions as interstitial telomeric-DNA was absent. Some repeat elements (CR1-like, Gypsy) were enriched in the centromeres of five chromosomes, whereas others were widespread in the CPI genome. Bacterial artificial chromosome (BAC) clones were hybridized to 18 of the 25 CPI chromosomes and anchored to a G-banded ideogram. Several CPI sex-determining genes mapped to five chromosomes, and homology was detected between yet other CPI autosomes and the globally nonhomologous sex chromosomes of chicken, other turtles, and squamates, underscoring the independent evolution of vertebrate sex-determining mechanisms.