Improved genome assembly of American alligator genome reveals conserved architecture of estrogen signaling

Evolution of ancient satellite DNAs in extant alligators and caimans (Crocodylia, Reptilia)

BACKGROUND

Crocodilians are one of the oldest extant vertebrate lineages, exhibiting a combination of evolutionary success and morphological resilience that has persisted throughout the history of life on Earth. This ability to endure over such a long geological time span is of great evolutionary importance. Here, we have utilized the combination of genomic and chromosomal data to identify and compare the full catalogs of satellite DNA families (satDNAs, i.e., the satellitomes) of 5 out of the 8 extant Alligatoridae species. As crocodilian genomes reveal ancestral patterns of evolution, by employing this multispecies data collection, we can investigate and assess how satDNA families evolve over time.

RESULTS

Alligators and caimans displayed a small number of satDNA families, ranging from 3 to 13 satDNAs in A. sinensis and C. latirostris, respectively. Together with little variation both within and between species it highlighted long-term conservation of satDNA elements throughout evolution. Furthermore, we traced the origin of the ancestral forms of all satDNAs belonging to the common ancestor of Caimaninae and Alligatorinae. Fluorescence in situ experiments showed distinct hybridization patterns for identical orthologous satDNAs, indicating their dynamic genomic placement.

CONCLUSIONS

Alligators and caimans possess one of the smallest satDNA libraries ever reported, comprising only four sets of satDNAs that are shared by all species. Besides, our findings indicated limited intraspecific variation in satellite DNA, suggesting that the majority of new satellite sequences likely evolved from pre-existing ones.

Convergent genomic signatures associated with vertebrate viviparity

BACKGROUND

Viviparity-live birth-is a complex and innovative mode of reproduction that has evolved repeatedly across the vertebrate Tree of Life. Viviparous species exhibit remarkable levels of reproductive diversity, both in the amount of care provided by the parent during gestation, and the ways in which that care is delivered. The genetic basis of viviparity has garnered increasing interest over recent years; however, such studies are often undertaken on small evolutionary timelines, and thus are not able to address changes occurring on a broader scale. Using whole genome data, we investigated the molecular basis of this innovation across the diversity of vertebrates to answer a long held question in evolutionary biology: is the evolution of convergent traits driven by convergent genomic changes?

RESULTS

We reveal convergent changes in protein family sizes, protein-coding regions, introns, and untranslated regions (UTRs) in a number of distantly related viviparous lineages. Specifically, we identify 15 protein families showing evidence of contraction or expansion associated with viviparity. We additionally identify elevated substitution rates in both coding and noncoding sequences in several viviparous lineages. However, we did not find any convergent changes-be it at the nucleotide or protein level-common to all viviparous lineages.

CONCLUSIONS

Our results highlight the value of macroevolutionary comparative genomics in determining the genomic basis of complex evolutionary transitions. While we identify a number of convergent genomic changes that may be associated with the evolution of viviparity in vertebrates, there does not appear to be a convergent molecular signature shared by all viviparous vertebrates. Ultimately, our findings indicate that a complex trait such as viviparity likely evolves with changes occurring in multiple different pathways.

Sex Determination and Ovarian Development in Reptiles and Amphibians: From Genetic Pathways to Environmental Influences

BACKGROUND

Reptiles and amphibians provide untapped potential for discovering how a diversity of genetic pathways and environmental conditions are incorporated into developmental processes that can lead to similar functional outcomes. These groups display a multitude of reproductive strategies, and whereas many attributes are conserved within groups and even across vertebrates, several aspects of sexual development show considerable variation.

SUMMARY

In this review, we focus our attention on the development of the reptilian and amphibian ovary. First, we review and describe the events leading to ovarian development, including sex determination and ovarian maturation, through a comparative lens. We then describe how these events are influenced by environmental factors, focusing on temperature and exposure to anthropogenic chemicals. Lastly, we identify critical knowledge gaps and future research directions that will be crucial to moving forward in our understanding of ovarian development and the influences of the environment in reptiles and amphibians.

