Extant and extinct bilby genomes combined with Indigenous knowledge improve conservation of a unique Australian marsupial.

Ninu (greater bilby, Macrotis lagotis) are desert-dwelling, culturally and ecologically important marsupials. In collaboration with Indigenous rangers and conservation managers, we generated the Ninu chromosome-level genome assembly (3.66 Gbp) and genome sequences for the extinct Yallara (lesser bilby, Macrotis leucura). We developed and tested a scat single-nucleotide polymorphism panel to inform current and future conservation actions, undertake ecological assessments and improve our understanding of Ninu genetic diversity in managed and wild populations. We also assessed the beneficial impact of translocations in the metapopulation (N = 363 Ninu). Resequenced genomes (temperate Ninu, 6; semi-arid Ninu, 6; and Yallara, 4) revealed two major population crashes during global cooling events for both species and differences in Ninu genes involved in anatomical and metabolic pathways. Despite their 45-year captive history, Ninu have fewer long runs of homozygosity than other larger mammals, which may be attributable to their boom-bust life history. Here we investigated the unique Ninu biology using 12 tissue transcriptomes revealing expression of all 115 conserved eutherian chorioallantoic placentation genes in the uterus, an XY1Y2 sex chromosome system and olfactory receptor gene expansions. Together, we demonstrate the holistic value of genomics in improving key conservation actions, understanding unique biological traits and developing tools for Indigenous rangers to monitor remote wild populations.

Gene expression of male pathway genes sox9 and amh during early sex differentiation in a reptile departs from the classical amniote model.

BACKGROUND: Sex determination is the process whereby the bipotential embryonic gonads become committed to differentiate into testes or ovaries. In genetic sex determination (GSD), the sex determining trigger is encoded by a gene on the sex chromosomes, which activates a network of downstream genes; in mammals these include SOX9, AMH and DMRT1 in the male pathway, and FOXL2 in the female pathway. Although mammalian and avian GSD systems have been well studied, few data are available for reptilian GSD systems. RESULTS: We conducted an unbiased transcriptome-wide analysis of gonad development throughout differentiation in central bearded dragon (Pogona vitticeps) embryos with GSD. We found that sex differentiation of transcriptomic profiles occurs at a very early stage, before the gonad consolidates as a body distinct from the gonad-kidney complex. The male pathway genes dmrt1 and amh and the female pathway gene foxl2 play a key role in early sex differentiation in P. vitticeps, but the central player of the mammalian male trajectory, sox9, is not differentially expressed in P. vitticeps at the bipotential stage. The most striking difference from GSD systems of other amniotes is the high expression of the male pathway genes amh and sox9 in female gonads during development. We propose that a default male trajectory progresses if not repressed by a W-linked dominant gene that tips the balance of gene expression towards the female trajectory. Further, weighted gene expression correlation network analysis revealed novel candidates for male and female sex differentiation. CONCLUSION: Our data reveal that interpretation of putative mechanisms of GSD in reptiles cannot solely depend on lessons drawn from mammals.

Polygenic sex determination in vertebrates - is there any such thing?

Genetic sex determination (SD) in most vertebrates is controlled by a single master sex gene, which ensures a 1:1 sex ratio. However, more complex systems abound, and several have been ascribed to polygenic SD (PSD), in which many genes at different loci interact to produce the sexual phenotype. Here we examine claims for PSD in vertebrates, finding that most constitute transient states during sex chromosome turnover, or aberrant systems in species hybrids. To avoid confusion about terminology, we propose a consistent nomenclature for genetic SD systems.

Truncated jarid2 and kdm6b transcripts are associated with temperature-induced sex reversal during development in a dragon lizard.

Sex determination and differentiation in reptiles are complex. In the model species, Pogona vitticeps, high incubation temperature can cause male to female sex reversal. To elucidate the epigenetic mechanisms of thermolabile sex, we used an unbiased genome-wide assessment of intron retention during sex reversal. The previously implicated chromatin modifiers (jarid2 and kdm6b) were two of three genes to display sex reversal-specific intron retention. In these species, embryonic intron retention resulting in C-terminally truncated jarid2 and kdm6b isoforms consistently occurs at low temperatures. High-temperature sex reversal is uniquely characterized by a high prevalence of N-terminally truncated isoforms of jarid2 and kdm6b, which are not present at low temperatures, or in two other reptiles with temperature-dependent sex determination. This work verifies that chromatin-modifying genes are involved in highly conserved temperature responses and can also be transcribed into isoforms with new sex-determining roles.

