Reconstruction and evolutionary history of eutherian chromosomes

Genome architecture of the allotetraploid wild grass Aegilops ventricosa reveals its evolutionary history and contributions to wheat improvement

The allotetraploid wild grass Aegilops ventricosa (2n = 4x = 28, genome DvDvNvNv) has been recognized as an important germplasm resource for wheat improvement owing to its ability to tolerate biotic stresses. In particular, the 2NvS segment from Ae. ventricosa, as a stable and effective resistance source, has contributed greatly to wheat improvement. The 2NvS/2AS translocation is a prevalent chromosomal translocation between common wheat and wild relatives, ranking just behind the 1B/1R translocation in importance for modern wheat breeding. Here, we assembled a high-quality chromosome-level reference genome of Ae. ventricosa RM271 with a total length of 8.67 Gb. Phylogenomic analyses revealed that the progenitor of the Dv subgenome of Ae. ventricosa is Ae. tauschii ssp. tauschii (genome DD); by contrast, the progenitor of the D subgenome of bread wheat (Triticum aestivum L.) is Ae. tauschii ssp. strangulata (genome DD). The oldest polyploidization time of Ae. ventricosa occurred ∼0.7 mya. The Dv subgenome of Ae. ventricosa is less conserved than the D subgenome of bread wheat. Construction of a graph-based pangenome of 2AS/6NvL (originally known as 2NvS) segments from Ae. ventricosa and other genomes in the Triticeae enabled us to identify candidate resistance genes sourced from Ae. ventricosa. We identified 12 nonredundant introgressed segments from the Dv and Nv subgenomes using a large winter wheat collection representing the full diversity of the European wheat genetic pool, and 29.40% of European wheat varieties inherit at least one of these segments. The high-quality RM271 reference genome will provide a basis for cloning key genes, including the Yr17-Lr37-Sr38-Cre5 resistance gene cluster in Ae. ventricosa, and facilitate the full use of elite wild genetic resources to accelerate wheat improvement.

Compensation of gene dosage on the mammalian X

Changes in gene dosage can have tremendous evolutionary potential (e.g. whole-genome duplications), but without compensatory mechanisms, they can also lead to gene dysregulation and pathologies. Sex chromosomes are a paradigmatic example of naturally occurring gene dosage differences and their compensation. In species with chromosome-based sex determination, individuals within the same population necessarily show 'natural' differences in gene dosage for the sex chromosomes. In this Review, we focus on the mammalian X chromosome and discuss recent new insights into the dosage-compensation mechanisms that evolved along with the emergence of sex chromosomes, namely X-inactivation and X-upregulation. We also discuss the evolution of the genetic loci and molecular players involved, as well as the regulatory diversity and potentially different requirements for dosage compensation across mammalian species.

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

Counting Sorting Scenarios and Intermediate Genomes for the Rank Distance

An important problem in genome comparison is the genome sorting problem, that is, the problem of finding a sequence of basic operations that transforms one genome into another whose length (possibly weighted) equals the distance between them. These sequences are called optimal sorting scenarios. However, there is usually a large number of such scenarios, and a naïve algorithm is very likely to be biased towards a specific type of scenario, impairing its usefulness in real-world applications. One way to go beyond the traditional sorting algorithms is to explore all possible solutions, looking at all the optimal sorting scenarios instead of just an arbitrary one. Another related approach is to analyze all the intermediate genomes, that is, all the genomes that can occur in an optimal sorting scenario. In this article, we show how to enumerate the optimal sorting scenarios and the intermediate genomes between any two given genomes, under the rank distance.

A conserved NR5A1-responsive enhancer regulates SRY in testis-determination

The Y-linked SRY gene initiates mammalian testis-determination. However, how the expression of SRY is regulated remains elusive. Here, we demonstrate that a conserved steroidogenic factor-1 (SF-1)/NR5A1 binding enhancer is required for appropriate SRY expression to initiate testis-determination in humans. Comparative sequence analysis of SRY 5' regions in mammals identified an evolutionary conserved SF-1/NR5A1-binding motif within a 250 bp region of open chromatin located 5 kilobases upstream of the SRY transcription start site. Genomic analysis of 46,XY individuals with disrupted testis-determination, including a large multigenerational family, identified unique single-base substitutions of highly conserved residues within the SF-1/NR5A1-binding element. In silico modelling and in vitro assays demonstrate the enhancer properties of the NR5A1 motif. Deletion of this hemizygous element by genome-editing, in a novel in vitro cellular model recapitulating human Sertoli cell formation, resulted in a significant reduction in expression of SRY. Therefore, human NR5A1 acts as a regulatory switch between testis and ovary development by upregulating SRY expression, a role that may predate the eutherian radiation. We show that disruption of an enhancer can phenocopy variants in the coding regions of SRY that cause human testis dysgenesis. Since disease causing variants in enhancers are currently rare, the regulation of gene expression in testis-determination offers a paradigm to define enhancer activity in a key developmental process.

A Bird's-Eye View of Chromosomic Evolution in the Class Aves

Birds (Aves) are the most speciose of terrestrial vertebrates, displaying Class-specific characteristics yet incredible external phenotypic diversity. Critical to agriculture and as model organisms, birds have adapted to many habitats. The only extant examples of dinosaurs, birds emerged ~150 mya and >10% are currently threatened with extinction. This review is a comprehensive overview of avian genome ("chromosomic") organization research based mostly on chromosome painting and BAC-based studies. We discuss traditional and contemporary tools for reliably generating chromosome-level assemblies and analyzing multiple species at a higher resolution and wider phylogenetic distance than previously possible. These results permit more detailed investigations into inter- and intrachromosomal rearrangements, providing unique insights into evolution and speciation mechanisms. The 'signature' avian karyotype likely arose ~250 mya and remained largely unchanged in most groups including extinct dinosaurs. Exceptions include Psittaciformes, Falconiformes, Caprimulgiformes, Cuculiformes, Suliformes, occasional Passeriformes, Ciconiiformes, and Pelecaniformes. The reasons for this remarkable conservation may be the greater diploid chromosome number generating variation (the driver of natural selection) through a greater possible combination of gametes and/or an increase in recombination rate. A deeper understanding of avian genomic structure permits the exploration of fundamental biological questions pertaining to the role of evolutionary breakpoint regions and homologous synteny blocks.

A chromosome-level genome of electric catfish (Malapterurus electricus) provided new insights into order Siluriformes evolution