KEY MESSAGES

Reptiles and amphibians provide excellent models for understanding the diversity of sex determination strategies and reproductive development. However, a greater understanding of the basic biology of these systems is necessary for deciphering the adaptive and potentially disruptive implications of embryo-by-environment interactions in a rapidly changing world.

Genome-wide DNA methylation patterns harbour signatures of hatchling sex and past incubation temperature in a species with environmental sex determination

Conservation of thermally sensitive species depends on monitoring organismal and population-level responses to environmental change in real time. Epigenetic processes are increasingly recognized as key integrators of environmental conditions into developmentally plastic responses, and attendant epigenomic data sets hold potential for revealing cryptic phenotypes relevant to conservation efforts. Here, we demonstrate the utility of genome-wide DNA methylation (DNAm) patterns in the face of climate change for a group of especially vulnerable species, those with temperature-dependent sex determination (TSD). Due to their reliance on thermal cues during development to determine sexual fate, contemporary shifts in temperature are predicted to skew offspring sex ratios and ultimately destabilize sensitive populations. Using reduced-representation bisulphite sequencing, we profiled the DNA methylome in blood cells of hatchling American alligators (Alligator mississippiensis), a TSD species lacking reliable markers of sexual dimorphism in early life stages. We identified 120 sex-associated differentially methylated cytosines (DMCs; FDR < 0.1) in hatchlings incubated under a range of temperatures, as well as 707 unique temperature-associated DMCs. We further developed DNAm-based models capable of predicting hatchling sex with 100% accuracy (in 20 training samples and four test samples) and past incubation temperature with a mean absolute error of 1.2°C (in four test samples) based on the methylation status of 20 and 24 loci, respectively. Though largely independent of epigenomic patterning occurring in the embryonic gonad during TSD, DNAm patterns in blood cells may serve as nonlethal markers of hatchling sex and past incubation conditions in conservation applications. These findings also raise intriguing questions regarding tissue-specific epigenomic patterning in the context of developmental plasticity.

Evolution of hormone-phenotype couplings and hormone-genome interactions

When selection favors a new relationship between a cue and a hormonally mediated response, adaptation can proceed by altering the hormonal signal that is produced or by altering the phenotypic response to the hormonal signal. The field of evolutionary endocrinology has made considerable progress toward understanding the evolution of hormonal signals, but we know much less about the evolution of hormone-phenotype couplings, particularly at the hormone-genome interface. We briefly review and classify the mechanisms through which these hormone-phenotype couplings likely evolve, using androgens and their receptors and genomic response elements to illustrate our view. We then present two empirical studies of hormone-phenotype couplings, one rooted in evolutionary quantitative genetics and another in comparative transcriptomics, each focused on the regulation of sexually dimorphic phenotypes by testosterone (T) in the brown anole lizard (Anolis sagrei). First, we illustrate the potential for hormone-phenotype couplings to evolve by showing that coloration of the dewlap (an ornament used in behavioral displays) exhibits significant heritability in its responsiveness to T, implying that anoles harbor genetic variance in the architecture of hormonal pleiotropy. Second, we combine T manipulations with analyses of the liver transcriptome to ask whether and how statistical methods for characterizing modules of co-expressed genes and in silico techniques for identifying androgen response elements (AREs) can improve our understanding of hormone-genome interactions. We conclude by emphasizing important avenues for future work at the hormone-genome interface, particularly those conducted in a comparative evolutionary framework.

The rattlesnake W chromosome: A GC-rich retroelement refugium with retained gene function across ancient evolutionary strata

