Semi-automated assembly of high-quality diploid human reference genomes.

The current human reference genome, GRCh38, represents over 20 years of effort to generate a high-quality assembly, which has benefitted society1,2. However, it still has many gaps and errors, and does not represent a biological genome as it is a blend of multiple individuals3,4. Recently, a high-quality telomere-to-telomere reference, CHM13, was generated with the latest long-read technologies, but it was derived from a hydatidiform mole cell line with a nearly homozygous genome5. To address these limitations, the Human Pangenome Reference Consortium formed with the goal of creating high-quality, cost-effective, diploid genome assemblies for a pangenome reference that represents human genetic diversity6. Here, in our first scientific report, we determined which combination of current genome sequencing and assembly approaches yield the most complete and accurate diploid genome assembly with minimal manual curation. Approaches that used highly accurate long reads and parent-child data with graph-based haplotype phasing during assembly outperformed those that did not. Developing a combination of the top-performing methods, we generated our first high-quality diploid reference assembly, containing only approximately four gaps per chromosome on average, with most chromosomes within ±1% of the length of CHM13. Nearly 48% of protein-coding genes have non-synonymous amino acid changes between haplotypes, and centromeric regions showed the highest diversity. Our findings serve as a foundation for assembling near-complete diploid human genomes at scale for a pangenome reference to capture global genetic variation from single nucleotides to structural rearrangements.

Evolutionary and biomedical insights from a marmoset diploid genome assembly.

The accurate and complete assembly of both haplotype sequences of a diploid organism is essential to understanding the role of variation in genome functions, phenotypes and diseases1. Here, using a trio-binning approach, we present a high-quality, diploid reference genome, with both haplotypes assembled independently at the chromosome level, for the common marmoset (Callithrix jacchus), an primate model system that is widely used in biomedical research2,3. The full spectrum of heterozygosity between the two haplotypes involves 1.36% of the genome-much higher than the 0.13% indicated by the standard estimation based on single-nucleotide heterozygosity alone. The de novo mutation rate is 0.43 × 10-8 per site per generation, and the paternal inherited genome acquired twice as many mutations as the maternal. Our diploid assembly enabled us to discover a recent expansion of the sex-differentiation region and unique evolutionary changes in the marmoset Y chromosome. In addition, we identified many genes with signatures of positive selection that might have contributed to the evolution of Callithrix biological features. Brain-related genes were highly conserved between marmosets and humans, although several genes experienced lineage-specific copy number variations or diversifying selection, with implications for the use of marmosets as a model system.

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.

The genome of Litomosoides sigmodontis illuminates the origins of Y chromosomes in filarial nematodes.

Heteromorphic sex chromosomes are usually thought to have originated from a pair of autosomes that acquired a sex-determining locus and subsequently stopped recombining, leading to degeneration of the sex-limited chromosome. The majority of nematode species lack heteromorphic sex chromosomes and determine sex using an X-chromosome counting mechanism, with males being hemizygous for one or more X chromosomes (XX/X0). Some filarial nematode species, including important parasites of humans, have heteromorphic XX/XY karyotypes. It has been assumed that sex is determined by a Y-linked locus in these species. However, karyotypic analyses suggested that filarial Y chromosomes are derived from the unfused homologue of an autosome involved in an X-autosome fusion event. Here, we generated a chromosome-level reference genome for Litomosoides sigmodontis, a filarial nematode with the ancestral filarial karyotype and sex determination mechanism (XX/X0). By mapping the assembled chromosomes to the rhabditid nematode ancestral linkage (or Nigon) elements, we infer that the ancestral filarial X chromosome was the product of a fusion between NigonX (the ancestrally X-linked element) and NigonD (ancestrally autosomal). In the two filarial lineages with XY systems, there have been two independent X-autosome chromosome fusion events involving different autosomal Nigon elements. In both lineages, the region shared by the neo-X and neo-Y chromosomes is within the ancestrally autosomal portion of the X, confirming that the filarial Y chromosomes are derived from the unfused homologue of the autosome. Sex determination in XY filarial nematodes therefore likely continues to operate via the ancestral X-chromosome counting mechanism, rather than via a Y-linked sex-determining locus.