The electric catfish (Malapterurus electricus), belonging to the family Malapteruridae, order Siluriformes (Actinopterygii: Ostariophysi), is one of the six branches that has independently evolved electrical organs. We assembled a 796.75 Mb M. electricus genome and anchored 88.72% sequences into 28 chromosomes. Gene family analysis revealed 295 expanded gene families that were enriched on functions related to glutamate receptors. Convergent evolutionary analyses of electric organs among different lineage of electric fishes further revealed that the coding gene of rho guanine nucleotide exchange factor 4-like (arhgef4), which is associated with G-protein coupled receptor (GPCR) signaling pathway, underwent adaptive parallel evolution. Gene identification suggests visual degradation in catfishes, and an important role for taste in environmental adaptation. Our findings fill in the genomic data for a branch of electric fish and provide a relevant genetic basis for the adaptive evolution of Siluriformes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00197-8.

OMA orthology in 2024: improved prokaryote coverage, ancestral and extant GO enrichment, a revamped synteny viewer and more in the OMA Ecosystem

In this update paper, we present the latest developments in the OMA browser knowledgebase, which aims to provide high-quality orthology inferences and facilitate the study of gene families, genomes and their evolution. First, we discuss the addition of new species in the database, particularly an expanded representation of prokaryotic species. The OMA browser now offers Ancestral Genome pages and an Ancestral Gene Order viewer, allowing users to explore the evolutionary history and gene content of ancestral genomes. We also introduce a revamped Local Synteny Viewer to compare genomic neighborhoods across both extant and ancestral genomes. Hierarchical Orthologous Groups (HOGs) are now annotated with Gene Ontology annotations, and users can easily perform extant or ancestral GO enrichments. Finally, we recap new tools in the OMA Ecosystem, including OMAmer for proteome mapping, OMArk for proteome quality assessment, OMAMO for model organism selection and Read2Tree for phylogenetic species tree construction from reads. These new features provide exciting opportunities for orthology analysis and comparative genomics. OMA is accessible at https://omabrowser.org.

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.

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/.

Emergence and influence of sequence bias in evolutionarily malleable, mammalian tandem arrays

BACKGROUND

The radiation of mammals at the extinction of the dinosaurs produced a plethora of new forms-as diverse as bats, dolphins, and elephants-in only 10-20 million years. Behind the scenes, adaptation to new niches is accompanied by extensive innovation in large families of genes that allow animals to contact the environment, including chemosensors, xenobiotic enzymes, and immune and barrier proteins. Genes in these "outward-looking" families are allelically diverse among humans and exhibit tissue-specific and sometimes stochastic expression.

RESULTS

Here, we show that these tandem arrays of outward-looking genes occupy AT-biased isochores and comprise the "tissue-specific" gene class that lack CpG islands in their promoters. Models of mammalian genome evolution have not incorporated the sharply different functions and transcriptional patterns of genes in AT- versus GC-biased regions. To examine the relationship between gene family expansion, sequence content, and allelic diversity, we use population genetic data and comparative analysis. First, we find that AT bias can emerge during evolutionary expansion of gene families in cis. Second, human genes in AT-biased isochores or with GC-poor promoters experience relatively low rates of de novo point mutation today but are enriched for non-synonymous variants. Finally, we find that isochores containing gene clusters exhibit low rates of recombination.

CONCLUSIONS

Our analyses suggest that tolerance of non-synonymous variation and low recombination are two forces that have produced the depletion of GC bases in outward-facing gene arrays. In turn, high AT content exerts a profound effect on their chromatin organization and transcriptional regulation.

Establishment of Primary Adult Skin Fibroblast Cell Lines from African Savanna Elephants (Loxodonta africana)

Following population declines of the African savanna elephant (Loxodonta africana) across the African continent, the establishment of primary cell lines of endangered wildlife species is paramount for the preservation of their genetic resources. In addition, it allows molecular and functional studies on the cancer suppression mechanisms of elephants, which have previously been linked to a redundancy of tumor suppressor gene TP53. This methodology describes the establishment of primary elephant dermal fibroblast (EDF) cell lines from skin punch biopsy samples (diameter: ±4 mm) of African savanna elephants (n = 4, 14-35 years). The applied tissue collection technique is minimally invasive and paves the way for future remote biopsy darting. On average, the first explant outgrowth was observed after 15.75 ± 6.30 days. The average doubling time (Td) was 93.02 ± 16.94 h and 52.39 ± 0.46 h at passage 1 and 4, respectively. Metaphase spreads confirmed the diploid number of 56 chromosomes. The successful establishment of EDF cell lines allows for future elephant cell characterization studies and for research on the cancer resistance mechanisms of elephants, which can be harnessed for human cancer prevention and treatment and contributes to the conservation of their genetic material.

TruEst: a better estimator of evolutionary distance under the INFER model

Genome rearrangements are evolutionary events that shuffle genomic architectures. The number of genome rearrangements that happened between two genomes is often used as the evolutionary distance between these species. This number is often estimated as the minimum number of genome rearrangements required to transform one genome into another which are only reliable for closely-related genomes. These estimations often underestimate the evolutionary distance for genomes that have substantially evolved from each other, and advanced statistical methods can be used to improve accuracy. Several statistical estimators have been developed, under various evolutionary models, of which the most complete one, INFER, takes into account different degrees of genome fragility. We present TruEst-an efficient tool that estimates the evolutionary distance between the genomes under the INFER model of genome rearrangements. We apply our method to both simulated and real data. It shows high accuracy on the simulated data. On the real datasets of mammal genomes the method found several pairs of genomes for which the estimated distances are in high consistency with the previous ancestral reconstruction studies.

Large-scale snake genome analyses provide insights into vertebrate development

Snakes are a remarkable squamate lineage with unique morphological adaptations, especially those related to the evolution of vertebrate skeletons, organs, and sensory systems. To clarify the genetic underpinnings of snake phenotypes, we assembled and analyzed 14 de novo genomes from 12 snake families. We also investigated the genetic basis of the morphological characteristics of snakes using functional experiments. We identified genes, regulatory elements, and structural variations that have potentially contributed to the evolution of limb loss, an elongated body plan, asymmetrical lungs, sensory systems, and digestive adaptations in snakes. We identified some of the genes and regulatory elements that might have shaped the evolution of vision, the skeletal system and diet in blind snakes, and thermoreception in infrared-sensitive snakes. Our study provides insights into the evolution and development of snakes and vertebrates.

A genomic timescale for placental mammal evolution