Sex chromosomes diverge after the establishment of recombination suppression, resulting in differential sex-linkage of genes involved in genetic sex determination and dimorphic traits. This process produces systems of male or female heterogamety wherein the Y and W chromosomes are only present in one sex and are often highly degenerated. Sex-limited Y and W chromosomes contain valuable information about the evolutionary transition from autosomes to sex chromosomes, yet detailed characterizations of the structure, composition, and gene content of sex-limited chromosomes are lacking for many species. In this study, we characterize the female-specific W chromosome of the prairie rattlesnake (Crotalus viridis) and evaluate how recombination suppression and other processes have shaped sex chromosome evolution in ZW snakes. Our analyses indicate that the rattlesnake W chromosome is over 80% repetitive and that an abundance of GC-rich mdg4 elements has driven an overall high degree of GC-richness despite a lack of recombination. The W chromosome is also highly enriched for repeat sequences derived from endogenous retroviruses and likely acts as a “refugium” for these and other retroelements. We annotated 219 putatively functional W-linked genes across at least two evolutionary strata identified based on estimates of sequence divergence between Z and W gametologs. The youngest of these strata is relatively gene-rich, however gene expression across strata suggests retained gene function amidst a greater degree of degeneration following ancient recombination suppression. Functional annotation of W-linked genes indicates a specialization of the W chromosome for reproductive and developmental function since recombination suppression from the Z chromosome.

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 Reference Genome Assembly for the American Pika (Ochotona princeps)

The American pika (Ochotona princeps) is an alpine lagomorph found throughout western North America. Primarily inhabiting talus slopes at higher elevations (>2000 m), American pikas are well adapted to cold, montane environments. Warming climates on both historical and contemporary scales have contributed to population declines in American pikas, positioning them as a focal mammalian species for investigating the ecological effects of climate change. To support and expand ongoing research efforts, here, we present a highly contiguous and annotated reference genome assembly for the American pika (OchPri4.0). This assembly was produced using Dovetail de novo proximity ligation methods and annotated through the NCBI Eukaryotic Genome Annotation pipeline. The resulting assembly was chromosome- scale, with a total length of 2.23 Gb across 9350 scaffolds and a scaffold N50 of 75.8 Mb. The vast majority (>97%) of the total assembly length was found within 36 large scaffolds; 33 of these scaffolds correlated to whole autosomes, while the X chromosome was covered by 3 large scaffolds. Additionally, we identified 17 enriched gene ontology terms among American pika-specific genes putatively related to adaptation to high-elevation environments. This high-quality genome assembly will serve as a springboard for exploring the evolutionary underpinnings of behavioral, ecological, and taxonomic diversification in pikas as well as broader-scale eco-evolutionary questions pertaining to cold-adapted species in general.

Effects of global warming on species with temperature-dependent sex determination: Bridging the gap between empirical research and management

Global warming could threaten over 400 species with temperature-dependent sex determination (TSD) worldwide, including all species of sea turtle. During embryonic development, rising temperatures might lead to the overproduction of one sex and, in turn, could bias populations' sex ratios to an extent that threatens their persistence. If climate change predictions are correct, and biased sex ratios reduce population viability, species with TSD may go rapidly extinct unless adaptive mechanisms, whether behavioural, physiological or molecular, exist to buffer these temperature-driven effects. Here, we summarize the discovery of the TSD phenomenon and its still elusive evolutionary significance. We then review the molecular pathways underpinning TSD in model species, along with the hormonal mechanisms that interact with temperatures to determine an individual's sex. To illustrate evolutionary mechanisms that can affect sex determination, we focus on sea turtle biology, discussing both the adaptive potential of this threatened TSD taxon, and the risks associated with conservation mismanagement.

Genome reconstruction and haplotype phasing using chromosome conformation capture methodologies