Telomere-to-telomere assembly of a complete human X chromosome.

After two decades of improvements, the current human reference genome (GRCh38) is the most accurate and complete vertebrate genome ever produced. However, no single chromosome has been finished end to end, and hundreds of unresolved gaps persist1,2. Here we present a human genome assembly that surpasses the continuity of GRCh382, along with a gapless, telomere-to-telomere assembly of a human chromosome. This was enabled by high-coverage, ultra-long-read nanopore sequencing of the complete hydatidiform mole CHM13 genome, combined with complementary technologies for quality improvement and validation. Focusing our efforts on the human X chromosome3, we reconstructed the centromeric satellite DNA array (approximately 3.1 Mb) and closed the 29 remaining gaps in the current reference, including new sequences from the human pseudoautosomal regions and from cancer-testis ampliconic gene families (CT-X and GAGE). These sequences will be integrated into future human reference genome releases. In addition, the complete chromosome X, combined with the ultra-long nanopore data, allowed us to map methylation patterns across complex tandem repeats and satellite arrays. Our results demonstrate that finishing the entire human genome is now within reach, and the data presented here will facilitate ongoing efforts to complete the other human chromosomes.

Divergent sensory and immune gene evolution in sea turtles with contrasting demographic and life histories.

Sea turtles represent an ancient lineage of marine vertebrates that evolved from terrestrial ancestors over 100 Mya. The genomic basis of the unique physiological and ecological traits enabling these species to thrive in diverse marine habitats remains largely unknown. Additionally, many populations have drastically declined due to anthropogenic activities over the past two centuries, and their recovery is a high global conservation priority. We generated and analyzed high-quality reference genomes for the leatherback (Dermochelys coriacea) and green (Chelonia mydas) turtles, representing the two extant sea turtle families. These genomes are highly syntenic and homologous, but localized regions of noncollinearity were associated with higher copy numbers of immune, zinc-finger, and olfactory receptor (OR) genes in green turtles, with ORs related to waterborne odorants greatly expanded in green turtles. Our findings suggest that divergent evolution of these key gene families may underlie immunological and sensory adaptations assisting navigation, occupancy of neritic versus pelagic environments, and diet specialization. Reduced collinearity was especially prevalent in microchromosomes, with greater gene content, heterozygosity, and genetic distances between species, supporting their critical role in vertebrate evolutionary adaptation. Finally, diversity and demographic histories starkly contrasted between species, indicating that leatherback turtles have had a low yet stable effective population size, exhibit extremely low diversity compared with other reptiles, and harbor a higher genetic load compared with green turtles, reinforcing concern over their persistence under future climate scenarios. These genomes provide invaluable resources for advancing our understanding of evolution and conservation best practices in an imperiled vertebrate lineage.

A High-Quality Blue Whale Genome, Segmental Duplications, and Historical Demography.

The blue whale, Balaenoptera musculus, is the largest animal known to have ever existed, making it an important case study in longevity and resistance to cancer. To further this and other blue whale-related research, we report a reference-quality, long-read-based genome assembly of this fascinating species. We assembled the genome from PacBio long reads and utilized Illumina/10×, optical maps, and Hi-C data for scaffolding, polishing, and manual curation. We also provided long read RNA-seq data to facilitate the annotation of the assembly by NCBI and Ensembl. Additionally, we annotated both haplotypes using TOGA and measured the genome size by flow cytometry. We then compared the blue whale genome with other cetaceans and artiodactyls, including vaquita (Phocoena sinus), the world's smallest cetacean, to investigate blue whale's unique biological traits. We found a dramatic amplification of several genes in the blue whale genome resulting from a recent burst in segmental duplications, though the possible connection between this amplification and giant body size requires further study. We also discovered sites in the insulin-like growth factor-1 gene correlated with body size in cetaceans. Finally, using our assembly to examine the heterozygosity and historical demography of Pacific and Atlantic blue whale populations, we found that the genomes of both populations are highly heterozygous and that their genetic isolation dates to the last interglacial period. Taken together, these results indicate how a high-quality, annotated blue whale genome will serve as an important resource for biology, evolution, and conservation research.