Two transcriptionally distinct pathways drive female development in a reptile with both genetic and temperature dependent sex determination.

How temperature determines sex remains unknown. A recent hypothesis proposes that conserved cellular mechanisms (calcium and redox; 'CaRe' status) sense temperature and identify genes and regulatory pathways likely to be involved in driving sexual development. We take advantage of the unique sex determining system of the model organism, Pogona vitticeps, to assess predictions of this hypothesis. P. vitticeps has ZZ male: ZW female sex chromosomes whose influence can be overridden in genetic males by high temperatures, causing male-to-female sex reversal. We compare a developmental transcriptome series of ZWf females and temperature sex reversed ZZf females. We demonstrate that early developmental cascades differ dramatically between genetically driven and thermally driven females, later converging to produce a common outcome (ovaries). We show that genes proposed as regulators of thermosensitive sex determination play a role in temperature sex reversal. Our study greatly advances the search for the mechanisms by which temperature determines sex.

Chromosomics: Bridging the Gap between Genomes and Chromosomes.

The recent advances in DNA sequencing technology are enabling a rapid increase in the number of genomes being sequenced. However, many fundamental questions in genome biology remain unanswered, because sequence data alone is unable to provide insight into how the genome is organised into chromosomes, the position and interaction of those chromosomes in the cell, and how chromosomes and their interactions with each other change in response to environmental stimuli or over time. The intimate relationship between DNA sequence and chromosome structure and function highlights the need to integrate genomic and cytogenetic data to more comprehensively understand the role genome architecture plays in genome plasticity. We propose adoption of the term 'chromosomics' as an approach encompassing genome sequencing, cytogenetics and cell biology, and present examples of where chromosomics has already led to novel discoveries, such as the sex-determining gene in eutherian mammals. More importantly, we look to the future and the questions that could be answered as we enter into the chromosomics revolution, such as the role of chromosome rearrangements in speciation and the role more rapidly evolving regions of the genome, like centromeres, play in genome plasticity. However, for chromosomics to reach its full potential, we need to address several challenges, particularly the training of a new generation of cytogeneticists, and the commitment to a closer union among the research areas of genomics, cytogenetics, cell biology and bioinformatics. Overcoming these challenges will lead to ground-breaking discoveries in understanding genome evolution and function.

Stress, novel sex genes, and epigenetic reprogramming orchestrate socially controlled sex change.

Bluehead wrasses undergo dramatic, socially cued female-to-male sex change. We apply transcriptomic and methylome approaches in this wild coral reef fish to identify the primary trigger and subsequent molecular cascade of gonadal metamorphosis. Our data suggest that the environmental stimulus is exerted via the stress axis and that repression of the aromatase gene (encoding the enzyme converting androgens to estrogens) triggers a cascaded collapse of feminizing gene expression and identifies notable sex-specific gene neofunctionalization. Furthermore, sex change involves distinct epigenetic reprogramming and an intermediate state with altered epigenetic machinery expression akin to the early developmental cells of mammals. These findings reveal at a molecular level how a normally committed developmental process remains plastic and is reversed to completely alter organ structures.

The Methylome of Vertebrate Sex Chromosomes.

DNA methylation is a key epigenetic modification in vertebrate genomes known to be involved in the regulation of gene expression, X chromosome inactivation, genomic imprinting, chromatin structure, and control of transposable elements. DNA methylation is common to all eukaryote genomes, but we still lack a complete understanding of the variation in DNA methylation patterns on sex chromosomes and between the sexes in diverse species. To better understand sex chromosome DNA methylation patterns between different amniote vertebrates, we review literature that has analyzed the genome-wide distribution of DNA methylation in mammals and birds. In each system, we focus on DNA methylation patterns on the autosomes versus the sex chromosomes.

Landscape of DNA Methylation on the Marsupial X.