The precise pattern and timing of speciation events that gave rise to all living placental mammals remain controversial. We provide a comprehensive phylogenetic analysis of genetic variation across an alignment of 241 placental mammal genome assemblies, addressing prior concerns regarding limited genomic sampling across species. We compared neutral genome-wide phylogenomic signals using concatenation and coalescent-based approaches, interrogated phylogenetic variation across chromosomes, and analyzed extensive catalogs of structural variants. Interordinal relationships exhibit relatively low rates of phylogenomic conflict across diverse datasets and analytical methods. Conversely, X-chromosome versus autosome conflicts characterize multiple independent clades that radiated during the Cenozoic. Genomic time trees reveal an accumulation of cladogenic events before and immediately after the Cretaceous-Paleogene (K-Pg) boundary, implying important roles for Cretaceous continental vicariance and the K-Pg extinction in the placental radiation.

Experimental Setups for In Vitro Studies on Radon Exposure in Mammalian Cells-A Critical Overview

Naturally occurring radon and its short lived progeny are the second leading cause of lung cancer after smoking, and the main risk factor for non-smokers. The radon progeny, mainly Polonium-218 (²¹⁸Po) and Polonium-214 (²¹⁴Po), are responsible for the highest dose deposition in the bronchial epithelium via alpha-decay. These alpha-particles release a large amount of energy over a short penetration range, which results in severe and complex DNA damage. In order to unravel the underlying biological mechanisms which are triggered by this complex DNA damage and eventually give rise to carcinogenesis, in vitro radiobiology experiments on mammalian cells have been performed using radon exposure setups, or radon analogues, which mimic alpha-particle exposure. This review provides an overview of the different experimental setups, which have been developed and used over the past decades for in vitro radon experiments. In order to guarantee reliable results, the design and dosimetry of these setups require careful consideration, which will be emphasized in this work. Results of these in vitro experiments, particularly on bronchial epithelial cells, can provide valuable information on biomarkers, which can assist to identify exposures, as well as to study the effects of localized high dose depositions and the heterogeneous dose distribution of radon.

Reconstruction of hundreds of reference ancestral genomes across the eukaryotic kingdom

Ancestral sequence reconstruction is a fundamental aspect of molecular evolution studies and can trace small-scale sequence modifications through the evolution of genomes and species. In contrast, fine-grained reconstructions of ancestral genome organizations are still in their infancy, limiting our ability to draw comprehensive views of genome and karyotype evolution. Here we reconstruct the detailed gene contents and organizations of 624 ancestral vertebrate, plant, fungi, metazoan and protist genomes, 183 of which are near-complete chromosomal gene order reconstructions. Reconstructed ancestral genomes are similar to their descendants in terms of gene content as expected and agree precisely with reference cytogenetic and in silico reconstructions when available. By comparing successive ancestral genomes along the phylogenetic tree, we estimate the intra- and interchromosomal rearrangement history of all major vertebrate clades at high resolution. This freely available resource introduces the possibility to follow evolutionary processes at genomic scales in chronological order, across multiple clades and without relying on a single extant species as reference.

Chromosome-level genome assembly for takin (Budorcas taxicolor) provides insights into its taxonomic status and genetic diversity

The takin (Budorcas taxicolor) is one of the largest bovid herbivores in the subfamily Caprinae. The takin is at high risk of extinction, but its taxonomic status and genetic diversity remain unclear. In this study, we constructed the first reference genome of Bu. taxicolor using PacBio long High-Fidelity reads and Hi-C technology. The assembled genome is ~2.95 Gb with a contig N50 of 68.05 Mb, which were anchored onto 25+XY chromosomes. We found that the takin was more closely related to muskox than to other Caprinae species. Compared to the common ancestral karyotype of bovidae (2n = 60), we found the takin (2n = 52) experienced four chromosome fusions and one large translocation. Furthermore, we resequenced nine golden takins from the main distribution area, the Qinling Mountains, and identified 3.3 million single nucleotide polymorphisms. The genetic diversity of takin was very low (θπ = 0.00028 and heterozygosity =0.00038), among the lowest detected in domestic and wild mammals. Takin genomes showed a high inbreeding coefficient (FROH =0.217), suggesting severe inbreeding depression. The demographic history showed that the effective population size of takins declined significantly from ~100,000 years ago. Our results provide valuable information for protection of takins and insights into their evolution.

Maps of Constitutive-Heterochromatin Distribution for Four Martes Species (Mustelidae, Carnivora, Mammalia) Show the Formative Role of Macrosatellite Repeats in Interspecific Variation of Chromosome Structure

Constitutive-heterochromatin placement in the genome affects chromosome structure by occupying centromeric areas and forming large blocks. To investigate the basis for heterochromatin variation in the genome, we chose a group of species with a conserved euchromatin part: the genus Martes [stone marten (M. foina, 2n = 38), sable (M. zibellina, 2n = 38), pine marten (M. martes, 2n = 38), and yellow-throated marten (M. flavigula, 2n = 40)]. We mined the stone marten genome for the most abundant tandem repeats and selected the top 11 macrosatellite repetitive sequences. Fluorescent in situ hybridization revealed distributions of the tandemly repeated sequences (macrosatellites, telomeric repeats, and ribosomal DNA). We next characterized the AT/GC content of constitutive heterochromatin by CDAG (Chromomycin A3-DAPI-after G-banding). The euchromatin conservatism was shown by comparative chromosome painting with stone marten probes in newly built maps of the sable and pine marten. Thus, for the four Martes species, we mapped three different types of tandemly repeated sequences critical for chromosome structure. Most macrosatellites are shared by the four species with individual patterns of amplification. Some macrosatellites are specific to a species, autosomes, or the X chromosome. The variation of core macrosatellites and their prevalence in a genome are responsible for the species-specific variation of the heterochromatic blocks.

Genome-wide identification and expression analysis of diacylglycerol acyltransferase genes in soybean (Glycine max)

Background

Soybean (Glycine max) is a major protein and vegetable oil source. In plants, diacylglycerol acyltransferase (DGAT) can exert strong flux control, which is rate-limiting for triacylglycerol biosynthesis in seed oil formation.

Methods

Here, we identified soybean DGAT genes via a bioinformatics method, thereby laying a solid foundation for further research on their function. Based on our bioinformatics analyses, including gene structure, protein domain characteristics, and phylogenetic analysis, 26 DGAT putative gene family members unevenly distributed on 12 of the 20 soybean chromosomes were identified and divided into the following four groups: DGAT1, DGAT2, WS/DGAT, and cytoplasmic DGAT.

Results

The Ka/Ks ratio of most of these genes indicated a significant positive selection pressure. DGAT genes exhibited characteristic expression patterns in soybean tissues. The differences in the structure and expression of soybean DGAT genes revealed the diversity of their functions and the complexity of soybean fatty acid metabolism. Our findings provide important information for research on the fatty acid metabolism pathway in soybean. Furthermore, our results will help identify candidate genes for potential fatty acid-profile modifications to improve soybean seed oil content.