Genomic analysis of individuals or organisms is predicated on the availability of high-quality reference and genotype information. With the rapidly dropping costs of high-throughput DNA sequencing, this is becoming readily available for diverse organisms and for increasingly large populations of individuals. Despite these advances, there are still aspects of genome sequencing that remain challenging for existing sequencing methods. This includes the generation of long-range contiguity during genome assembly, identification of structural variants in both germline and somatic tissues, the phasing of haplotypes in diploid organisms and the resolution of genome sequence for organisms derived from complex samples. These types of information are valuable for understanding the role of genome sequence and genetic variation on genome function, and numerous approaches have been developed to address them. Recently, chromosome conformation capture (3C) experiments, such as the Hi-C assay, have emerged as powerful tools to aid in these challenges for genome reconstruction. We will review the current use of Hi-C as a tool for aiding in genome sequencing, addressing the applications, strengths, limitations and potential future directions for the use of 3C data in genome analysis. We argue that unique features of Hi-C experiments make this data type a powerful tool to address challenges in genome sequencing, and that future integration of Hi-C data with alternative sequencing assays will facilitate the continuing revolution in genomic analysis and genome sequencing.

A new emu genome illuminates the evolution of genome configuration and nuclear architecture of avian chromosomes

Emu and other ratites are more informative than any other birds in reconstructing the evolution of the ancestral avian or vertebrate karyotype because of their much slower rate of genome evolution. Here, we generated a new chromosome-level genome assembly of a female emu, and estimated the tempo of chromosome evolution across major avian phylogenetic branches, by comparing it to chromosome-level genome assemblies of 11 other bird and one turtle species. We found ratites exhibited the lowest numbers of intra- and inter-chromosomal changes among birds since their divergence with turtles. The small-sized and gene-rich emu microchromosomes have frequent inter-chromosomal contacts that are associated with housekeeping genes, which appears to be driven by clustering their centromeres in the nuclear interior, away from the macrochromosomes in the nuclear periphery. Unlike nonratite birds, only less than one-third of the emu W Chromosome regions have lost homologous recombination and diverged between the sexes. The emu W is demarcated into a highly heterochromatic region (WS0) and another recently evolved region (WS1) with only moderate sequence divergence with the Z Chromosome. WS1 has expanded its inactive chromatin compartment, increased chromatin contacts within the region, and decreased contacts with the nearby regions, possibly influenced by the spreading of heterochromatin from WS0. These patterns suggest that alteration of chromatin conformation comprises an important early step of sex chromosome evolution. Overall, our results provide novel insights into the evolution of avian genome structure and sex chromosomes in three-dimensional space.

Accurate assembly of the olive baboon (Papio anubis) genome using long-read and Hi-C data

BACKGROUND

Baboons are a widely used nonhuman primate model for biomedical, evolutionary, and basic genetics research. Despite this importance, the genomic resources for baboons are limited. In particular, the current baboon reference genome Panu_3.0 is a highly fragmented, reference-guided (i.e., not fully de novo) assembly, and its poor quality inhibits our ability to conduct downstream genomic analyses.

FINDINGS

Here we present a de novo genome assembly of the olive baboon (Papio anubis) that uses data from several recently developed single-molecule technologies. Our assembly, Panubis1.0, has an N50 contig size of ∼1.46 Mb (as opposed to 139 kb for Panu_3.0) and has single scaffolds that span each of the 20 autosomes and the X chromosome.

CONCLUSIONS

We highlight multiple lines of evidence (including Bionano Genomics data, pedigree linkage information, and linkage disequilibrium data) suggesting that there are several large assembly errors in Panu_3.0, which have been corrected in Panubis1.0.

Sex determination, gonadal sex differentiation, and plasticity in vertebrate species