DNA methylation plays a key role in maintaining transcriptional silence on the inactive X chromosome of eutherian mammals. Beyond eutherians, there are limited genome wide data on DNA methylation from other vertebrates. Previous studies of X borne genes in various marsupial models revealed no differential DNA methylation of promoters between the sexes, leading to the conclusion that CpG methylation plays no role in marsupial X-inactivation. Using reduced representation bisulfite sequencing, we generated male and female CpG methylation profiles in four representative vertebrates (mouse, gray short-tailed opossum, platypus, and chicken). A variety of DNA methylation patterns were observed. Platypus and chicken displayed no large-scale differential DNA methylation between the sexes on the autosomes or the sex chromosomes. As expected, a metagene analysis revealed hypermethylation at transcription start sites (TSS) of genes subject to X-inactivation in female mice. This contrasted with the opossum, in which metagene analysis did not detect differential DNA methylation between the sexes at TSSs of genes subject to X-inactivation. However, regions flanking TSSs of these genes were hypomethylated. Our data are the first to demonstrate that, for genes subject to X-inactivation in both eutherian and marsupial mammals, there is a consistent difference between DNA methylation levels at TSSs and immediate flanking regions, which we propose has a silencing effect in both groups.

Differential intron retention in Jumonji chromatin modifier genes is implicated in reptile temperature-dependent sex determination.

In many vertebrates, sex of offspring is determined by external environmental cues rather than by sex chromosomes. In reptiles, for instance, temperature-dependent sex determination (TSD) is common. Despite decades of work, the mechanism by which temperature is converted into a sex-determining signal remains mysterious. This is partly because it is difficult to distinguish the primary molecular events of TSD from the confounding downstream signatures of sexual differentiation. We use the Australian central bearded dragon, in which chromosomal sex determination is overridden at high temperatures to produce sex-reversed female offspring, as a unique model to identify TSD-specific features of the transcriptome. We show that an intron is retained in mature transcripts from each of two Jumonji family genes, JARID2 and JMJD3, in female dragons that have been sex-reversed by temperature but not in normal chromosomal females or males. JARID2 is a component of the master chromatin modifier Polycomb Repressive Complex 2, and the mammalian sex-determining factor SRY is directly regulated by an independent but closely related Jumonji family member. We propose that the perturbation of JARID2/JMJD3 function by intron retention alters the epigenetic landscape to override chromosomal sex-determining cues, triggering sex reversal at extreme temperatures. Sex reversal may then facilitate a transition from genetic sex determination to TSD, with JARID2/JMJD3 intron retention preserved as the decisive regulatory signal. Significantly, we also observe sex-associated differential retention of the equivalent introns in JARID2/JMJD3 transcripts expressed in embryonic gonads from TSD alligators and turtles, indicative of a reptile-wide mechanism controlling TSD.

Sex reversal triggers the rapid transition from genetic to temperature-dependent sex.

Sex determination in animals is amazingly plastic. Vertebrates display contrasting strategies ranging from complete genetic control of sex (genotypic sex determination) to environmentally determined sex (for example, temperature-dependent sex determination). Phylogenetic analyses suggest frequent evolutionary transitions between genotypic and temperature-dependent sex determination in environmentally sensitive lineages, including reptiles. These transitions are thought to involve a genotypic system becoming sensitive to temperature, with sex determined by gene-environment interactions. Most mechanistic models of transitions invoke a role for sex reversal. Sex reversal has not yet been demonstrated in nature for any amniote, although it occurs in fish and rarely in amphibians. Here we make the first report of reptile sex reversal in the wild, in the Australian bearded dragon (Pogona vitticeps), and use sex-reversed animals to experimentally induce a rapid transition from genotypic to temperature-dependent sex determination. Controlled mating of normal males to sex-reversed females produces viable and fertile offspring whose phenotypic sex is determined solely by temperature (temperature-dependent sex determination). The W sex chromosome is eliminated from this lineage in the first generation. The instantaneous creation of a lineage of ZZ temperature-sensitive animals reveals a novel, climate-induced pathway for the rapid transition between genetic and temperature-dependent sex determination, and adds to concern about adaptation to rapid global climate change.

Major Histocompatibility Complex Genes Map to Two Chromosomes in an Evolutionarily Ancient Reptile, the Tuatara Sphenodon punctatus.