Chasing perfection: validation and polishing strategies for telomere-to-telomere genome assemblies.

Advances in long-read sequencing technologies and genome assembly methods have enabled the recent completion of the first telomere-to-telomere human genome assembly, which resolves complex segmental duplications and large tandem repeats, including centromeric satellite arrays in a complete hydatidiform mole (CHM13). Although derived from highly accurate sequences, evaluation revealed evidence of small errors and structural misassemblies in the initial draft assembly. To correct these errors, we designed a new repeat-aware polishing strategy that made accurate assembly corrections in large repeats without overcorrection, ultimately fixing 51% of the existing errors and improving the assembly quality value from 70.2 to 73.9 measured from PacBio high-fidelity and Illumina k-mers. By comparing our results to standard automated polishing tools, we outline common polishing errors and offer practical suggestions for genome projects with limited resources. We also show how sequencing biases in both high-fidelity and Oxford Nanopore Technologies reads cause signature assembly errors that can be corrected with a diverse panel of sequencing technologies.

Caecilian Genomes Reveal the Molecular Basis of Adaptation and Convergent Evolution of Limblessness in Snakes and Caecilians.

We present genome sequences for the caecilians Geotrypetes seraphini (3.8 Gb) and Microcaecilia unicolor (4.7 Gb), representatives of a limbless, mostly soil-dwelling amphibian clade with reduced eyes, and unique putatively chemosensory tentacles. More than 69% of both genomes are composed of repeats, with retrotransposons being the most abundant. We identify 1,150 orthogroups that are unique to caecilians and enriched for functions in olfaction and detection of chemical signals. There are 379 orthogroups with signatures of positive selection on caecilian lineages with roles in organ development and morphogenesis, sensory perception, and immunity amongst others. We discover that caecilian genomes are missing the zone of polarizing activity regulatorysequence (ZRS) enhancer of Sonic Hedgehog which is also mutated in snakes. In vivo deletions have shown ZRS is required for limb development in mice, thus, revealing a shared molecular target implicated in the independent evolution of limblessness in snakes and caecilians.

Do brushings taken during bronchoscopy for suspected primary lung cancer increase diagnostic yield? An audit of real-world clinical practice.

We reviewed clinical records to determine whether the use of bronchial brushings improved diagnostic yield in a setting where bronchoscopy for suspected primary lung cancer is routinely guided by prior chest computed tomography but endobronchial ultrasound-guided sampling is unavailable. For 29% of cases who had brushings and at least one other test taken (bronchial biopsies or washings), the histological diagnosis was made solely on the basis of samples obtained by brushings.

Surfaces Containing Sharp Nanostructures Enhance Antibiotic Efficacy.

The ever-increasing rate of medical device implantations is met by a proportionately high burden of implant-associated infections. To mitigate this threat, much research has been directed toward the development of antibacterial surface modifications by various means. One recent approach involves surfaces containing sharp nanostructures capable of killing bacteria upon contact. Herein, we report that the mechanical interaction between Staphylococcus aureus and such surface nanostructures leads to a sensitization of the pathogen to the glycopeptide antibiotic vancomycin. We demonstrate that this is due to cell wall damage and impeded bacterial defenses against reactive oxygen species. The results of this study promise to be impactful in the clinic, as a combination of nanostructured antibacterial surfaces and antibiotics commonly used in hospitals may improve antimicrobial therapy strategies, helping clinicians to prevent and treat implant-associated infections using reduced antibiotic concentrations instead of relying on invasive revision surgeries with often poor outcomes.

Prioritizing Endangered Species in Genome Sequencing: Conservation Genomics in Action with the First Platinum-Standard Reference-Quality Genome of the Critically Endangered European Mink Mustela lutreola L., 1761.