Conclusions

This is the first time that in silico studies have been used to report the genomic and proteomic characteristics of DGAT in soybean and the effect of its specific expression on organs, age, and stages.

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.

SYNPHONI: scale-free and phylogeny-aware reconstruction of synteny conservation and transformation across animal genomes

SUMMARY

Current approaches detect conserved genomic order either at chromosomal (macrosynteny) or at subchromosomal scales (microsynteny). The latter generally requires collinearity and hard thresholds on syntenic region size, thus excluding a major proportion of syntenies with recent expansions or minor rearrangements. 'SYNPHONI' bridges the gap between micro- and macrosynteny detection, providing detailed information on both synteny conservation and transformation throughout the evolutionary history of animal genomes.

AVAILABILITY AND IMPLEMENTATION

Source code is freely available at https://github.com/nsmro/SYNPHONI, implemented in Python 3.9.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

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.

The genes from the pseudoautosomal region 1 (PAR1) of the mammalian sex chromosomes: Synteny, phylogeny and selection

Sex chromosomes recombine restrictly in their homologous area, the pseudoautosomal region (PAR), represented by PAR1 and PAR2, which behave like an autosome in both pairing and recombination. The PAR1, common to most of the eutherian mammals, is located at the terminus of the sex chromosomes short arm and exhibit recombination rates ~20 times higher than the autosomes. Here, we assessed the interspecific evolutionary genomic dynamics of 15 genes of the PAR1 across 41 mammalian genera (representing six orders). The strong negative selection detected in most of the assessed groups reinforces the presence of evolutionary constraints, imposed by the important function of the PAR1 genes. Indeed, mutations in these genes are associated with various diseases in humans, including stature problems (Klinefelter Syndrome), leukemia and mental diseases. Yet, a few genes exhibiting positive selection (ω-value >1) were depicted in Rodentia (ASMT and ZBED1) and Primates (CRLF2 and CSF2RA). Rodents have the smallest described PAR1, while that of simian primates/humans underwent a 3 to 5 fold size reduction. The assessment of the PAR1 genes synteny revealed differences among the mammalian species, especially in the Rodentia order where chromosomic translocations from the sex chromosomes to the autosomes were observed. Such syntenic changes may be an evidence of the rapid evolution in rodents, as previous referred in other papers, also depicted by their increased branch lengths in the phylogenetic analyses. Concluding, we suggest that genome migration is an important factor influencing the evolution of mammals and may result in changes of the selective pressures operating on the genome.

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.

The role of genome folding in the heritability and evolvability of structural variations is not well understood. Here the authors investigate the impact of the three-dimensional genome topology of germ cells in the formation and transmission of gross structural genomic changes detected from comparing whole-genome sequences of 14 rodent species.

Conservation of chromatin conformation in carnivores

We found the three-dimensional (3D) structure of chromatin at the chromosome level to be highly conserved for more than 50 million y of carnivore evolution. Intrachromosomal contacts were maintained even after chromosome rearrangements within carnivore lineages, demonstrating that the maintenance of 3D chromatin architecture is essential for conserved genome functions. These discoveries enabled the identification of orthologous chromosomal DNA segments among related species, a method we call 3D comparative scaffotyping. The method has application for putative chromosome assignment of chromosome-scale DNA sequence scaffolds produced by de novo genome sequencing. Broadly applied to biodiversity genome sequencing efforts, the approach can reduce costs associated with karyotyping and the physical mapping of DNA segments to chromosomes.

High throughput chromatin conformation capture (Hi-C) of leukocyte DNA was used to investigate the evolutionary stability of chromatin conformation at the chromosomal level in 11 species from three carnivore families: Felidae, Canidae, and Ursidae. Chromosome-scale scaffolds (C-scaffolds) of each species were initially used for whole-genome alignment to a reference genome within each family. This approach established putative orthologous relationships between C-scaffolds among the different species. Hi-C contact maps for all C-scaffolds were then visually compared and found to be distinct for a given reference chromosome or C-scaffold within a species and indistinguishable for orthologous C-scaffolds having a 1:1 relationship within a family. The visual patterns within families were strongly supported by eigenvectors from the Hi-C contact maps. Analysis of Hi-C contact maps and eigenvectors across the three carnivore families revealed that most cross-family orthologous subchromosomal fragments have a conserved three-dimensional (3D) chromatin structure and thus have been under strong evolutionary constraint for ∼54 My of carnivore evolution. The most pronounced differences in chromatin conformation were observed for the X chromosome and the red fox genome, whose chromosomes have undergone extensive rearrangements relative to other canids. We also demonstrate that Hi-C contact map pattern analysis can be used to accurately identify orthologous relationships between C-scaffolds and chromosomes, a method we termed “3D comparative scaffotyping.” This method provides a powerful means for estimating karyotypes in de novo sequenced species that have unknown karyotype and no physical mapping information.

Decoding sex: Elucidating sex determination and how high-quality genome assemblies are untangling the evolutionary dynamics of sex chromosomes

Sexual reproduction is a diverse and widespread process. In gonochoristic species, the differentiation of sexes occurs through diverse mechanisms, influenced by environmental and genetic factors. In most vertebrates, a master-switch gene is responsible for triggering a sex determination network. However, only a few genes have acquired master-switch functions, and this process is associated with the evolution of sex-chromosomes, which have a significant influence in evolution. Additionally, their highly repetitive regions impose challenges for high-quality sequencing, even using high-throughput, state-of-the-art techniques. Here, we review the mechanisms involved in sex determination and their role in the evolution of species, particularly vertebrates, focusing on sex chromosomes and the challenges involved in sequencing these genomic elements. We also address the improvements provided by the growth of sequencing projects, by generating a massive number of near-gapless, telomere-to-telomere, chromosome-level, phased assemblies, increasing the number and quality of sex-chromosome sequences available for further studies.

Reconstruction of the Doradinae (Siluriformes-Doradidae) ancestral diploid number and NOR pattern reveals new insights about the karyotypic diversification of the Neotropical thorny catfishes