Major histocompatibility complex (MHC) genes are a central component of the vertebrate immune system and usually exist in a single genomic region. However, considerable differences in MHC organization and size exist between different vertebrate lineages. Reptiles occupy a key evolutionary position for understanding how variation in MHC structure evolved in vertebrates, but information on the structure of the MHC region in reptiles is limited. In this study, we investigate the organization and cytogenetic location of MHC genes in the tuatara (Sphenodon punctatus), the sole extant representative of the early-diverging reptilian order Rhynchocephalia. Sequencing and mapping of 12 clones containing class I and II MHC genes from a bacterial artificial chromosome library indicated that the core MHC region is located on chromosome 13q. However, duplication and translocation of MHC genes outside of the core region was evident, because additional class I MHC genes were located on chromosome 4p. We found a total of seven class I sequences and 11 class II ß sequences, with evidence for duplication and pseudogenization of genes within the tuatara lineage. The tuatara MHC is characterized by high repeat content and low gene density compared with other species and we found no antigen processing or MHC framework genes on the MHC gene-containing clones. Our findings indicate substantial differences in MHC organization in tuatara compared with mammalian and avian MHCs and highlight the dynamic nature of the MHC. Further sequencing and annotation of tuatara and other reptile MHCs will determine if the tuatara MHC is representative of nonavian reptiles in general.

Pathogenesis and molecular biology of a transmissible tumor in the Tasmanian devil.

The emergence of a fatal transmissible cancer known as devil facial tumor disease (DFTD) is threatening the iconic Tasmanian devil with extinction in the wild within the next few decades. Since the first report of the disease in 1996, DFTD has spread to over 85% of the devils' distribution and dramatically reduced devil numbers. Research into DFTD has focused on gaining a deeper understanding of the disease on multiple levels, including an accurate assessment of the tissue origin of the tumor, elucidation of how the tumor evades immune detection, and determination of how the tumor is transmitted between individuals and how it is evolving as it spreads through the population. Knowledge gained from these studies has important implications for DFTD management and devil conservation.

Highly differentiated ZW sex microchromosomes in the Australian Varanus species evolved through rapid amplification of repetitive sequences.

Transitions between sex determination systems have occurred in many lineages of squamates and it follows that novel sex chromosomes will also have arisen multiple times. The formation of sex chromosomes may be reinforced by inhibition of recombination and the accumulation of repetitive DNA sequences. The karyotypes of monitor lizards are known to be highly conserved yet the sex chromosomes in this family have not been fully investigated. Here, we compare male and female karyotypes of three Australian monitor lizards, Varanus acanthurus, V. gouldii and V. rosenbergi, from two different clades. V. acanthurus belongs to the acanthurus clade and the other two belong to the gouldii clade. We applied C-banding and comparative genomic hybridization to reveal that these species have ZZ/ZW sex micro-chromosomes in which the W chromosome is highly differentiated from the Z chromosome. In combination with previous reports, all six Varanus species in which sex chromosomes have been identified have ZZ/ZW sex chromosomes, spanning several clades on the varanid phylogeny, making it likely that the ZZ/ZW sex chromosome is ancestral for this family. However, repetitive sequences of these ZW chromosome pairs differed among species. In particular, an (AAT)n microsatellite repeat motif mapped by fluorescence in situ hybridization on part of W chromosome in V. acanthurus only, whereas a (CGG)n motif mapped onto the W chromosomes of V. gouldii and V. rosenbergi. Furthermore, the W chromosome probe for V. acanthurus produced hybridization signals only on the centromeric regions of W chromosomes of the other two species. These results suggest that the W chromosome sequences were not conserved between gouldii and acanthurus clades and that these repetitive sequences have been amplified rapidly and independently on the W chromosome of the two clades after their divergence.

Sequence and gene content of a large fragment of a lizard sex chromosome and evaluation of candidate sex differentiating gene R-spondin 1.

