syntenyPlotteR: a user-friendly R package to visualize genome synteny, ideal for both experienced and novice bioinformaticians.

Motivation: The rapid increase in the number of chromosome-scale genome assemblies has renewed interest in chromosome evolution studies. The visualization of syntenic relationships between genomes is a crucial initial step in the study of chromosome rearrangements and evolution. There are few tools available that serve this purpose, and they can be difficult to learn. Moreover, these tools are limited in the number of species comparisons that can be visualized and the size of chromosome rearrangements identified. Thus, the development of novel visualization tools is in strong need. Results: Here, we present syntenyPlotteR, an R package developed to visualize homologous synteny blocks in a pairwise or multispecies manner. This package contains three functions that allow users to generate publication-quality representations of syntenic relationships easily and quickly between genomes of interest. Availability and implementation: SyntenyPlotteR can be installed from CRAN with the documentation found in https://farre-lab.github.io/syntenyPlotteR/.

Patterns of chromosome evolution in ruminants.

Studying when and where gross genomic rearrangements occurred during evolution is key to understanding changes in genome structure with functional consequences that might eventually lead to speciation. Here we identified chromosome rearrangements in ruminants, a clade characterized by large chromosome differences. Using 26 genome assemblies, we reconstructed five ancestral karyotypes and classified the rearrangement events occurring in each lineage. With these reconstructions, we then identified evolutionary breakpoints regions (EBRs) and synteny fragments. Ruminant karyotype evolution is characterized by inversions, while interchromosomal rearrangements occurred preferentially in the oldest ancestor of ruminants. We found that EBRs are depleted of protein coding genes, including housekeeping genes. Similarly, EBRs are not enriched in high GC regions, suggesting that meiotic double strand breaks might not be their origin. Overall, our results characterize at fine detail the location of chromosome rearrangements in ruminant evolution and provide new insights into the formation of EBRs.

A shared 'vulnerability code' underpins varying sources of DNA damage throughout paternal germline transmission in mouse.

During mammalian spermatogenesis, the paternal genome is extensively remodelled via replacement of histones with protamines forming the highly compact mature sperm nucleus. Compaction occurs in post-meiotic spermatids and is accompanied by extensive double strand break (DSB) formation. We investigate the epigenomic and genomic context of mouse spermatid DSBs, identifying primary sequence motifs, secondary DNA structures and chromatin contexts associated with this damage. Consistent with previously published results we find spermatid DSBs positively associated with short tandem repeats and LINE elements. We further show spermatid DSBs preferentially occur in association with (CA)n, (NA)n and (RY)n repeats, in predicted Z-DNA, are not associated with G-quadruplexes, are preferentially found in regions of low histone mark coverage and engage the remodelling/NHEJ factor BRD4. Locations incurring DSBs in spermatids also show distinct epigenetic profiles throughout later developmental stages: regions retaining histones in mature sperm, regions susceptible to oxidative damage in mature sperm, and fragile two-cell like embryonic stem cell regions bound by ZSCAN4 all co-localise with spermatid DSBs and with each other. Our results point to a common 'vulnerability code' unifying several types of DNA damage occurring on the paternal genome during reproduction, potentially underpinned by torsional changes during sperm chromatin remodelling.

Principles of 3D chromosome folding and evolutionary genome reshuffling in mammals.

Studying the similarities and differences in genomic interactions between species provides fertile grounds for determining the evolutionary dynamics underpinning genome function and speciation. Here, we describe the principles of 3D genome folding in vertebrates and show how lineage-specific patterns of genome reshuffling can result in different chromatin configurations. We (1) identified different patterns of chromosome folding in across vertebrate species (centromere clustering versus chromosomal territories); (2) reconstructed ancestral marsupial and afrotherian genomes analyzing whole-genome sequences of species representative of the major therian phylogroups; (3) detected lineage-specific chromosome rearrangements; and (4) identified the dynamics of the structural properties of genome reshuffling through therian evolution. We present evidence of chromatin configurational changes that result from ancestral inversions and fusions/fissions. We catalog the close interplay between chromatin higher-order organization and therian genome evolution and introduce an interpretative hypothesis that explains how chromatin folding influences evolutionary patterns of genome reshuffling.

Evolution of the ancestral mammalian karyotype and syntenic regions.