A diverse array of sex determination (SD) mechanisms, encompassing environmental to genetic, have been found to exist among vertebrates, covering a spectrum from fixed SD mechanisms (mammals) to functional sex change in fishes (sequential hermaphroditic fishes). A major landmark in vertebrate SD was the discovery of the SRY gene in 1990. Since that time, many attempts to clone an SRY ortholog from nonmammalian vertebrates remained unsuccessful, until 2002, when DMY/dmrt1by was discovered as the SD gene of a small fish, medaka. Surprisingly, however, DMY/dmrt1by was found in only 2 species among more than 20 species of medaka, suggesting a large diversity of SD genes among vertebrates. Considerable progress has been made over the last 3 decades, such that it is now possible to formulate reasonable paradigms of how SD and gonadal sex differentiation may work in some model vertebrate species. This review outlines our current understanding of vertebrate SD and gonadal sex differentiation, with a focus on the molecular and cellular mechanisms involved. An impressive number of genes and factors have been discovered that play important roles in testicular and ovarian differentiation. An antagonism between the male and female pathway genes exists in gonads during both sex differentiation and, surprisingly, even as adults, suggesting that, in addition to sex-changing fishes, gonochoristic vertebrates including mice maintain some degree of gonadal sexual plasticity into adulthood. Importantly, a review of various SD mechanisms among vertebrates suggests that this is the ideal biological event that can make us understand the evolutionary conundrums underlying speciation and species diversity.

Assessing the Ability of Developmentally Precocious Estrogen Signaling to Recapitulate Ovarian Transcriptomes and Follicle Dynamics in Alligators from a Contaminated Lake

BACKGROUND

Concern has grown in recent decades over anthropogenic contaminants that interfere with the functioning of endocrine hormones. However, mechanisms connecting developmental processes to pathologies associated with endocrine-disrupting chemical (EDC) exposure are poorly understood in naturally exposed populations.

OBJECTIVES

We sought to a) characterize divergence in ovarian transcriptomic and follicular profiles between alligators originating from a historically EDC-contaminated site, Lake Apopka, and a reference site; b) test the ability of developmentally precocious estrogen exposure to recapitulate site-associated patterns of divergence; and c) test whether treatment with exogenous follicle-stimulating hormone (FSH) is capable of rescuing phenotypes associated with contaminant exposure and/or embryonic estrogen treatment.

METHODS

Alligators eggs were collected from a contaminated site and a reference site, and a subset of eggs from the reference site were treated with estradiol (E2) during embryonic development prior to gonadal differentiation. After hatching, alligators were raised under controlled laboratory settings for 5 months. Juveniles from both sites were divided and treated with exogenous FSH. Histological analyses and RNA-sequencing were conducted to characterize divergence in ovarian follicle dynamics and transcriptomes between sites, between reference and E2-treated animals, and between FSH-treated and nontreated animals.

RESULTS

We observed broad site-of-origin divergence in ovarian transcriptomes and reductions in ovarian follicle density between juvenile alligators from Lake Apopka and the reference site. Treating embryos from the reference site with E2 overwhelmingly recapitulated transcriptional and histological alterations observed in Lake Apopka juveniles. Ovarian phenotypes observed in Lake Apopka alligators or resulting from estrogen treatment were only partially rescued by treatment with exogenous FSH.

DISCUSSION

Recapitulation of ovarian abnormalities by precocious E2 revealed a relatively simple mechanism underlying contaminant-induced pathologies in a historical example of environmental endocrine disruption. Findings reported here support a model where the developmental timing of estrogen signaling has the potential to permanently alter ovarian organization and function. https://doi.org/10.1289/EHP6627.

A High-Quality Reference Genome Assembly of the Saltwater Crocodile, Crocodylus porosus, Reveals Patterns of Selection in Crocodylidae

Crocodilians are an economically, culturally, and biologically important group. To improve researchers’ ability to study genome structure, evolution, and gene regulation in the clade, we generated a high-quality de novo genome assembly of the saltwater crocodile, Crocodylus porosus, from Illumina short read data from genomic libraries and in vitro proximity-ligation libraries. The assembled genome is 2,123.5 Mb, with N50 scaffold size of 17.7 Mb and N90 scaffold size of 3.8 Mb. We then annotated this new assembly, increasing the number of annotated genes by 74%. In total, 96% of 23,242 annotated genes were associated with a functional protein domain. Furthermore, multiple noncoding functional regions and mappable genetic markers were identified. Upon analysis and overlapping the results of branch length estimation and site selection tests for detecting potential selection, we found 16 putative genes under positive selection in crocodilians, 10 in C. porosus and 6 in Alligator mississippiensis. The annotated C. porosus genome will serve as an important platform for osmoregulatory, physiological, and sex determination studies, as well as an important reference in investigating the phylogenetic relationships of crocodilians, birds, and other tetrapods.

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.