BACKGROUND: Scant genomic information from non-avian reptile sex chromosomes is available, and for only a few lizards, several snakes and one turtle species, and it represents only a small fraction of the total sex chromosome sequences in these species. RESULTS: We report a 352 kb of contiguous sequence from the sex chromosome of a squamate reptile, Pogona vitticeps, with a ZZ/ZW sex microchromosome system. This contig contains five protein coding genes (oprd1, rcc1, znf91, znf131, znf180), and major families of repetitive sequences with a high number of copies of LTR and non-LTR retrotransposons, including the CR1 and Bov-B LINEs. The two genes, oprd1 and rcc1 are part of a homologous syntenic block, which is conserved among amniotes. While oprd1 and rcc1 have no known function in sex determination or differentiation in amniotes, this homologous syntenic block in mammals and chicken also contains R-spondin 1 (rspo1), the ovarian differentiating gene in mammals. In order to explore the probability that rspo1 is sex determining in dragon lizards, genomic BAC and cDNA clones were mapped using fluorescence in situ hybridisation. Their location on an autosomal microchromosome pair, not on the ZW sex microchromosomes, eliminates rspo1 as a candidate sex determining gene in P. vitticeps. CONCLUSION: Our study has characterized the largest contiguous stretch of physically mapped sex chromosome sequence (352 kb) from a ZZ/ZW lizard species. Although this region represents only a small fraction of the sex chromosomes of P. vitticeps, it has revealed several features typically associated with sex chromosomes including the accumulation of large blocks of repetitive sequences.

Independent evolution of transcriptional inactivation on sex chromosomes in birds and mammals.

X chromosome inactivation in eutherian mammals has been thought to be tightly controlled, as expected from a mechanism that compensates for the different dosage of X-borne genes in XX females and XY males. However, many X genes escape inactivation in humans, inactivation of the X in marsupials is partial, and the unrelated sex chromosomes of monotreme mammals have incomplete and gene-specific inactivation of X-linked genes. The bird ZW sex chromosome system represents a third independently evolved amniote sex chromosome system with dosage compensation, albeit partial and gene-specific, via an unknown mechanism (i.e. upregulation of the single Z in females, down regulation of one or both Zs in males, or a combination). We used RNA-fluorescent in situ hybridization (RNA-FISH) to demonstrate, on individual fibroblast cells, inactivation of 11 genes on the chicken Z and 28 genes on the X chromosomes of platypus. Each gene displayed a reproducible frequency of 1Z/1X-active and 2Z/2X-active cells in the homogametic sex. Our results indicate that the probability of inactivation is controlled on a gene-by-gene basis (or small domains) on the chicken Z and platypus X chromosomes. This regulatory mechanism must have been exapted independently to the non-homologous sex chromosomes in birds and mammals in response to an over-expressed Z or X in the homogametic sex, highlighting the universal importance that (at least partial) silencing plays in the evolution on amniote dosage compensation and, therefore, the differentiation of sex chromosomes.

A first-generation integrated tammar wallaby map and its use in creating a tammar wallaby first-generation virtual genome map.

BACKGROUND: The limited (2X) coverage of the tammar wallaby (Macropus eugenii) genome sequence dataset currently presents a challenge for assembly and anchoring onto chromosomes. To provide a framework for this assembly, it would be a great advantage to have a dense map of the tammar wallaby genome. However, only limited mapping data are available for this non-model species, comprising a physical map and a linkage map. RESULTS: We combined all available tammar wallaby mapping data to create a tammar wallaby integrated map, using the Location DataBase (LDB) strategy. This first-generation integrated map combines all available information from the second-generation tammar wallaby linkage map with 148 loci, and extensive FISH mapping data for 492 loci, especially for genes likely to be located at the ends of wallaby chromosomes or at evolutionary breakpoints inferred from comparative information. For loci whose positions are only approximately known, their location in the integrated map was refined on the basis of comparative information from opossum (Monodelphis domestica) and human. Interpolation of segments from the opossum and human assemblies into the integrated map enabled the subsequent construction of a tammar wallaby first-generation virtual genome map, which comprises 14336 markers, including 13783 genes recruited from opossum and human assemblies. Both maps are freely available at http://compldb.angis.org.au. CONCLUSIONS: The first-generation integrated map and the first-generation virtual genome map provide a backbone for the chromosome assembly of the tammar wallaby genome sequence. For example, 78% of the 10257 gene-scaffolds in the Ensembl annotation of the tammar wallaby genome sequence (including 10522 protein-coding genes) can now be given a chromosome location in the tammar wallaby virtual genome map.

A second-generation anchored genetic linkage map of the tammar wallaby (Macropus eugenii).