Decrypting the rearrangements that drive mammalian chromosome evolution is critical to understanding the molecular bases of speciation, adaptation, and disease susceptibility. Using 8 scaffolded and 26 chromosome-scale genome assemblies representing 23/26 mammal orders, we computationally reconstructed ancestral karyotypes and syntenic relationships at 16 nodes along the mammalian phylogeny. Three different reference genomes (human, sloth, and cattle) representing phylogenetically distinct mammalian superorders were used to assess reference bias in the reconstructed ancestral karyotypes and to expand the number of clades with reconstructed genomes. The mammalian ancestor likely had 19 pairs of autosomes, with nine of the smallest chromosomes shared with the common ancestor of all amniotes (three still conserved in extant mammals), demonstrating a striking conservation of synteny for ∼320 My of vertebrate evolution. The numbers and types of chromosome rearrangements were classified for transitions between the ancestral mammalian karyotype, descendent ancestors, and extant species. For example, 94 inversions, 16 fissions, and 14 fusions that occurred over 53 My differentiated the therian from the descendent eutherian ancestor. The highest breakpoint rate was observed between the mammalian and therian ancestors (3.9 breakpoints/My). Reconstructed mammalian ancestor chromosomes were found to have distinct evolutionary histories reflected in their rates and types of rearrangements. The distributions of genes, repetitive elements, topologically associating domains, and actively transcribed regions in multispecies homologous synteny blocks and evolutionary breakpoint regions indicate that purifying selection acted over millions of years of vertebrate evolution to maintain syntenic relationships of developmentally important genes and regulatory landscapes of gene-dense chromosomes.

3D chromatin remodelling in the germ line modulates genome evolutionary plasticity.

Chromosome folding has profound impacts on gene regulation, whose evolutionary consequences are far from being understood. Here we explore the relationship between 3D chromatin remodelling in mouse germ cells and evolutionary changes in genome structure. Using a comprehensive integrative computational analysis, we (i) reconstruct seven ancestral rodent genomes analysing whole-genome sequences of 14 species representatives of the major phylogroups, (ii) detect lineage-specific chromosome rearrangements and (iii) identify the dynamics of the structural and epigenetic properties of evolutionary breakpoint regions (EBRs) throughout mouse spermatogenesis. Our results show that EBRs are devoid of programmed meiotic DNA double-strand breaks (DSBs) and meiotic cohesins in primary spermatocytes, but are associated in post-meiotic cells with sites of DNA damage and functional long-range interaction regions that recapitulate ancestral chromosomal configurations. Overall, we propose a model that integrates evolutionary genome reshuffling with DNA damage response mechanisms and the dynamic spatial genome organisation of germ cells.

Preimplantation Genetic Testing for Aneuploidy Improves Live Birth Rates with In Vitro Produced Bovine Embryos: A Blind Retrospective Study.

Approximately one million in vitro produced (IVP) cattle embryos are transferred worldwide each year as a way to improve the rates of genetic gain. The most advanced programmes also apply genomic selection at the embryonic stage by SNP genotyping and the calculation of genomic estimated breeding values (GEBVs). However, a high proportion of cattle embryos fail to establish a pregnancy. Here, we demonstrate that further interrogation of the SNP data collected for GEBVs can effectively remove aneuploid embryos from the pool, improving live births per embryo transfer (ET). Using three preimplantation genetic testing for aneuploidy (PGT-A) approaches, we assessed 1713 cattle blastocysts in a blind, retrospective analysis. Our findings indicate aneuploid embryos have a 5.8% chance of establishing a pregnancy and a 5.0% chance of given rise to a live birth. This compares to 59.6% and 46.7% for euploid embryos (p < 0.0001). PGT-A improved overall pregnancy and live birth rates by 7.5% and 5.8%, respectively (p < 0.0001). More detailed analyses revealed donor, chromosome, stage, grade, and sex-specific rates of error. Notably, we discovered a significantly higher incidence of aneuploidy in XY embryos and, as in humans, detected a preponderance of maternal meiosis I errors. Our data strongly support the use of PGT-A in cattle IVP programmes.

Identification of sequence changes in myosin II that adjust muscle contraction velocity.