The European mink Mustela lutreola (Mustelidae) ranks among the most endangered mammalian species globally, experiencing a rapid and severe decline in population size, density, and distribution. Given the critical need for effective conservation strategies, understanding its genomic characteristics becomes paramount. To address this challenge, the platinum-quality, chromosome-level reference genome assembly for the European mink was successfully generated under the project of the European Mink Centre consortium. Leveraging PacBio HiFi long reads, we obtained a 2586.3 Mbp genome comprising 25 scaffolds, with an N50 length of 154.1 Mbp. Through Hi-C data, we clustered and ordered the majority of the assembly (>99.9%) into 20 chromosomal pseudomolecules, including heterosomes, ranging from 6.8 to 290.1 Mbp. The newly sequenced genome displays a GC base content of 41.9%. Additionally, we successfully assembled the complete mitochondrial genome, spanning 16.6 kbp in length. The assembly achieved a BUSCO (Benchmarking Universal Single-Copy Orthologs) completeness score of 98.2%. This high-quality reference genome serves as a valuable genomic resource for future population genomics studies concerning the European mink and related taxa. Furthermore, the newly assembled genome holds significant potential in addressing key conservation challenges faced by M. lutreola. Its applications encompass potential revision of management units, assessment of captive breeding impacts, resolution of phylogeographic questions, and facilitation of monitoring and evaluating the efficiency and effectiveness of dedicated conservation strategies for the European mink. This species serves as an example that highlights the paramount importance of prioritizing endangered species in genome sequencing projects due to the race against time, which necessitates the comprehensive exploration and characterization of their genomic resources before their populations face extinction.

Identifying and removing haplotypic duplication in primary genome assemblies.

MOTIVATION: Rapid development in long-read sequencing and scaffolding technologies is accelerating the production of reference-quality assemblies for large eukaryotic genomes. However, haplotype divergence in regions of high heterozygosity often results in assemblers creating two copies rather than one copy of a region, leading to breaks in contiguity and compromising downstream steps such as gene annotation. Several tools have been developed to resolve this problem. However, they either focus only on removing contained duplicate regions, also known as haplotigs, or fail to use all the relevant information and hence make errors. RESULTS: Here we present a novel tool, purge_dups, that uses sequence similarity and read depth to automatically identify and remove both haplotigs and heterozygous overlaps. In comparison with current tools, we demonstrate that purge_dups can reduce heterozygous duplication and increase assembly continuity while maintaining completeness of the primary assembly. Moreover, purge_dups is fully automatic and can easily be integrated into assembly pipelines. AVAILABILITY AND IMPLEMENTATION: The source code is written in C and is available at https://github.com/dfguan/purge_dups. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Oligodendrocyte pathology exceeds axonal pathology in white matter in human amyotrophic lateral sclerosis.

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease. The majority of cases are sporadic (sALS), while the most common inherited form is due to C9orf72 mutation (C9ALS). A high burden of inclusion pathology is seen in glia (including oligodendrocytes) in ALS, especially in C9ALS. Myelin basic protein (MBP) messenger RNA (mRNA) must be transported to oligodendrocyte processes for myelination, a possible vulnerability for normal function. TDP43 is found in pathological inclusions in ALS and is a component of mRNA transport granules. Thus, TDP43 aggregation could lead to MBP loss. Additionally, the hexanucleotide expansion of mutant C9ALS binds hnRNPA2/B1, a protein essential for mRNA transport, causing potential further impairment of hnRNPA2/B1 function, and thus myelination. Using immunohistochemistry for p62 and TDP43 in human post-mortem tissue, we found a high burden of glial inclusions in the prefrontal cortex, precentral gyrus, and spinal cord in ALS, which was greater in C9ALS than in sALS cases. Double staining demonstrated that the majority of these inclusions were in oligodendrocytes. Using immunoblotting, we demonstrated reduced MBP protein levels relative to PLP (a myelin component that relies on protein not mRNA transport) and neurofilament protein (an axonal marker) in the spinal cord. This MBP loss was disproportionate to the level of PLP and axonal loss, suggesting that impaired mRNA transport may be partly responsible. Finally, we show that in C9ALS cases, the level of oligodendroglial inclusions correlates inversely with levels of hnRNPA2/B1 and the number of oligodendrocyte precursor cells. We conclude that there is considerable oligodendrocyte pathology in ALS, which at least partially reflects impairment of mRNA transport. © 2020 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Perspective: child abuse in the military.