Doradinae (Siluriformes: Doradidae) is the most species-rich subfamily among thorny catfishes, encompassing over 77 valid species, found mainly in Amazon and Platina hydrographic basins. Here, we analyzed seven Doradinae species using combined methods (e.g., cytogenetic tools and Mesquite ancestral reconstruction software) in order to scrutinize the processes that mediated the karyotype diversification in this subfamily. Our ancestral reconstruction recovered that 2n=58 chromosomes and simple nucleolar organizer regions (NOR) are ancestral features only for Wertheimerinae and the most clades of Doradinae. Some exceptions were found in Trachydoras paraguayensis (2n=56), Trachydoras steindachneri (2n=60), Ossancora punctata (2n=66) and Platydoras hancockii whose karyotypes showed a multiple NOR system. The large thorny catfishes, such as Pterodoras granulosus, Oxydoras niger and Centrodoras brachiatus share several karyotype features, with subtle variations only regarding their heterochromatin distribution. On the other hand, a remarkable karyotypic variability has been reported in the fimbriate barbells thorny catfishes. These two contrasting karyoevolution trajectories emerged from a complex interaction between chromosome rearrangements (e.g., inversions and Robertsonian translocations) and mechanisms of heterochromatin dispersion. Moreover, we believe that biological features, such as microhabitats preferences, populational size, low vagility and migratory behavior played a key role during the origin and maintenance of chromosome diversity in Doradinae subfamily.

Microchromosomes are building blocks of bird, reptile, and mammal chromosomes

Microchromosomes, once considered unimportant shreds of the chicken genome, are gene-rich elements with a high GC content and few transposable elements. Their origin has been debated for decades. We used cytological and whole-genome sequence comparisons, and chromosome conformation capture, to trace their origin and fate in genomes of reptiles, birds, and mammals. We find that microchromosomes as well as macrochromosomes are highly conserved across birds and share synteny with single small chromosomes of the chordate amphioxus, attesting to their origin as elements of an ancient animal genome. Turtles and squamates (snakes and lizards) share different subsets of ancestral microchromosomes, having independently lost microchromosomes by fusion with other microchromosomes or macrochromosomes. Patterns of fusions were quite different in different lineages. Cytological observations show that microchromosomes in all lineages are spatially separated into a central compartment at interphase and during mitosis and meiosis. This reflects higher interaction between microchromosomes than with macrochromosomes, as observed by chromosome conformation capture, and suggests some functional coherence. In highly rearranged genomes fused microchromosomes retain most ancestral characteristics, but these may erode over evolutionary time; surprisingly, de novo microchromosomes have rapidly adopted high interaction. Some chromosomes of early-branching monotreme mammals align to several bird microchromosomes, suggesting multiple microchromosome fusions in a mammalian ancestor. Subsequently, multiple rearrangements fueled the extraordinary karyotypic diversity of therian mammals. Thus, microchromosomes, far from being aberrant genetic elements, represent fundamental building blocks of amniote chromosomes, and it is mammals, rather than reptiles and birds, that are atypical.

What Have We Learned from the First 500 Avian Genomes?

The increased capacity of DNA sequencing has significantly advanced our understanding of the phylogeny of birds and the proximate and ultimate mechanisms molding their genomic diversity. In less than a decade, the number of available avian reference genomes has increased to over 500—approximately 5% of bird diversity—placing birds in a privileged position to advance the fields of phylogenomics and comparative, functional, and population genomics. Whole-genome sequence data, as well as indels and rare genomic changes, are further resolving the avian tree of life. The accumulation of bird genomes, increasingly with long-read sequence data, greatly improves the resolution of genomic features such as germline-restricted chromosomes and the W chromosome, and is facilitating the comparative integration of genotypes and phenotypes. Community-based initiatives such as the Bird 10,000 Genomes Project and Vertebrate Genome Project are playing a fundamental role in amplifying and coalescing a vibrant international program in avian comparative genomics.

Genome collinearity analysis illuminates the evolution of donkey chromosome 1 and horse chromosome 5 in perissodactyls: A comparative study

BACKGROUND

It is important to resolve the evolutionary history of species genomes as it has affected both genome organization and chromosomal architecture. The rapid innovation in sequencing technologies and the improvement in assembly algorithms have enabled the creation of highly contiguous genomes. DNA Zoo, a global organization dedicated to animal conservation, offers more than 150 chromosome-length genome assemblies. This database has great potential in the comparative genomics field.

RESULTS

Using the donkey (Equus asinus asinus, EAS) genome provided by DNA Zoo as an example, the scaffold N50 length and Benchmarking Universal Single-Copy Ortholog score reached 95.5 Mb and 91.6%, respectively. We identified the cytogenetic nomenclature, corrected the direction of the chromosome-length sequence of the donkey genome, analyzed the genome-wide chromosomal rearrangements between the donkey and horse, and illustrated the evolution of the donkey chromosome 1 and horse chromosome 5 in perissodactyls.

CONCLUSIONS

The donkey genome provided by DNA Zoo has relatively good continuity and integrity. Sequence-based comparative genomic analyses are useful for chromosome evolution research. Several previously published chromosome painting results can be used to identify the cytogenetic nomenclature and correct the direction of the chromosome-length sequence of new assemblies. Compared with the horse genome, the donkey chromosomes 1, 4, 20, and X have several obvious inversions, consistent with the results of previous studies. A 4.8 Mb inverted structure was first discovered in the donkey chromosome 25 and plains zebra chromosome 11. We speculate that the inverted structure and the tandem fusion of horse chromosome 31 and 4 are common features of non-caballine equids, which supports the correctness of the existing Equus phylogeny to an extent.

Comparative genomics provides insights into the aquatic adaptations of mammals

Divergent lineages can respond to common environmental factors through convergent processes involving shared genomic components or pathways, but the molecular mechanisms are poorly understood. Here, we provide genomic resources and insights into the evolution of mammalian lineages adapting to aquatic life. Our data suggest convergent evolution, for example, in association with thermoregulation through genes associated with a surface heat barrier (NFIA) and internal heat exchange (SEMA3E). Combined with the support of previous reports showing that the UCP1 locus has been lost in many marine mammals independently, our results suggest that the thermostatic strategy of marine mammals shifted from enhancing heat production to limiting heat loss.

The ancestors of marine mammals once roamed the land and independently committed to an aquatic lifestyle. These macroevolutionary transitions have intrigued scientists for centuries. Here, we generated high-quality genome assemblies of 17 marine mammals (11 cetaceans and six pinnipeds), including eight assemblies at the chromosome level. Incorporating previously published data, we reconstructed the marine mammal phylogeny and population histories and identified numerous idiosyncratic and convergent genomic variations that possibly contributed to the transition from land to water in marine mammal lineages. Genes associated with the formation of blubber (NFIA), vascular development (SEMA3E), and heat production by brown adipose tissue (UCP1) had unique changes that may contribute to marine mammal thermoregulation. We also observed many lineage-specific changes in the marine mammals, including genes associated with deep diving and navigation. Our study advances understanding of the timing, pattern, and molecular changes associated with the evolution of mammalian lineages adapting to aquatic life.