The speed of muscle contraction is related to body size; muscles in larger species contract at slower rates. Since contraction speed is a property of the myosin isoform expressed in a muscle, we investigated how sequence changes in a range of muscle myosin II isoforms enable this slower rate of muscle contraction. We considered 798 sequences from 13 mammalian myosin II isoforms to identify any adaptation to increasing body mass. We identified a correlation between body mass and sequence divergence for the motor domain of the 4 major adult myosin II isoforms (ß/Type I, IIa, IIb, and IIx), suggesting that these isoforms have adapted to increasing body mass. In contrast, the non-muscle and developmental isoforms show no correlation of sequence divergence with body mass. Analysis of the motor domain sequence of ß-myosin (predominant myosin in Type I/slow and cardiac muscle) from 67 mammals from 2 distinct clades identifies 16 sites, out of 800, associated with body mass (padj < 0.05) but not with the clade (padj > 0.05). Both clades change the same small set of amino acids, in the same order from small to large mammals, suggesting a limited number of ways in which contraction velocity can be successfully manipulated. To test this relationship, the 9 sites that differ between human and rat were mutated in the human ß-myosin to match the rat sequence. Biochemical analysis revealed that the rat-human ß-myosin chimera functioned like the native rat myosin with a 2-fold increase in both motility and in the rate of ADP release from the actin-myosin crossbridge (the step that limits contraction velocity). Thus, these sequence changes indicate adaptation of ß-myosin as species mass increased to enable a reduced contraction velocity and heart rate.

Mitochondrial Genome Evolution, Genetic Diversity, and Population Structure in British Water Voles (Arvicola amphibius).

The European water vole (Arvicola amphibius) is a rodent within the subfamily Arvicolinae. In Britain, water voles have declined rapidly during the last century, making them a conservation priority. The relationship of Arvicola to other genera within Arvicolinae remains debated. Additionally, the impact that captive breeding programs in Britain are having on the genetic diversity of water voles is unknown. We use available mitochondrial genomes to construct the phylogeny of species within Arvicolinae, followed by sequencing the mitochondrial DNA control region of 17 individuals from a captive population of water voles in Britain to assess their genetic diversity and population structure. Our study first provides an updated phylogenetic tree of Arvicolinae using the mitochondrial genome of 31 species. Second, our results show considerable genetic diversity in the captive population of water voles, when compared with natural populations in Britain. We confirm the grouping of British water voles into two clades, with all captive individuals found in the English/Welsh clade. Moreover, captive water voles clustered closely with populations in the South East and East of England. The mitochondrial genome provides a useful marker to study the phylogenetics of this rodent clade and in addition, our study provides support for the breeding program at Wildwood Trust and provides a framework for future conservation genetics studies in this species.

Validity and reliability of the Spanish caregiver reaction assessment scale for caregivers of people with dementia.

AIM: Of the few instruments available to measure the important, positive aspects of caring, the Caregiver Reaction Assessment is regarded as the most suitable, but there is no validated Spanish version. The aim of this study was to translate the Caregiver Reaction Assessment into Spanish and assess its psychometric properties. METHODS: The Caregiver Reaction Assessment was translated into Spanish and then back translated. Validation included construct validity, criterion validity and reliability: 287 interviews were conducted with informal caregivers of patients with dementia between November 2010 and April 2012. RESULTS: Principal component analysis confirmed the original instrument's five subscales. Criterion validity showed a moderate negative correlation between the impact on health and the EQ-5D (r = -.43), that is, a greater impact on health correlated with lower health-related quality of life. The Caregiver Reaction Assessment showed good internal consistency, with a Cronbach alpha of .804, and good temporal stability for the distinct subscales, with intraclass correlation coefficients varying from .683 to .729 (p < .001). CONCLUSION: The Caregiver Reaction Assessment is a reliable, valid instrument for the measurement of the reactions of informal caregivers of patients with dementia, with good psychometric properties.

The Plasticity of Genome Architecture.

Understanding the origin of species and their adaptability to new environments is one of the main questions in biology [...].

Comparative Chromosome Mapping of Musk Ox and the X Chromosome among Some Bovidae Species.