Child maltreatment is an unfortunate aspect of our society, afflicting civilian and military families alike. However, unlike their civilian counterparts, military families bear additional burdens inherent to military service that can exacerbate some of the root causes of child abuse. For this reason, the U.S. Department of Defense is committed to ensuring not only a highly disciplined and ready force, but also a healthy force - the foundation of which is healthy families. Therefore, understanding the military health care system, how it functions and how it collects data is a necessary first step in evaluating the efficacy of current programs and identifying opportunities for improvement. Moving beyond treatment and prevention, the military also boasts an independent judicial system designed to promote the dual interests of justice and good order as well as discipline in the armed forces, and this also contributes to a distinct culture. These two independent systems, often viewed as having diametrically opposed interests, can work together synergistically to promote the ultimate goal of fewer instances of child maltreatment in the military.

High-resolution µCT of a mouse embryo using a compact laser-driven X-ray betatron source.

In the field of X-ray microcomputed tomography (µCT) there is a growing need to reduce acquisition times at high spatial resolution (approximate micrometers) to facilitate in vivo and high-throughput operations. The state of the art represented by synchrotron light sources is not practical for certain applications, and therefore the development of high-brightness laboratory-scale sources is crucial. We present here imaging of a fixed embryonic mouse sample using a compact laser-plasma-based X-ray light source and compare the results to images obtained using a commercial X-ray µCT scanner. The radiation is generated by the betatron motion of electrons inside a dilute and transient plasma, which circumvents the flux limitations imposed by the solid or liquid anodes used in conventional electron-impact X-ray tubes. This X-ray source is pulsed (duration <30 fs), bright (>1010 photons per pulse), small (diameter <1 µm), and has a critical energy >15 keV. Stable X-ray performance enabled tomographic imaging of equivalent quality to that of the µCT scanner, an important confirmation of the suitability of the laser-driven source for applications. The X-ray flux achievable with this approach scales with the laser repetition rate without compromising the source size, which will allow the recording of high-resolution µCT scans in minutes.

Use-dependent plasticity explains aftereffects in visually guided locomotor learning of a novel step length asymmetry.

Studies of upper extremity reaching show that use-dependent plasticity, or learning from repetition, plays an important role in shaping motor behaviors. Yet the impact of repetition on locomotor learning is unclear, despite the fact that gait is developed and practiced over millions of repetitions. To test whether repetition alone can induce storage of a novel walking pattern, we instructed two groups of young healthy subjects to learn an asymmetric walking pattern through two distinct learning paradigms. The first group learned a new pattern through an established visual distortion paradigm, which provided both sensory prediction error and repetition of movement patterns to induce walking aftereffects, and the second received veridical feedback with a target change, which provided only repetition (use-dependent plasticity) to induce aftereffects. When feedback was removed, both groups demonstrated aftereffects in the primary outcome, step asymmetry index. Surprisingly, despite the different task demands, both groups produced similar aftereffect magnitudes, which also had similar rates of decay, suggesting that the addition of sensory prediction errors did not improve storage of learning beyond that induced by the use-dependent process alone. To further characterize the use-dependent process, we conducted a second experiment to quantify aftereffect size in a third group who practiced double the asymmetry magnitude. This new group showed a proportionately greater magnitude of the use-dependent aftereffect. Together, these findings show that the primary driver of storage of a new step length asymmetry during visually guided locomotor learning is repetition, not sensory prediction error, and this effect scales with the learning magnitude.NEW & NOTEWORTHY Use-dependent plasticity, or learning from repetition, is an important process for upper extremity reaching tasks, but its contribution to walking is not well established. Here, we demonstrate the existence of a dose-dependent, use-dependent process during visually guided treadmill walking. We also show that sensory prediction errors, previously thought to drive aftereffects in similar locomotor learning paradigms, do not appear to play a significant role in visually driven learning of a novel step asymmetry during treadmill walking.

Evaluation of GRCh38 and de novo haploid genome assemblies demonstrates the enduring quality of the reference assembly.

