The origins and evolution of chromosomes, dosage compensation, and mechanisms underlying venom regulation in snakes.

Here we use a chromosome-level genome assembly of a prairie rattlesnake (Crotalus viridis), together with Hi-C, RNA-seq, and whole-genome resequencing data, to study key features of genome biology and evolution in reptiles. We identify the rattlesnake Z Chromosome, including the recombining pseudoautosomal region, and find evidence for partial dosage compensation driven by an evolutionary accumulation of a female-biased up-regulation mechanism. Comparative analyses with other amniotes provide new insight into the origins, structure, and function of reptile microchromosomes, which we demonstrate have markedly different structure and function compared to macrochromosomes. Snake microchromosomes are also enriched for venom genes, which we show have evolved through multiple tandem duplication events in multiple gene families. By overlaying chromatin structure information and gene expression data, we find evidence for venom gene-specific chromatin contact domains and identify how chromatin structure guides precise expression of multiple venom gene families. Further, we find evidence for venom gland-specific transcription factor activity and characterize a complement of mechanisms underlying venom production and regulation. Our findings reveal novel and fundamental features of reptile genome biology, provide insight into the regulation of snake venom, and broadly highlight the biological insight enabled by chromosome-level genome assemblies.

Transcriptomic Analysis Reveals Potential Candidate Pathways and Genes Involved in Toxin Biosynthesis in True Toads.

Synthesized chemical defenses have broadly evolved across countless taxa and are important in shaping evolutionary and ecological interactions within ecosystems. However, the underlying genomic mechanisms by which these organisms synthesize and utilize their toxins are relatively unknown. Herein, we use comparative transcriptomics to uncover potential toxin synthesizing genes and pathways, as well as interspecific patterns of toxin synthesizing genes across 10 species of North American true toads (Bufonidae). Upon assembly and annotation of the 10 transcriptomes, we explored patterns of relative gene expression and possible protein-protein interactions across the species to determine what genes and/or pathways may be responsible for toxin synthesis. We also tested our transcriptome dataset for signatures of positive selection to reveal how selection may be acting upon potential toxin producing genes. We assembled high-quality transcriptomes of the bufonid parotoid gland, a tissue not often investigated in other bufonid-related RNAseq studies. We found several genes involved in metabolic and biosynthetic pathways (e.g., steroid biosynthesis, terpenoid backbone biosynthesis, isoquinoline biosynthesis, and glucosinolate biosynthesis) that were functionally enriched and/or relatively expressed across the 10 focal species that may be involved in the synthesis of alkaloid and steroid toxins, as well as other small metabolic compounds that cause distastefulness in bufonids. We hope that our study lays a foundation for future studies to explore the genomic underpinnings and specific pathways of toxin synthesis in toads, as well as at the macroevolutionary scale across numerous taxa that produce their own defensive toxins.

Divergence time estimation of genus Tribolium by extensive sampling of highly conserved orthologs.

Tribolium castaneum, the red flour beetle, is among the most well-studied eukaryotic genetic model organisms. Tribolium often serves as a comparative bridge from highly derived Drosophila traits to other organisms. Simultaneously, as a member of the most diverse order of metazoans, Coleoptera, Tribolium informs us about innovations that accompany hyper diversity. However, understanding the tempo and mode of evolutionary innovation requires well-resolved, time-calibrated phylogenies, which are not available for Tribolium. The most recent effort to understand Tribolium phylogenetics used two mitochondrial and three nuclear markers. The study concluded that the genus may be paraphyletic and reported a broad range for divergence time estimates. Here we employ recent advances in Bayesian methods to estimate the relationships and divergence times among Tribolium castaneum, T. brevicornis, T. confusum, T. freemani, and Gnatocerus cornutus using 1368 orthologs conserved across all five species and an independent substitution rate estimate. We find that the most basal split within Tribolium occurred ~86 Mya [95% HPD 85.90-87.04 Mya] and that the most recent split was between T. freemani and T. castaneum at ~14 Mya [95% HPD 13.55-14.00]. Our results are consistent with broader phylogenetic analyses of insects and suggest that Cenozoic climate changes played a role in the Tribolium diversification.

scRNA-seq Reveals Novel Genetic Pathways and Sex Chromosome Regulation in Tribolium Spermatogenesis.

Spermatogenesis is critical to sexual reproduction yet evolves rapidly in many organisms. High-throughput single-cell transcriptomics promises unparalleled insight into this important process but understanding can be impeded in nonmodel systems by a lack of known genes that can reliably demarcate biologically meaningful cell populations. Tribolium castaneum, the red flour beetle, lacks known markers for spermatogenesis found in insect species like Drosophila melanogaster. Using single-cell sequencing data collected from adult beetle testes, we implement a strategy for elucidating biologically meaningful cell populations by using transient expression stage identification markers, weighted principal component clustering, and SNP-based haploid/diploid phasing. We identify populations that correspond to observable points in sperm differentiation and find species specific markers for each stage. Our results indicate that molecular pathways underlying spermatogenesis in Coleoptera are substantially diverged from those in Diptera. We also show that most genes on the X chromosome experience meiotic sex chromosome inactivation. Temporal expression of Drosophila MSL complex homologs coupled with spatial analysis of potential chromatin entry sites further suggests that the dosage compensation machinery may mediate escape from meiotic sex chromosome inactivation and postmeiotic reactivation of the X chromosome.

Improvements to the Gulf pipefish Syngnathus scovelli genome.

The Gulf pipefish Syngnathus scovelli has emerged as an important species for studying sexual selection, development, and physiology. Comparative evolutionary genomics research involving fishes from Syngnathidae depends on having a high-quality genome assembly and annotation. However, the first S. scovelli genome assembled using short-read sequences and a smaller RNA-sequence dataset has limited contiguity and a relatively poor annotation. Here, using PacBio long-read high-fidelity sequences and a proximity ligation library, we generate an improved assembly to obtain 22 chromosome-level scaffolds. Compared to the first assembly, the gaps in the improved assembly are smaller, the N75 is larger, and our genome is ~95% BUSCO complete. Using a large body of RNA-Seq reads from different tissue types and NCBI's Eukaryotic Annotation Pipeline, we discovered 28,162 genes, of which 8,061 are non-coding genes. Our new genome assembly and annotation are tagged as a RefSeq genome by NCBI and provide enhanced resources for research work involving S. scovelli..