: Bovidae, the largest family in Pecora infraorder, are characterized by a striking variability in diploid number of chromosomes between species and among individuals within a species. The bovid X chromosome is also remarkably variable, with several morphological types in the family. Here we built a detailed chromosome map of musk ox (Ovibosmoschatus), a relic species originating from Pleistocene megafauna, with dromedary and human probes using chromosome painting. We trace chromosomal rearrangements during Bovidae evolution by comparing species already studied by chromosome painting. The musk ox karyotype differs from the ancestral pecoran karyotype by six fusions, one fission, and three inversions. We discuss changes in pecoran ancestral karyotype in the light of new painting data. Variations in the X chromosome structure of four bovid species nilgai bull (Boselaphustragocamelus), saola (Pseudoryxnghetinhensis), gaur (Bosgaurus), and Kirk's Dikdik (Madoquakirkii) were further analyzed using 26 cattle BAC-clones. We found the duplication on the X in saola. We show main rearrangements leading to the formation of four types of bovid X: Bovinae type with derived cattle subtype formed by centromere reposition and Antilopinae type with Caprini subtype formed by inversion in XSB3.

Extreme heterogeneity in sex chromosome differentiation and dosage compensation in livebearers.

Once recombination is halted between the X and Y chromosomes, sex chromosomes begin to differentiate and transition to heteromorphism. While there is a remarkable variation across clades in the degree of sex chromosome divergence, far less is known about the variation in sex chromosome differentiation within clades. Here, we combined whole-genome and transcriptome sequencing data to characterize the structure and conservation of sex chromosome systems across Poeciliidae, the livebearing clade that includes guppies. We found that the Poecilia reticulata XY system is much older than previously thought, being shared not only with its sister species, Poecilia wingei, but also with Poecilia picta, which diverged roughly 20 million years ago. Despite the shared ancestry, we uncovered an extreme heterogeneity across these species in the proportion of the sex chromosome with suppressed recombination, and the degree of Y chromosome decay. The sex chromosomes in P. reticulata and P. wingei are largely homomorphic, with recombination in the former persisting over a substantial fraction. However, the sex chromosomes in P. picta are completely nonrecombining and strikingly heteromorphic. Remarkably, the profound degradation of the ancestral Y chromosome in P. picta is counterbalanced by the evolution of functional chromosome-wide dosage compensation in this species, which has not been previously observed in teleost fish. Our results offer important insight into the initial stages of sex chromosome evolution and dosage compensation.

An integrated chromosome-scale genome assembly of the Masai giraffe (Giraffa camelopardalis tippelskirchi).

BACKGROUND: The Masai giraffe (Giraffa camelopardalis tippelskirchi) is the largest-bodied giraffe and the world's tallest terrestrial animal. With its extreme size and height, the giraffe's unique anatomical and physiological adaptations have long been of interest to diverse research fields. Giraffes are also critical to ecosystems of sub-Saharan Africa, with their long neck serving as a conduit to food sources not shared by other herbivores. Although the genome of a Masai giraffe has been sequenced, the assembly was highly fragmented and suboptimal for genome analysis. Herein we report an improved giraffe genome assembly to facilitate evolutionary analysis of the giraffe and other ruminant genomes. FINDINGS: Using SOAPdenovo2 and 170 Gbp of Illumina paired-end and mate-pair reads, we generated a 2.6-Gbp male Masai giraffe genome assembly, with a scaffold N50 of 3 Mbp. The incorporation of 114.6 Gbp of Chicago library sequencing data resulted in a HiRise SOAPdenovo + Chicago assembly with an N50 of 48 Mbp and containing 95% of expected genes according to BUSCO analysis. Using the Reference-Assisted Chromosome Assembly tool, we were able to order and orient scaffolds into 42 predicted chromosome fragments (PCFs). Using fluorescence in situ hybridization, we placed 153 cattle bacterial artificial chromosomes onto giraffe metaphase spreads to assess and assign the PCFs on 14 giraffe autosomes and the X chromosome resulting in the final assembly with an N50 of 177.94 Mbp. In this assembly, 21,621 protein-coding genes were identified using both de novo and homology-based predictions. CONCLUSIONS: We have produced the first chromosome-scale genome assembly for a Giraffidae species. This assembly provides a valuable resource for the study of artiodactyl evolution and for understanding the molecular basis of the unique adaptive traits of giraffes. In addition, the assembly will provide a powerful resource to assist conservation efforts of Masai giraffe, whose population size has declined by 52% in recent years.

The Red Fox Y-Chromosome in Comparative Context.