The human reference genome assembly plays a central role in nearly all aspects of today's basic and clinical research. GRCh38 is the first coordinate-changing assembly update since 2009; it reflects the resolution of roughly 1000 issues and encompasses modifications ranging from thousands of single base changes to megabase-scale path reorganizations, gap closures, and localization of previously orphaned sequences. We developed a new approach to sequence generation for targeted base updates and used data from new genome mapping technologies and single haplotype resources to identify and resolve larger assembly issues. For the first time, the reference assembly contains sequence-based representations for the centromeres. We also expanded the number of alternate loci to create a reference that provides a more robust representation of human population variation. We demonstrate that the updates render the reference an improved annotation substrate, alter read alignments in unchanged regions, and impact variant interpretation at clinically relevant loci. We additionally evaluated a collection of new de novo long-read haploid assemblies and conclude that although the new assemblies compare favorably to the reference with respect to continuity, error rate, and gene completeness, the reference still provides the best representation for complex genomic regions and coding sequences. We assert that the collected updates in GRCh38 make the newer assembly a more robust substrate for comprehensive analyses that will promote our understanding of human biology and advance our efforts to improve health.

Disrupted-in-Schizophrenia-1 is essential for normal hypothalamic-pituitary-interrenal (HPI) axis function.

Psychiatric disorders arise due to an interplay of genetic and environmental factors, including stress. Studies in rodents have shown that mutants for Disrupted-In-Schizophrenia-1 (DISC1), a well-accepted genetic risk factor for mental illness, display abnormal behaviours in response to stress, but the mechanisms through which DISC1 affects stress responses remain poorly understood. Using two lines of zebrafish homozygous mutant for disc1, we investigated behaviour and functioning of the hypothalamic-pituitary-interrenal (HPI) axis, the fish equivalent of the hypothalamic-pituitary-adrenal (HPA) axis. Here, we show that the role of DISC1 in stress responses is evolutionarily conserved and that DISC1 is essential for normal functioning of the HPI axis. Adult zebrafish homozygous mutant for disc1 show aberrant behavioural responses to stress. Our studies reveal that in the embryo, disc1 is expressed in neural progenitor cells of the hypothalamus, a conserved region of the vertebrate brain that centrally controls responses to environmental stressors. In disc1 mutant embryos, proliferating rx3+ hypothalamic progenitors are not maintained normally and neuronal differentiation is compromised: rx3-derived ff1b+ neurons, implicated in anxiety-related behaviours, and corticotrophin releasing hormone (crh) neurons, key regulators of the stress axis, develop abnormally, and rx3-derived pomc+ neurons are disorganised. Abnormal hypothalamic development is associated with dysfunctional behavioural and neuroendocrine stress responses. In contrast to wild type siblings, disc1 mutant larvae show altered crh levels, fail to upregulate cortisol levels when under stress and do not modulate shoal cohesion, indicative of abnormal social behaviour. These data indicate that disc1 is essential for normal development of the hypothalamus and for the correct functioning of the HPA/HPI axis.

The pig X and Y Chromosomes: structure, sequence, and evolution.

We have generated an improved assembly and gene annotation of the pig X Chromosome, and a first draft assembly of the pig Y Chromosome, by sequencing BAC and fosmid clones from Duroc animals and incorporating information from optical mapping and fiber-FISH. The X Chromosome carries 1033 annotated genes, 690 of which are protein coding. Gene order closely matches that found in primates (including humans) and carnivores (including cats and dogs), which is inferred to be ancestral. Nevertheless, several protein-coding genes present on the human X Chromosome were absent from the pig, and 38 pig-specific X-chromosomal genes were annotated, 22 of which were olfactory receptors. The pig Y-specific Chromosome sequence generated here comprises 30 megabases (Mb). A 15-Mb subset of this sequence was assembled, revealing two clusters of male-specific low copy number genes, separated by an ampliconic region including the HSFY gene family, which together make up most of the short arm. Both clusters contain palindromes with high sequence identity, presumably maintained by gene conversion. Many of the ancestral X-related genes previously reported in at least one mammalian Y Chromosome are represented either as active genes or partial sequences. This sequencing project has allowed us to identify genes--both single copy and amplified--on the pig Y Chromosome, to compare the pig X and Y Chromosomes for homologous sequences, and thereby to reveal mechanisms underlying pig X and Y Chromosome evolution.