Consequence of Paradigm Shift with Repeat Landscapes in Reptiles: Powerful Facilitators of Chromosomal Rearrangements for Diversity and Evolution

Reptiles are notable for the extensive genomic diversity and species richness among amniote classes, but there is nevertheless a need for detailed genome-scale studies. Although the monophyletic amniotes have recently been a focus of attention through an increasing number of genome sequencing projects, the abundant repetitive portion of the genome, termed the "repeatome", remains poorly understood across different lineages. Consisting predominantly of transposable elements or mobile and satellite sequences, these repeat elements are considered crucial in causing chromosomal rearrangements that lead to genomic diversity and evolution. Here, we propose major repeat landscapes in representative reptilian species, highlighting their evolutionary dynamics and role in mediating chromosomal rearrangements. Distinct karyotype variability, which is typically a conspicuous feature of reptile genomes, is discussed, with a particular focus on rearrangements correlated with evolutionary reorganization of micro- and macrochromosomes and sex chromosomes. The exceptional karyotype variation and extreme genomic diversity of reptiles are used to test several hypotheses concerning genomic structure, function, and evolution.

Identification and characterization of microRNAs (miRNAs) and their transposable element origins in the saltwater crocodile, Crocodylus porosus

MicroRNAs (miRNAs) are 18-24 nucleotide regulatory RNAs. They are involved in the regulation of genetic and biological pathways through post transcriptional gene silencing and/or translational repression. Data suggests a slow evolutionary rate for the saltwater crocodile (Crocodylus porosus) over the past several million years when compared to birds, the closest extant relatives of crocodilians. Understanding gene regulation in the saltwater crocodile in the context of relatively slow genomic change thus holds potential for the investigation of genomics, evolution, and adaptation. Utilizing eleven tissue types and sixteen small RNA libraries, we report 644 miRNAs in the saltwater crocodile with >78% of miRNAs being novel to crocodilians. We also identified potential targets for the miRNAs and analyzed the relationship of the miRNA repertoire to transposable elements (TEs). Results suggest an increased association of DNA transposons with miRNAs when compared to retrotransposons. This work reports the first comprehensive analysis of miRNAs in Crocodylus porosus and addresses the potential impacts of miRNAs in regulating the genome in the saltwater crocodile. In addition, the data suggests a supporting role of TEs as a source for miRNAs, adding to the increasing evidence that TEs play a significant role in the evolution of gene regulation.

Role for androgens in determination of ovarian fate in the common snapping turtle, Chelydra serpentina

Sex steroids are involved in sex determination in almost all vertebrates, including species with temperature-dependent sex determination (TSD). It is well established that aromatase and estrogens are involved in ovary determination in TSD species. In contrast, the role of non-aromatizable androgens in TSD is less clear. In this study, we used dihydrotestosterone (DHT) and an antagonist of the mammalian androgen receptor (flutamide) to examine the impact of androgens on sex determination in the snapping turtle. We incubated eggs at a male-producing temperature and treated embryos with drug delivery vehicle (5 L ethanol), DHT in vehicle, or flutamide in vehicle during the sex-determining period. We then measured expression of markers for ovarian and testicular development and genes involved in steroidogenesis. A subset of embryos and hatchlings were collected for histological analysis of gonad differentiation and sex determination. DHT and flutamide both induced ovarian development: 100% of vehicle-treated hatchlings had testes, while 60% of DHT-treated and 32% flutamide-treated hatchlings had ovaries. DHT and flutamide treatments also had feminizing effects on gene expression patterns and the structure of embryonic gonads. DHT treatment increased expression of FoxL2, androgen receptor, aromatase and several steroidogenic genes. Flutamide produced a similar, but weaker, pattern of gene expression. Genes involved in testis development (Sox9 and Amh) were influenced by flutamide treatment. Our findings support the hypothesis that androgens and the androgen receptor are involved in ovary determination in the common snapping turtle.