While the number of mammalian genome assemblies has proliferated, Y-chromosome assemblies have lagged behind. This discrepancy is caused by biological features of the Y-chromosome, such as its high repeat content, that present challenges to assembly with short-read, next-generation sequencing technologies. Partial Y-chromosome assemblies have been developed for the cat (Feliscatus), dog (Canislupusfamiliaris), and grey wolf (Canislupuslupus), providing the opportunity to examine the red fox (Vulpesvulpes) Y-chromosome in the context of closely related species. Here we present a data-driven approach to identifying Y-chromosome sequence among the scaffolds that comprise the short-read assembled red fox genome. First, scaffolds containing genes found on the Y-chromosomes of cats, dogs, and wolves were identified. Next, analysis of the resequenced genomes of 15 male and 15 female foxes revealed scaffolds containing male-specific k-mers and patterns of inter-sex copy number variation consistent with the heterogametic chromosome. Analyzing variation across these two metrics revealed 171 scaffolds containing 3.37 Mbp of putative Y-chromosome sequence. The gene content of these scaffolds is consistent overall with that of the Y-chromosome in other carnivore species, though the red fox Y-chromosome carries more copies of BCORY2 and UBE1Y than has been reported in related species and fewer copies of SRY than in other canids. The assignment of these scaffolds to the Y-chromosome serves to further characterize the content of the red fox draft genome while providing resources for future analyses of canid Y-chromosome evolution.

Evolution of gene regulation in ruminants differs between evolutionary breakpoint regions and homologous synteny blocks.

The role of chromosome rearrangements in driving evolution has been a long-standing question of evolutionary biology. Here we focused on ruminants as a model to assess how rearrangements may have contributed to the evolution of gene regulation. Using reconstructed ancestral karyotypes of Cetartiodactyls, Ruminants, Pecorans, and Bovids, we traced patterns of gross chromosome changes. We found that the lineage leading to the ruminant ancestor after the split from other cetartiodactyls was characterized by mostly intrachromosomal changes, whereas the lineage leading to the pecoran ancestor (including all livestock ruminants) included multiple interchromosomal changes. We observed that the liver cell putative enhancers in the ruminant evolutionary breakpoint regions are highly enriched for DNA sequences under selective constraint acting on lineage-specific transposable elements (TEs) and a set of 25 specific transcription factor (TF) binding motifs associated with recently active TEs. Coupled with gene expression data, we found that genes near ruminant breakpoint regions exhibit more divergent expression profiles among species, particularly in cattle, which is consistent with the phylogenetic origin of these breakpoint regions. This divergence was significantly greater in genes with enhancers that contain at least one of the 25 specific TF binding motifs and located near bovidae-to-cattle lineage breakpoint regions. Taken together, by combining ancestral karyotype reconstructions with analysis of cis regulatory element and gene expression evolution, our work demonstrated that lineage-specific regulatory elements colocalized with gross chromosome rearrangements may have provided valuable functional modifications that helped to shape ruminant evolution.

A Near Chromosome Assembly of the Dromedary Camel Genome.

The dromedary camel is an economically and socially important species of livestock in many parts of the world, being used for transport and the production of milk and meat. Much like cattle and horses, the camel may be found in industrial farming conditions as well as used in sporting. Camel racing is a multi-million dollar industry, with some specimens being valued at upward of 9.5 million USD. Despite its apparent value to humans, the dromedary camel is a neglected species in genomics. While cattle and other domesticated species have had much attention in terms of genome assembly, the camel has only been assembled to scaffold level, which does not give a clear indication of the order or chromosomal location of sequenced fragments. In this study, the Reference Assistant Chromosome Assembly (RACA) algorithm was implemented to use read-pair information of camel scaffolds, aligned with the cattle and human genomes in order to organize and orient these scaffolds in a near-chromosome level assembly. This method generated 72 large size fragments (N50 54.36 Mb). These predicted chromosome fragments (PCFs) were then compared with comparative maps of camel and cytogenetic map of alpaca chromosomes, allowing us to further upgrade the assembly. This dromedary camel assembly will be an invaluable tool to verify future camel assemblies generated with chromatin conformation or/and long read technologies. This study provides the first near-chromosome assembly of the dromedary camel, thus adding this economically important species to a growing pool of knowledge regarding the genome structure of domesticated livestock.

A near-chromosome-scale genome assembly of the gemsbok (Oryx gazella): an iconic antelope of the Kalahari desert.