Genomic architecture constrained placental mammal X Chromosome evolution

Susumu Ohno proposed that the gene content of the mammalian X Chromosome should remain highly conserved due to dosage compensation. X Chromosome linkage (gene order) conservation is widespread in placental mammals but does not fall within the scope of Ohno's prediction and may be an indirect result of selection on gene content or selection against rearrangements that might disrupt X-Chromosome inactivation (XCI). Previous comparisons between the human and mouse X Chromosome sequences have suggested that although single-copy X Chromosome genes are conserved between species, most ampliconic genes were independently acquired. To better understand the evolutionary and functional constraints on X-linked gene content and linkage conservation in placental mammals, we aligned a new, high-quality, long-read X Chromosome reference assembly from the domestic cat (incorporating 19.3 Mb of targeted BAC clone sequence) to the pig, human, and mouse assemblies. A comprehensive analysis of annotated X-linked orthologs in public databases demonstrated that the majority of ampliconic gene families were present on the ancestral placental X Chromosome. We generated a domestic cat Hi-C contact map from an F1 domestic cat/Asian leopard cat hybrid and demonstrated the formation of the bipartite structure found in primate and rodent inactivated X Chromosomes. Conservation of gene order and recombination patterns is attributable to strong selective constraints on three-dimensional genomic architecture necessary for superloop formation. Species with rearranged X Chromosomes retain the ancestral order and relative spacing of loci critical for superloop formation during XCI, with compensatory inversions evolving to maintain these long-range physical interactions.

A telomere-to-telomere assembly of Oscheius tipulae and the evolution of rhabditid nematode chromosomes

Eukaryotic chromosomes have phylogenetic persistence. In many taxa, each chromosome has a single functional centromere with essential roles in spindle attachment and segregation. Fusion and fission can generate chromosomes with no or multiple centromeres, leading to genome instability. Groups with holocentric chromosomes (where centromeric function is distributed along each chromosome) might be expected to show karyotypic instability. This is generally not the case, and in Caenorhabditis elegans, it has been proposed that the role of maintenance of a stable karyotype has been transferred to the meiotic pairing centers, which are found at one end of each chromosome. Here, we explore the phylogenetic stability of nematode chromosomes using a new telomere-to-telomere assembly of the rhabditine nematode Oscheius tipulae generated from nanopore long reads. The 60-Mb O. tipulae genome is resolved into six chromosomal molecules. We find the evidence of specific chromatin diminution at all telomeres. Comparing this chromosomal O. tipulae assembly with chromosomal assemblies of diverse rhabditid nematodes, we identify seven ancestral chromosomal elements (Nigon elements) and present a model for the evolution of nematode chromosomes through rearrangement and fusion of these elements. We identify frequent fusion events involving NigonX, the element associated with the rhabditid X chromosome, and thus sex chromosome-associated gene sets differ markedly between species. Despite the karyotypic stability, gene order within chromosomes defined by Nigon elements is not conserved. Our model for nematode chromosome evolution provides a platform for investigation of the tensions between local genome rearrangement and karyotypic evolution in generating extant genome architectures.

Two-phased evolution: Genome chaos-mediated information creation and maintenance

Cancer is traditionally labeled a "cellular growth problem." However, it is fundamentally an issue of macroevolution where new systems emerge from tissue by breaking various constraints. To study this process, we used experimental platforms to "watch evolution in action" by comparing the profiles of karyotypes, transcriptomes, and cellular phenotypes longitudinally before, during, and after key phase transitions. This effort, alongside critical rethinking of current gene-based genomic and evolutionary theory, led to the development of the Genome Architecture Theory. Following a brief historical review, we present four case studies and their takeaways to describe the pattern of genome-based cancer evolution. Our discoveries include 1. The importance of non-clonal chromosome aberrations or NCCAs; 2. Two-phased cancer evolution, comprising a punctuated phase and a gradual phase, dominated by karyotype changes and gene mutation/epigenetic alterations, respectively; 3. How the karyotype codes system inheritance, which organizes gene interactions and provides the genomic basis for physiological regulatory networks; and 4. Stress-induced genome chaos, which creates genomic information by reorganizing chromosomes for macroevolution. Together, these case studies redefine the relationship between cellular macro- and microevolution: macroevolution does not equal microevolution + time. Furthermore, we incorporate genome chaos and gene mutation in a general model: genome reorganization creates new karyotype coding, then diverse cancer gene mutations can promote the dominance of tumor cell populations. Finally, we call for validation of the Genome Architecture Theory of cancer and organismal evolution, as well as the systematic study of genomic information flow in evolutionary processes.

Platypus and echidna genomes reveal mammalian biology and evolution

Egg-laying mammals (monotremes) are the only extant mammalian outgroup to therians (marsupial and eutherian animals) and provide key insights into mammalian evolution1,2. Here we generate and analyse reference genomes of the platypus (Ornithorhynchus anatinus) and echidna (Tachyglossus aculeatus), which represent the only two extant monotreme lineages. The nearly complete platypus genome assembly has anchored almost the entire genome onto chromosomes, markedly improving the genome continuity and gene annotation. Together with our echidna sequence, the genomes of the two species allow us to detect the ancestral and lineage-specific genomic changes that shape both monotreme and mammalian evolution. We provide evidence that the monotreme sex chromosome complex originated from an ancestral chromosome ring configuration. The formation of such a unique chromosome complex may have been facilitated by the unusually extensive interactions between the multi-X and multi-Y chromosomes that are shared by the autosomal homologues in humans. Further comparative genomic analyses unravel marked differences between monotremes and therians in haptoglobin genes, lactation genes and chemosensory receptor genes for smell and taste that underlie the ecological adaptation of monotremes.

Towards complete and error-free genome assemblies of all vertebrate species

High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species1-4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.

A towering genome: Experimentally validated adaptations to high blood pressure and extreme stature in the giraffe

The suite of adaptations associated with the extreme stature of the giraffe has long interested biologists and physiologists. By generating a high-quality chromosome-level giraffe genome and a comprehensive comparison with other ruminant genomes, we identified a robust catalog of giraffe-specific mutations. These are primarily related to cardiovascular, bone growth, vision, hearing, and circadian functions. Among them, the giraffe FGFRL1 gene is an outlier with seven unique amino acid substitutions not found in any other ruminant. Gene-edited mice with the giraffe-type FGFRL1 show exceptional hypertension resistance and higher bone mineral density, both of which are tightly connected with giraffe adaptations to high stature. Our results facilitate a deeper understanding of the molecular mechanism underpinning distinct giraffe traits, and may provide insights into the study of hypertension in humans.