Mechanisms related to sexual determination by temperature in reptiles

A number of strategies have emerged that appear to relate to the evolution of mechanisms for sexual determination in vertebrates, among which are genetic sex determination caused by sex chromosomes and environmental sex determination, where environmental factors influence the phenotype of the sex of an individual. Within the reptile group, some orders such as: Chelonia, Crocodylia, Squamata and Rhynchocephalia, manifest one of the most intriguing and exciting environmental sexual determination mechanisms that exists, comprising temperature-dependent sex determination (TSD), where the temperature of incubation that the embryo experiences during its development is fundamental to establishing the sex of the individual. This makes them an excellent model for the study of sexual determination at the molecular, cellular and physiological level, as well as in terms of their implications at an evolutionary and ecological level. There are different hypotheses concerning how this process is triggered and this review aims to describe any new contributions to particular TSD hypotheses, analyzing them from the "eco-evo-devo" perspective.

Integrating natural history collections and comparative genomics to study the genetic architecture of convergent evolution

Evolutionary convergence has been long considered primary evidence of adaptation driven by natural selection and provides opportunities to explore evolutionary repeatability and predictability. In recent years, there has been increased interest in exploring the genetic mechanisms underlying convergent evolution, in part, owing to the advent of genomic techniques. However, the current 'genomics gold rush' in studies of convergence has overshadowed the reality that most trait classifications are quite broadly defined, resulting in incomplete or potentially biased interpretations of results. Genomic studies of convergence would be greatly improved by integrating deep 'vertical', natural history knowledge with 'horizontal' knowledge focusing on the breadth of taxonomic diversity. Natural history collections have and continue to be best positioned for increasing our comprehensive understanding of phenotypic diversity, with modern practices of digitization and databasing of morphological traits providing exciting improvements in our ability to evaluate the degree of morphological convergence. Combining more detailed phenotypic data with the well-established field of genomics will enable scientists to make progress on an important goal in biology: to understand the degree to which genetic or molecular convergence is associated with phenotypic convergence. Although the fields of comparative biology or comparative genomics alone can separately reveal important insights into convergent evolution, here we suggest that the synergistic and complementary roles of natural history collection-derived phenomic data and comparative genomics methods can be particularly powerful in together elucidating the genomic basis of convergent evolution among higher taxa. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Improving Illumina assemblies with Hi-C and long reads: An example with the North African dromedary

Researchers have assembled thousands of eukaryotic genomes using Illumina reads, but traditional mate‐pair libraries cannot span all repetitive elements, resulting in highly fragmented assemblies. However, both chromosome conformation capture techniques, such as Hi‐C and Dovetail Genomics Chicago libraries and long‐read sequencing, such as Pacific Biosciences and Oxford Nanopore, help span and resolve repetitive regions and therefore improve genome assemblies. One important livestock species of arid regions that does not have a high‐quality contiguous reference genome is the dromedary (Camelus dromedarius). Draft genomes exist but are highly fragmented, and a high‐quality reference genome is needed to understand adaptation to desert environments and artificial selection during domestication. Dromedaries are among the last livestock species to have been domesticated, and together with wild and domestic Bactrian camels, they are the only representatives of the Camelini tribe, which highlights their evolutionary significance. Here we describe our efforts to improve the North African dromedary genome. We used Chicago and Hi‐C sequencing libraries from Dovetail Genomics to resolve the order of previously assembled contigs, producing almost chromosome‐level scaffolds. Remaining gaps were filled with Pacific Biosciences long reads, and then scaffolds were comparatively mapped to chromosomes. Long reads added 99.32 Mbp to the total length of the new assembly. Dovetail Chicago and Hi‐C libraries increased the longest scaffold over 12‐fold, from 9.71 Mbp to 124.99 Mbp and the scaffold N50 over 50‐fold, from 1.48 Mbp to 75.02 Mbp. We demonstrate that Illumina de novo assemblies can be substantially upgraded by combining chromosome conformation capture and long‐read sequencing.

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.