Background: The gemsbok (Oryx gazella) is one of the largest antelopes in Africa. Gemsbok are heterothermic and thus highly adapted to live in the desert, changing their feeding behavior when faced with extreme drought and heat. A high-quality genome sequence of this species will assist efforts to elucidate these and other important traits of gemsbok and facilitate research on conservation efforts. Findings: Using 180 Gbp of Illumina paired-end and mate-pair reads, a 2.9 Gbp assembly with scaffold N50 of 1.48 Mbp was generated using SOAPdenovo. Scaffolds were extended using Chicago library sequencing, which yielded an additional 114.7 Gbp of DNA sequence. The HiRise assembly using SOAPdenovo + Chicago library sequencing produced a scaffold N50 of 47 Mbp and a final genome size of 2.9 Gbp, representing 90.6% of the estimated genome size and including 93.2% of expected genes according to Benchmarking Universal Single-Copy Orthologs analysis. The Reference-Assisted Chromosome Assembly tool was used to generate a final set of 47 predicted chromosome fragments with N50 of 86.25 Mbp and containing 93.8% of expected genes. A total of 23,125 protein-coding genes and 1.14 Gbp of repetitive sequences were annotated using de novo and homology-based predictions. Conclusions: Our results provide the first high-quality, chromosome-scale genome sequence assembly for gemsbok, which will be a valuable resource for studying adaptive evolution of this species and other ruminants.

Reconstruction of avian ancestral karyotypes reveals differences in the evolutionary history of macro- and microchromosomes.

BACKGROUND: Reconstruction of ancestral karyotypes is critical for our understanding of genome evolution, allowing for the identification of the gross changes that shaped extant genomes. The identification of such changes and their time of occurrence can shed light on the biology of each species, clade and their evolutionary history. However, this is impeded by both the fragmented nature of the majority of genome assemblies and the limitations of the available software to work with them. These limitations are particularly apparent in birds, with only 10 chromosome-level assemblies reported thus far. Algorithmic approaches applied to fragmented genome assemblies can nonetheless help define patterns of chromosomal change in defined taxonomic groups. RESULTS: Here, we make use of the DESCHRAMBLER algorithm to perform the first large-scale study of ancestral chromosome structure and evolution in birds. This algorithm allows us to reconstruct the overall genome structure of 14 key nodes of avian evolution from the Avian ancestor to the ancestor of the Estrildidae, Thraupidae and Fringillidae families. CONCLUSIONS: Analysis of these reconstructions provides important insights into the variability of rearrangement rates during avian evolution and allows the detection of patterns related to the chromosome distribution of evolutionary breakpoint regions. Moreover, the inclusion of microchromosomes in our reconstructions allows us to provide novel insights into the evolution of these avian chromosomes, specifically.

Chromosome-level assembly reveals extensive rearrangement in saker falcon and budgerigar, but not ostrich, genomes.

BACKGROUND: The number of de novo genome sequence assemblies is increasing exponentially; however, relatively few contain one scaffold/contig per chromosome. Such assemblies are essential for studies of genotype-to-phenotype association, gross genomic evolution, and speciation. Inter-species differences can arise from chromosomal changes fixed during evolution, and we previously hypothesized that a higher fraction of elements under negative selection contributed to avian-specific phenotypes and avian genome organization stability. The objective of this study is to generate chromosome-level assemblies of three avian species (saker falcon, budgerigar, and ostrich) previously reported as karyotypically rearranged compared to most birds. We also test the hypothesis that the density of conserved non-coding elements is associated with the positions of evolutionary breakpoint regions. RESULTS: We used reference-assisted chromosome assembly, PCR, and lab-based molecular approaches, to generate chromosome-level assemblies of the three species. We mapped inter- and intrachromosomal changes from the avian ancestor, finding no interchromosomal rearrangements in the ostrich genome, despite it being previously described as chromosomally rearranged. We found that the average density of conserved non-coding elements in evolutionary breakpoint regions is significantly reduced. Fission evolutionary breakpoint regions have the lowest conserved non-coding element density, and intrachromomosomal evolutionary breakpoint regions have the highest. CONCLUSIONS: The tools used here can generate inexpensive, efficient chromosome-level assemblies, with > 80% assigned to chromosomes, which is comparable to genomes assembled using high-density physical or genetic mapping. Moreover, conserved non-coding elements are important factors in defining where rearrangements, especially interchromosomal, are fixed during evolution without deleterious effects.