Scalable multiple whole-genome alignment and locally collinear block construction with SibeliaZ

Multiple whole-genome alignment is a challenging problem in bioinformatics. Despite many successes, current methods are not able to keep up with the growing number, length, and complexity of assembled genomes, especially when computational resources are limited. Approaches based on compacted de Bruijn graphs to identify and extend anchors into locally collinear blocks have potential for scalability, but current methods do not scale to mammalian genomes. We present an algorithm, SibeliaZ-LCB, for identifying collinear blocks in closely related genomes based on analysis of the de Bruijn graph. We further incorporate this into a multiple whole-genome alignment pipeline called SibeliaZ. SibeliaZ shows run-time improvements over other methods while maintaining accuracy. On sixteen recently-assembled strains of mice, SibeliaZ runs in under 16 hours on a single machine, while other tools did not run to completion for eight mice within a week. SibeliaZ makes a significant step towards improving scalability of multiple whole-genome alignment and collinear block reconstruction algorithms on a single machine.

Karyotype Evolution in 10 Pinniped Species: Variability of Heterochromatin versus High Conservatism of Euchromatin as Revealed by Comparative Molecular Cytogenetics

Pinnipedia karyotype evolution was studied here using human, domestic dog, and stone marten whole-chromosome painting probes to obtain comparative chromosome maps among species of Odobenidae (Odobenus rosmarus), Phocidae (Phoca vitulina, Phoca largha, Phoca hispida, Pusa sibirica, Erignathus barbatus), and Otariidae (Eumetopias jubatus, Callorhinus ursinus, Phocarctos hookeri, and Arctocephalus forsteri). Structural and functional chromosomal features were assessed with telomere repeat and ribosomal-DNA probes and by CBG (C-bands revealed by barium hydroxide treatment followed by Giemsa staining) and CDAG (Chromomycin A3-DAPI after G-banding) methods. We demonstrated diversity of heterochromatin among pinniped karyotypes in terms of localization, size, and nucleotide composition. For the first time, an intrachromosomal rearrangement common for Otariidae and Odobenidae was revealed. We postulate that the order of evolutionarily conserved segments in the analyzed pinnipeds is the same as the order proposed for the ancestral Carnivora karyotype (2n = 38). The evolution of conserved genomes of pinnipeds has been accompanied by few fusion events (less than one rearrangement per 10 million years) and by novel intrachromosomal changes including the emergence of new centromeres and pericentric inversion/centromere repositioning. The observed interspecific diversity of pinniped karyotypes driven by constitutive heterochromatin variation likely has played an important role in karyotype evolution of pinnipeds, thereby contributing to the differences of pinnipeds’ chromosome sets.

Genome chaos: Creating new genomic information essential for cancer macroevolution

Cancer research has traditionally focused on the characterization of individual molecular mechanisms that can contribute to cancer. Due to the multiple levels of genomic and non-genomic heterogeneity, however, overwhelming molecular mechanisms have been identified, most with low clinical predictability. It is thus necessary to search for new concepts to unify these diverse mechanisms and develop better strategies to understand and treat cancer. In recent years, two-phased cancer evolution (comprised of the genome reorganization-mediated punctuated phase and gene mutation-mediated stepwise phase), initially described by tracing karyotype evolution, was confirmed by the Cancer Genome Project. In particular, genome chaos, the process of rapid and massive genome reorganization, has been commonly detected in various cancers-especially during key phase transitions, including cellular transformation, metastasis, and drug resistance-suggesting the importance of genome-level changes in cancer evolution. In this Perspective, genome chaos is used as a discussion point to illustrate new genome-mediated somatic evolutionary frameworks. By rephrasing cancer as a new system emergent from normal tissue, we present the multiple levels (or scales) of genomic and non-genomic information. Of these levels, evolutionary studies at the chromosomal level are determined to be of ultimate importance, since altered genomes change the karyotype coding and karyotype change is the key event for punctuated cellular macroevolution. Using this lens, we differentiate and analyze developmental processes and cancer evolution, as well as compare the informational relationship between genome chaos and its various subtypes in the context of macroevolution under crisis. Furthermore, the process of deterministic genome chaos is discussed to interpret apparently random events (including stressors, chromosomal variation subtypes, surviving cells with new karyotypes, and emergent stable cellular populations) as nonrandom patterns, which supports the new cancer evolutionary model that unifies genome and gene contributions during different phases of cancer evolution. Finally, the new perspective of using cancer as a model for organismal evolution is briefly addressed, emphasizing the Genome Theory as a new and necessary conceptual framework for future research and its practical implications, not only in cancer but evolutionary biology as a whole.

Vertebrate Chromosome Evolution

The study of chromosome evolution is undergoing a resurgence of interest owing to advances in DNA sequencing technology that facilitate the production of chromosome-scale whole-genome assemblies de novo. This review focuses on the history, methods, discoveries, and current challenges facing the field, with an emphasis on vertebrate genomes. A detailed examination of the literature on the biology of chromosome rearrangements is presented, specifically the relationship between chromosome rearrangements and phenotypic evolution, adaptation, and speciation. A critical review of the methods for identifying, characterizing, and visualizing chromosome rearrangements and computationally reconstructing ancestral karyotypes is presented. We conclude by looking to the future, identifying the enormous technical and scientific challenges presented by the accumulation of hundreds and eventually thousands of chromosome-scale assemblies.

A comparative genomics multitool for scientific discovery and conservation

The Zoonomia Project is investigating the genomics of shared and specialized traits in eutherian mammals. Here we provide genome assemblies for 131 species, of which all but 9 are previously uncharacterized, and describe a whole-genome alignment of 240 species of considerable phylogenetic diversity, comprising representatives from more than 80% of mammalian families. We find that regions of reduced genetic diversity are more abundant in species at a high risk of extinction, discern signals of evolutionary selection at high resolution and provide insights from individual reference genomes. By prioritizing phylogenetic diversity and making data available quickly and without restriction, the Zoonomia Project aims to support biological discovery, medical research and the conservation of biodiversity.

A whole-genome alignment of 240 phylogenetically diverse species of eutherian mammal—including 131 previously uncharacterized species—from the Zoonomia Project provides data that support biological discovery, medical research and conservation.

Marsupial chromosomics: bridging the gap between genomes and chromosomes