New Approaches for Genome Assembly and Scaffolding

Affordable, high-throughput DNA sequencing has accelerated the pace of genome assembly over the past decade. Genome assemblies from high-throughput, short-read sequencing, however, are often not as contiguous as the first generation of genome assemblies. Whereas early genome assembly projects were often aided by clone maps or other mapping data, many current assembly projects forego these scaffolding data and only assemble genomes into smaller segments. Recently, new technologies have been invented that allow chromosome-scale assembly at a lower cost and faster speed than traditional methods. Here, we give an overview of the problem of chromosome-scale assembly and traditional methods for tackling this problem. We then review new technologies for chromosome-scale assembly and recent genome projects that used these technologies to create highly contiguous genome assemblies at low cost.

Divergent evolution in the genomes of closely related lacertids, Lacerta viridis and L. bilineata, and implications for speciation

BACKGROUND

Lacerta viridis and Lacerta bilineata are sister species of European green lizards (eastern and western clades, respectively) that, until recently, were grouped together as the L. viridis complex. Genetic incompatibilities were observed between lacertid populations through crossing experiments, which led to the delineation of two separate species within the L. viridis complex. The population history of these sister species and processes driving divergence are unknown. We constructed the first high-quality de novo genome assemblies for both L. viridis and L. bilineata through Illumina and PacBio sequencing, with annotation support provided from transcriptome sequencing of several tissues. To estimate gene flow between the two species and identify factors involved in reproductive isolation, we studied their evolutionary history, identified genomic rearrangements, detected signatures of selection on non-coding RNA, and on protein-coding genes.

FINDINGS

Here we show that gene flow was primarily unidirectional from L. bilineata to L. viridis after their split at least 1.15 million years ago. We detected positive selection of the non-coding repertoire; mutations in transcription factors; accumulation of divergence through inversions; selection on genes involved in neural development, reproduction, and behavior, as well as in ultraviolet-response, possibly driven by sexual selection, whose contribution to reproductive isolation between these lacertid species needs to be further evaluated.

CONCLUSION

The combination of short and long sequence reads resulted in one of the most complete lizard genome assemblies. The characterization of a diverse array of genomic features provided valuable insights into the demographic history of divergence among European green lizards, as well as key species differences, some of which are candidates that could have played a role in speciation. In addition, our study generated valuable genomic resources that can be used to address conservation-related issues in lacertids.

Genome assembly with in vitro proximity ligation data and whole-genome triplication in lettuce

Lettuce (Lactuca sativa) is a major crop and a member of the large, highly successful Compositae family of flowering plants. Here we present a reference assembly for the species and family. This was generated using whole-genome shotgun Illumina reads plus in vitro proximity ligation data to create large superscaffolds; it was validated genetically and superscaffolds were oriented in genetic bins ordered along nine chromosomal pseudomolecules. We identify several genomic features that may have contributed to the success of the family, including genes encoding Cycloidea-like transcription factors, kinases, enzymes involved in rubber biosynthesis and disease resistance proteins that are expanded in the genome. We characterize 21 novel microRNAs, one of which may trigger phasiRNAs from numerous kinase transcripts. We provide evidence for a whole-genome triplication event specific but basal to the Compositae. We detect 26% of the genome in triplicated regions containing 30% of all genes that are enriched for regulatory sequences and depleted for genes involved in defence.