BACKGROUND: The tammar wallaby, Macropus eugenii, a small kangaroo used for decades for studies of reproduction and metabolism, is the model Australian marsupial for genome sequencing and genetic investigations. The production of a more comprehensive cytogenetically-anchored genetic linkage map will significantly contribute to the deciphering of the tammar wallaby genome. It has great value as a resource to identify novel genes and for comparative studies, and is vital for the ongoing genome sequence assembly and gene ordering in this species. RESULTS: A second-generation anchored tammar wallaby genetic linkage map has been constructed based on a total of 148 loci. The linkage map contains the original 64 loci included in the first-generation map, plus an additional 84 microsatellite loci that were chosen specifically to increase coverage and assist with the anchoring and orientation of linkage groups to chromosomes. These additional loci were derived from (a) sequenced BAC clones that had been previously mapped to tammar wallaby chromosomes by fluorescence in situ hybridization (FISH), (b) End sequence from BACs subsequently FISH-mapped to tammar wallaby chromosomes, and (c) tammar wallaby genes orthologous to opossum genes predicted to fill gaps in the tammar wallaby linkage map as well as three X-linked markers from a published study. Based on these 148 loci, eight linkage groups were formed. These linkage groups were assigned (via FISH-mapped markers) to all seven autosomes and the X chromosome. The sex-pooled map size is 1402.4 cM, which is estimated to provide 82.6% total coverage of the genome, with an average interval distance of 10.9 cM between adjacent markers. The overall ratio of female/male map length is 0.84, which is comparable to the ratio of 0.78 obtained for the first-generation map. CONCLUSIONS: Construction of this second-generation genetic linkage map is a significant step towards complete coverage of the tammar wallaby genome and considerably extends that of the first-generation map. It will be a valuable resource for ongoing tammar wallaby genetic research and assembling the genome sequence. The sex-pooled map is available online at http://compldb.angis.org.au/.

Evolutionary transitions between mechanisms of sex determination in vertebrates.

Sex in many organisms is a dichotomous phenotype--individuals are either male or female. The molecular pathways underlying sex determination are governed by the genetic contribution of parents to the zygote, the environment in which the zygote develops or interaction of the two, depending on the species. Systems in which multiple interacting influences or a continuously varying influence (such as temperature) determines a dichotomous outcome have at least one threshold. We show that when sex is viewed as a threshold trait, evolution in that threshold can permit novel transitions between genotypic and temperature-dependent sex determination (TSD) and remarkably, between male (XX/XY) and female (ZZ/ZW) heterogamety. Transitions are possible without substantive genotypic innovation of novel sex-determining mutations or transpositions, so that the master sex gene and sex chromosome pair can be retained in ZW-XY transitions. We also show that evolution in the threshold can explain all observed patterns in vertebrate TSD, when coupled with evolution in embryonic survivorship limits.

Non-homologous sex chromosomes of birds and snakes share repetitive sequences.

Snake sex chromosomes provided Susumo Ohno with the material on which he based his theory of how sex chromosomes differentiate from autosomal pairs. Like birds, snakes have a ZZ male/ZW female sex chromosome system, in which the snake Z is a macrochromosome much the same size as the bird Z. However, the gene content shows clearly that the snake and bird Z chromosomes are completely non-homologous. The molecular aspect of W chromosome degeneration in snakes remains largely unexplored. We used comparative genomic hybridization to identify the female-specific region of the W chromosome in representative species of Australian snakes. Using this approach, we show that an increasingly complex suite of repeats accompanies the evolution of W chromosome heteromorphy. In particular, we found that while the python Liasis fuscus exhibits no sex-specific repeats and indeed, no cytologically recognizable sex-specific region, the colubrid Stegonotus cucullatus shows a large domain on the short arm of the W chromosome that consists of female-specific repeats, and the large W of Notechis scutatus is composed almost entirely of repetitive sequences, including Bkm and 18S rDNA-related elements. FISH mapping of both simple and complex probes shows patterns of repeat amplification concordant with the size of the female-specific region in each species examined. Mapping of intronic sequences of genes that are sex-linked in both birds (DMRT1) and snakes (CTNNB1) reveals massive amplification in discrete domains on the W chromosome of the elapid N. scutatus. Using chicken W chromosome paint, we demonstrate that repetitive sequences are shared between the sex chromosomes of birds and derived snakes. This could be explained by ancestral but as yet undetected shared synteny of bird and snake sex chromosomes or may indicate functional homology of the repeats and suggests that degeneration is a convergent property of sex chromosome evolution. We also establish that synteny of snake Z-linked genes has been conserved for at least 166 million years and that the snake Z consists of two conserved blocks derived from the same ancestral vertebrate chromosome.