Marsupials have unique features that make them particularly interesting to study, and sequencing of marsupial genomes is helping to understand their evolution. A decade ago, it was a huge feat to sequence the first marsupial genome. Now, the advances in sequencing technology have made the sequencing of many more marsupial genomes possible. However, the DNA sequence is only one component of the structures it is packaged into: chromosomes. Knowing the arrangement of the DNA sequence on each chromosome is essential for a genome assembly to be used to its full potential. The importance of combining sequence information with cytogenetics has previously been demonstrated for rapidly evolving regions of the genome, such as the sex chromosomes, as well as for reconstructing the ancestral marsupial karyotype and understanding the chromosome rearrangements involved in the Tasmanian devil facial tumour disease. Despite the recent advances in sequencing technology assisting in genome assembly, physical anchoring of the sequence to chromosomes is required to achieve a chromosome-level assembly. Once chromosome-level assemblies are achieved for more marsupials, we will be able to investigate changes in the packaging and interactions between chromosomes to gain an understanding of the role genome architecture has played during marsupial evolution.

An Indo-Pacific Humpback Dolphin Genome Reveals Insights into Chromosome Evolution and the Demography of a Vulnerable Species

The Indo-Pacific humpback dolphin (Sousa chinensis) is a small inshore species of odontocete cetacean listed as Vulnerable on the IUCN Red List. Here, we report on the evolution of S. chinensis chromosomes from its cetruminant ancestor and elucidate the evolutionary history and population genetics of two neighboring S. chinensis populations. We found that breakpoints in ancestral chromosomes leading to S. chinensis could have affected the function of genes related to kidney filtration, body development, and immunity. Resequencing of individuals from two neighboring populations in the northwestern South China Sea, Leizhou Bay and Sanniang Bay, revealed genetic differentiation, low diversity, and small contemporary effective population sizes. Demographic analyses showed a marked decrease in the population size of the two investigated populations over the last ~4,000 years, possibly related to climatic oscillations. This study implies a high risk of extinction and strong conservation requirement for the Indo-Pacific humpback dolphin.

•

Deducing chromosome evolution from ancestral Cetruminantia and ancestral Odontoceti

•

Reconstructing the demographic history of Sousa chinensis

•

Implying high risk of extinction and strong conservation requirement for S. chinensis

Evolutionary rearrangements of X chromosomes in voles (Arvicolinae, Rodentia)

Euchromatic segments of the X chromosomes of placental mammals are the most conservative elements of the karyotype, only rarely subjected to either inter- or intrachromosomal rearrangements. Here, using microdissection-derived set of region-specific probes of Terricola savii we detailed the evolutionary rearrangements found in X chromosomes in 20 vole species (Arvicolinae, Rodentia). We show that the evolution of X chromosomes in this taxon was accompanied by multiple para- and pericentric inversions and centromere shifts. The contribution of intrachromosomal rearrangements to the karyotype evolution of Arvicolinae species was approximately equivalent in both the separate autosomal conserved segments and the X chromosomes. Intrachromosmal rearrangements and structural reorganization of the X chromosomes was likely accompanied by an accumulation, distribution, and evolution of repeated sequences.

Evolution of the Human Chromosome 13 Synteny: Evolutionary Rearrangements, Plasticity, Human Disease Genes and Cancer Breakpoints

The history of each human chromosome can be studied through comparative cytogenetic approaches in mammals which permit the identification of human chromosomal homologies and rearrangements between species. Comparative banding, chromosome painting, Bacterial Artificial Chromosome (BAC) mapping and genome data permit researchers to formulate hypotheses about ancestral chromosome forms. Human chromosome 13 has been previously shown to be conserved as a single syntenic element in the Ancestral Primate Karyotype; in this context, in order to study and verify the conservation of primate chromosomes homologous to human chromosome 13, we mapped a selected set of BAC probes in three platyrrhine species, characterised by a high level of rearrangements, using fluorescence in situ hybridisation (FISH). Our mapping data on Saguinus oedipus, Callithrix argentata and Alouatta belzebul provide insight into synteny of human chromosome 13 evolution in a comparative perspective among primate species, showing rearrangements across taxa. Furthermore, in a wider perspective, we have revised previous cytogenomic literature data on chromosome 13 evolution in eutherian mammals, showing a complex origin of the eutherian mammal ancestral karyotype which has still not been completely clarified. Moreover, we analysed biomedical aspects (the OMIM and Mitelman databases) regarding human chromosome 13, showing that this autosome is characterised by a certain level of plasticity that has been implicated in many human cancers and diseases.

Evolution of Marsupial Genomes

Marsupial genomes, which are packaged into large chromosomes, provide a powerful resource for studying the mechanisms of genome evolution. The extensive and valuable body of work on marsupial cytogenetics, combined more recently with genome sequence data, has enabled prediction of the 2n = 14 karyotype ancestral to all marsupial families. The application of both chromosome biology and genome sequencing, or chromosomics, has been a necessary approach for various aspects of mammalian genome evolution, such as understanding sex chromosome evolution and the origin and evolution of transmissible tumors in Tasmanian devils. The next phase of marsupial genome evolution research will employ chromosomics approaches to begin addressing fundamental questions in marsupial genome evolution and chromosome evolution more generally. The answers to these complex questions will impact our understanding across a broad range of fields, including the genetics of speciation, genome adaptation to environmental stressors, and species management.

Chromosome-level genome assembly for giant panda provides novel insights into Carnivora chromosome evolution

BACKGROUND

Chromosome evolution is an important driver of speciation and species evolution. Previous studies have detected chromosome rearrangement events among different Carnivora species using chromosome painting strategies. However, few of these studies have focused on chromosome evolution at a nucleotide resolution due to the limited availability of chromosome-level Carnivora genomes. Although the de novo genome assembly of the giant panda is available, current short read-based assemblies are limited to moderately sized scaffolds, making the study of chromosome evolution difficult.

RESULTS

Here, we present a chromosome-level giant panda draft genome with a total size of 2.29 Gb. Based on the giant panda genome and published chromosome-level dog and cat genomes, we conduct six large-scale pairwise synteny alignments and identify evolutionary breakpoint regions. Interestingly, gene functional enrichment analysis shows that for all of the three Carnivora genomes, some genes located in evolutionary breakpoint regions are significantly enriched in pathways or terms related to sensory perception of smell. In addition, we find that the sweet receptor gene TAS1R2, which has been proven to be a pseudogene in the cat genome, is located in an evolutionary breakpoint region of the giant panda, suggesting that interchromosomal rearrangement may play a role in the cat TAS1R2 pseudogenization.

CONCLUSIONS

We show that the combined strategies employed in this study can be used to generate efficient chromosome-level genome assemblies. Moreover, our comparative genomics analyses provide novel insights into Carnivora chromosome evolution, linking chromosome evolution to functional gene evolution.