Comparative analysis of mammalian Y chromosomes illuminates ancestral structure and lineage-specific evolution

The Red Fox Y-Chromosome in Comparative Context

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

Lysine demethylases KDM6A and UTY: The X and Y of histone demethylation

Histone demethylases remove transcriptional repressive marks from histones in the nucleus. KDM6A (also known as UTX) is a lysine demethylase which acts on the trimethylated lysine at position 27 in histone 3. The KDM6A gene is located on the X chromosome but escapes X inactivation even though it is not located in the pseudoautosomal region. There is a homologue of KDM6A on the Y chromosome, known as UTY. UTY was thought to have lost its demethylase activity and to represent a non-functional remnant of the ancestral KDM6A gene. However, results with knockout mice suggest that the gene is expressed and the protein performs some function within the cell. Female mice with homozygous deletion of Kdm6a do not survive, but hemizygous males are viable, attributed to the presence of the Uty gene. KDM6A is mutated in the human condition Kabuki syndrome type 2 (OMIM 300867) and in many cases of cancer. The amino acid sequence of KDM6A has been conserved across animal phyla, although it is only found on the X chromosome in eutherian mammals. In this review, we reanalyse existing data from various sources (protein sequence comparison, evolutionary genetics, transcription factor binding and gene expression analysis) to determine the function, expression and evolution of KDM6A and UTY and show that UTY has a functional role similar to KDM6A in metabolism and development.

Novel Y Chromosome Retrocopies in Canids Revealed through a Genome-Wide Association Study for Sex

The lack of an annotated reference sequence for the canine Y chromosome has limited evolutionary studies, as well as our understanding of the role of Y-linked sequences in phenotypes with a sex bias. In genome-wide association studies (GWASs), we observed spurious associations with autosomal SNPs when sex was unbalanced in case-control cohorts and hypothesized that a subset of SNPs mapped to autosomes are in fact sex-linked. Using the Illumina 230K CanineHD array in a GWAS for sex, we identified SNPs that amplify in both sexes but possess significant allele frequency differences between males and females. We found 48 SNPs mapping to 14 regions of eight autosomes and the X chromosome that are Y-linked, appearing heterozygous in males and monomorphic in females. Within these 14 regions are eight genes: three autosomal and five X-linked. We investigated the autosomal genes (MITF, PPP2CB, and WNK1) and determined that the SNPs are diverged nucleotides in retrocopies that have transposed to the Y chromosome. MITFY and WNK1Y are expressed and appeared recently in the Canidae lineage, whereas PPP2CBY represents a much older insertion with no evidence of expression in the dog. This work reveals novel canid Y chromosome sequences and provides evidence for gene transposition to the Y from autosomes and the X.

The evolutionary history of grey wolf Y chromosomes

Analyses of Y chromosome haplotypes uniquely provide a paternal picture of evolutionary histories and offer a very useful contrast to studies based on maternally inherited mitochondrial DNA (mtDNA). Here we used a bioinformatic approach based on comparison of male and female sequence coverage to identify 4.7 Mb from the grey wolf (Canis lupis) Y chromosome, probably representing most of the male-specific, nonampliconic sequence from the euchromatic part of the chromosome. We characterized this sequence and then identified ≈1,500 Y-linked single nucleotide polymorphisms in a sample of 145 resequenced male wolves, including 75 Finnish wolf genomes newly sequenced in this study, and in 24 dogs and eight other canids. We found 53 Y chromosome haplotypes, of which 26 were seen in grey wolves, that clustered in four major haplogroups. All four haplogroups were represented in samples of Finnish wolves, showing that haplogroup lineages were not partitioned on a continental scale. However, regional population structure was indicated because individual haplotypes were never shared between geographically distant areas, and genetically similar haplotypes were only found within the same geographical region. The deepest split between grey wolf haplogroups was estimated to have occurred 125,000 years ago, which is considerably older than recent estimates of the time of divergence of wolf populations. The distribution of dogs in a phylogenetic tree of Y chromosome haplotypes supports multiple domestication events, or wolf paternal introgression, starting 29,000 years ago. We also addressed the disputed origin of a recently founded population of Scandinavian wolves and observed that founding as well as most recent immigrant haplotypes were present in the neighbouring Finnish population, but not in sequenced wolves from elsewhere in the world, or in dogs.

Evolution and meiotic organization of heteromorphic sex chromosomes

The evolution of heteromorphic sex chromosomes has occurred independently many times in different lineages. The differentiation of sex chromosomes leads to dramatic changes in sequence composition and function and guides the evolutionary trajectory and utilization of genes in pivotal sex determination and reproduction roles. In addition, meiotic recombination and pairing mechanisms are key in orchestrating the resultant impact, retention and maintenance of heteromorphic sex chromosomes, as the resulting exposure of unpaired DNA at meiosis triggers ancient repair and checkpoint pathways. In this review, we summarize the different ways in which sex chromosome systems are organized at meiosis, how pairing is affected, and differences in unpaired DNA responses. We hypothesize that lineage specific differences in meiotic organization is not only a consequence of sex chromosome evolution, but that the establishment of epigenetic changes on sex chromosomes contributes toward their evolutionary conservation.

Evolutionary conservation of Y Chromosome ampliconic gene families despite extensive structural variation

Despite claims that the mammalian Y Chromosome is on a path to extinction, comparative sequence analysis of primate Y Chromosomes has shown the decay of the ancestral single-copy genes has all but ceased in this eutherian lineage. The suite of single-copy Y-linked genes is highly conserved among the majority of eutherian Y Chromosomes due to strong purifying selection to retain dosage-sensitive genes. In contrast, the ampliconic regions of the Y Chromosome, which contain testis-specific genes that encode the majority of the transcripts on eutherian Y Chromosomes, are rapidly evolving and are thought to undergo species-specific turnover. However, ampliconic genes are known from only a handful of species, limiting insights into their long-term evolutionary dynamics. We used a clone-based sequencing approach employing both long- and short-read sequencing technologies to assemble ∼2.4 Mb of representative ampliconic sequence dispersed across the domestic cat Y Chromosome, and identified the major ampliconic gene families and repeat units. We analyzed fluorescence in situ hybridization, qPCR, and whole-genome sequence data from 20 cat species and revealed that ampliconic gene families are conserved across the cat family Felidae but show high transcript diversity, copy number variation, and structural rearrangement. Our analysis of ampliconic gene evolution unveils a complex pattern of long-term gene content stability despite extensive structural variation on a nonrecombining background.

Genome-Wide Evolutionary Analysis of Natural History and Adaptation in the World's Tigers

No other species attracts more international resources, public attention, and protracted controversies over its intraspecific taxonomy than the tiger (Panthera tigris) [1, 2]. Today, fewer than 4,000 free-ranging tigers survive, covering only 7% of their historical range, and debates persist over whether they comprise six, five, or two subspecies [3-6]. The lack of consensus over the number of tiger subspecies has partially hindered the global effort to recover the species from the brink of extinction, as both captive breeding and landscape intervention of wild populations increasingly require an explicit delineation of the conservation management units [7]. The recent coalescence to a late Pleistocene bottleneck (circa 110 kya) [5, 8, 9] poses challenges for detecting tiger subspecific morphological traits, suggesting that elucidating intraspecific evolution in the tiger requires analyses at the genomic scale. Here, we present whole-genome sequencing analyses from 32 voucher specimens that resolve six statistically robust monophyletic clades corresponding to extant subspecies, including the recently recognized Malayan tiger (P. tigris jacksoni). The intersubspecies gene flow is very low, corroborating the recognized phylogeographic units. We identified multiple genomic regions that are candidates for identifying the adaptive divergence of subspecies. The body-size-related gene ADH7 appears to have been strongly selected in the Sumatran tiger, perhaps in association with adaptation to the tropical Sunda Islands. The identified genomic signatures provide a solid basis for recognizing appropriate conservation management units in the tiger and can benefit global conservation strategic planning for this charismatic megafauna icon.

Spermatogenesis and the Evolution of Mammalian Sex Chromosomes

Developmental constraint and sexual conflict shape the evolution of heteromorphic sex chromosomes. These contrasting forces are perhaps strongest during spermatogenesis in species with XY males. In this review, we consider how the unique regulatory environment and selective pressures of spermatogenesis interact to impact sex chromosome evolution in mammals. We explore how each developmental phase of spermatogenesis influences sex chromosome gene content, structure, and rate of molecular evolution, and how these attributes may contribute to speciation. We argue that a developmental context is fundamental to understanding sex chromosome evolution and that an evolutionary perspective can shed new light on our understanding of sperm development.

Deep Sequencing of Fosmid Clones Indicates Gene Conversion in the Male-Specific Region of the Giant Panda Y Chromosome

The giant panda (Ailuropoda melanoleuca) is popular around the world and is widely recognized as a symbol of nature conservation. A draft genome of the giant panda is now available, but its Y chromosome has not been sequenced. Y chromosome data are necessary for study of sex chromosome evolution, male development, and spermatogenesis. Thus, in the present study, we sequenced two parts of the giant panda Y chromosome utilizing a male giant panda fosmid library. The sequencing data were assembled into two contigs, each ∼100 kb in length with no gaps, providing high-quality resources for studying the giant panda Y chromosome. Annotation and transposable element comparison indicates varied evolutionary pressure in different regions of the Y chromosome. Two genes, zinc finger protein, Y-linked (ZFY) and lysine demethylase 5D (KDM5D), were annotated and gene conversion was observed for ZFY exon 7. Phylogenetic analysis also revealed that this gene conversion event happened independently in multiple mammalian lineages, indicating a putative mechanism to maintain the function of this particular gene on the Y chromosome. Furthermore, a transposition event, discovered through comparative alignment with the giant panda X chromosome sequence, may be involved in the process of gaining new genes on the Y chromosome. Thus, these newly obtained Y chromosome sequences provide valuable insights into the genomic patterns of the giant panda.

Red fox genome assembly identifies genomic regions associated with tame and aggressive behaviours

Strains of red fox (Vulpes vulpes) with markedly different behavioural phenotypes have been developed in the famous long-term selective breeding programme known as the Russian farm-fox experiment. Here we sequenced and assembled the red fox genome and re-sequenced a subset of foxes from the tame, aggressive and conventional farm-bred populations to identify genomic regions associated with the response to selection for behaviour. Analysis of the re-sequenced genomes identified 103 regions with either significantly decreased heterozygosity in one of the three populations or increased divergence between the populations. A strong positional candidate gene for tame behaviour was highlighted: SorCS1, which encodes the main trafficking protein for AMPA glutamate receptors and neurexins and suggests a role for synaptic plasticity in fox domestication. Other regions identified as likely to have been under selection in foxes include genes implicated in human neurological disorders, mouse behaviour and dog domestication. The fox represents a powerful model for the genetic analysis of affiliative and aggressive behaviours that can benefit genetic studies of behaviour in dogs and other mammals, including humans.

Horse Y chromosome assembly displays unique evolutionary features and putative stallion fertility genes

Dynamic evolutionary processes and complex structure make the Y chromosome among the most diverse and least understood regions in mammalian genomes. Here, we present an annotated assembly of the male specific region of the horse Y chromosome (eMSY), representing the first comprehensive Y assembly in odd-toed ungulates. The eMSY comprises single-copy, equine specific multi-copy, PAR transposed, and novel ampliconic sequence classes. The eMSY gene density approaches that of autosomes with the highest number of retained X-Y gametologs recorded in eutherians, in addition to novel Y-born and transposed genes. Horse, donkey and mule testis RNAseq reveals several candidate genes for stallion fertility. A novel testis-expressed XY ampliconic sequence class, ETSTY7, is shared with the parasite Parascaris genome, providing evidence for eukaryotic horizontal transfer and inter-chromosomal mobility. Our study highlights the dynamic nature of the Y and provides a reference sequence for improved understanding of equine male development and fertility.

Methodology for Y Chromosome Capture: A complete genome sequence of Y chromosome using flow cytometry, laser microdissection and magnetic streptavidin-beads

This study is a comparison of the efficiency of three technologies used for Y chromosome capture and the next-generation sequencing (NGS) technologies applied for determining its whole sequence. Our main findings disclose that streptavidin-biotin magnetic particle-based capture methodology offers better and a deeper sequence coverage for Y chromosome capture, compared to chromosome sorting and microdissection procedures. Moreover, this methodology is less time consuming and the most selective for capturing only Y chromosomal material, in contrast with other methodologies that result in considerable background material from other, non-targeted chromosomes. NGS results compared between two platforms, NextSeq 500 and SOLID 5500xl, produce the same coverage results. This is the first study to explore a methodological comparison of Y chromosome capture and genetic analysis. Our results indicate an improved strategy for Y chromosome research with applications in several scientific fields where this chromosome plays an important role, such as forensics, medical sciences, molecular anthropology and cancer sciences.

The Case of X and Y Localization of Nucleolus Organizer Regions (NORs) in Tragulus javanicus (Cetartiodactyla, Mammalia)

There are differences in number and localization of nucleolus organizer regions (NORs) in genomes. In mammalian genomes, NORs are located on autosomes, which are often situated on short arms of acrocentric chromosomes and more rarely in telomeric, pericentromeric, or interstitial regions. In this work, we report the unique case of active NORs located on gonоsomes of a eutherian mammal, the Javan mouse-deer (Tragulus javanicus). We have investigated the position of NORs by FISH experiments with ribosomal DNA (rDNA) sequences (18S, 5.8S, and 28S) and show the presence of a single NOR site on the X and Y chromosomes. The NOR is localized interstitially on the p-arm of the X chromosome in close proximity with prominent C-positive heterochromatin blocks and in the pericentromeric area of mostly heterochromatic Y. The NOR sites are active on both the X and Y chromosomes in the studied individual and surrounded by GC enriched heterochromatin. We hypothesize that the surrounding heterochromatin might have played a role in the transfer of NORs from autosomes to sex chromosomes during the karyotype evolution of the Javan mouse-deer.

Analysis of the canid Y-chromosome phylogeny using short-read sequencing data reveals the presence of distinct haplogroups among Neolithic European dogs

Background

Most genetic analyses of ancient and modern dogs have focused on variation in the autosomes or on the mitochondria. Mitochondrial DNA is more easily obtained from ancient samples than nuclear DNA and mitochondrial analyses have revealed important insights into the evolutionary history of canids. Utilizing a recently published dog Y-chromosome reference, we analyzed Y-chromosome sequence across a diverse collection of canids and determined the Y haplogroup of three ancient European dogs.

Results

We identified 1121 biallelic Y-chromosome SNVs using whole-genome sequences from 118 canids and defined variants diagnostic to distinct dog Y haplogroups. Similar to that of the mitochondria and previous more limited studies of Y diversity, we observe several deep splits in the Y-chromosome tree which may be the result of retained Y-chromosome diversity which predates dog domestication or post-domestication admixture with wolves. We find that Y-chromosomes from three ancient European dogs (4700–7000 years old) belong to distinct clades.

Conclusions

We estimate that the time to the most recent comment ancestor of dog Y haplogroups is 68–151 thousand years ago. Analysis of three Y-chromosomes from the Neolithic confirms long stranding population structure among European dogs.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4749-z) contains supplementary material, which is available to authorized users.

Cost-effective high-throughput single-haplotype iterative mapping and sequencing for complex genomic structures

The reference sequences of structurally complex regions can be obtained only through highly accurate clone-based approaches. We and others have successfully used single-haplotype iterative mapping and sequencing (SHIMS) 1.0 to assemble structurally complex regions across the sex chromosomes of several vertebrate species and to allow for targeted improvements to the reference sequences of human autosomes. However, SHIMS 1.0 is expensive and time consuming, requiring resources that only a genome center can provide. Here we introduce SHIMS 2.0, an improved SHIMS protocol that allows even a small laboratory to generate high-quality reference sequence from complex genomic regions. Using a streamlined and parallelized library-preparation protocol, and taking advantage of inexpensive high-throughput short-read-sequencing technologies, a small laboratory with both molecular biology and bioinformatics experience can sequence and assemble 192 large-insert bacterial artificial chromosome (BAC) or fosmid clones in 1 week. In SHIMS 2.0, in contrast to other pooling strategies, each clone is sequenced with a unique barcode, thus enabling clones containing nearly identical sequences to be multiplexed in a single sequencing run and assembled separately. Relative to SHIMS 1.0, SHIMS 2.0 decreases the required cost and time by two orders of magnitude while preserving high sequencing accuracy.

Genetic sex assignment in wild populations using genotyping-by-sequencing data: A statistical threshold approach

Establishing the sex of individuals in wild systems can be challenging and often requires genetic testing. Genotyping-by-sequencing (GBS) and other reduced-representation DNA sequencing (RRS) protocols (e.g., RADseq, ddRAD) have enabled the analysis of genetic data on an unprecedented scale. Here, we present a novel approach for the discovery and statistical validation of sex-specific loci in GBS data sets. We used GBS to genotype 166 New Zealand fur seals (NZFS, Arctocephalus forsteri) of known sex. We retained monomorphic loci as potential sex-specific markers in the locus discovery phase. We then used (i) a sex-specific locus threshold (SSLT) to identify significantly male-specific loci within our data set; and (ii) a significant sex-assignment threshold (SSAT) to confidently assign sex in silico the presence or absence of significantly male-specific loci to individuals in our data set treated as unknowns (98.9% accuracy for females; 95.8% for males, estimated via cross-validation). Furthermore, we assigned sex to 86 individuals of true unknown sex using our SSAT and assessed the effect of SSLT adjustments on these assignments. From 90 verified sex-specific loci, we developed a panel of three sex-specific PCR primers that we used to ascertain sex independently of our GBS data, which we show amplify reliably in at least two other pinniped species. Using monomorphic loci normally discarded from large SNP data sets is an effective way to identify robust sex-linked markers for nonmodel species. Our novel pipeline can be used to identify and statistically validate monomorphic and polymorphic sex-specific markers across a range of species and RRS data sets.

Y-Chromosome Markers for the Red Fox

The de novo assembly of the red fox (Vulpes vulpes) genome has facilitated the development of genomic tools for the species. Efforts to identify the population history of red foxes in North America have previously been limited by a lack of information about the red fox Y-chromosome sequence. However, a megabase of red fox Y-chromosome sequence was recently identified over 2 scaffolds in the reference genome. Here, these scaffolds were scanned for repeated motifs, revealing 194 likely microsatellites. Twenty-three of these loci were selected for primer development and, after testing, produced a panel of 11 novel markers that were analyzed alongside 2 markers previously developed for the red fox from dog Y-chromosome sequence. The markers were genotyped in 76 male red foxes from 4 populations: 7 foxes from Newfoundland (eastern Canada), 12 from Maryland (eastern United States), and 9 from the island of Great Britain, as well as 48 foxes of known North American origin maintained on an experimental farm in Novosibirsk, Russia. The full marker panel revealed 22 haplotypes among these red foxes, whereas the 2 previously known markers alone would have identified only 10 haplotypes. The haplotypes from the 4 populations clustered primarily by continent, but unidirectional gene flow from Great Britain and farm populations may influence haplotype diversity in the Maryland population. The development of new markers has increased the resolution at which red fox Y-chromosome diversity can be analyzed and provides insight into the contribution of males to red fox population diversity and patterns of phylogeography.

Abundance and Characterization of Perfect Microsatellites on the Cattle Y Chromosome

Microsatellites or simple sequence repeats (SSRs) are found in most organisms and play an important role in genomic organization and function. To characterize the abundance of SSRs (1-6 base-pairs [bp]) on the cattle Y chromsome, the relative frequency and density of perfect or uninterrupted SSRs based on the published Y chromosome sequence were examined. A total of 17,273 perfect SSRs were found, with total length of 324.78 kb, indicating that approximately 0.75% of the cattle Y chromosome sequence (43.30 Mb) comprises perfect SSRs, with an average length of 18.80 bp. The relative frequency and density were 398.92 loci/Mb and 7500.62 bp/Mb, respectively. The proportions of the six classes of perfect SSRs were highly variable on the cattle Y chromosome. Mononucleotide repeats had a total number of 8073 (46.74%) and an average length of 15.45 bp, and were the most abundant SSRs class, while the percentages of di-, tetra-, tri-, penta-, and hexa-nucleotide repeats were 22.86%, 11.98%, 11.58%, 6.65%, and 0.19%, respectively. Different classes of SSRs varied in their repeat number, with the highest being 42 for dinucleotides. Results reveal that repeat categories A, AC, AT, AAC, AGC, GTTT, CTTT, ATTT, and AACTG predominate on the Y chromosome. This study provides insight into the organization of cattle Y chromosome repetitive DNA, as well as information useful for developing more polymorphic cattle Y-chromosome-specific SSRs.

Highly dynamic temporal changes of TSPY gene copy number in aging bulls

The Y-chromosomal TSPY gene is one of the highest copy number mammalian protein coding gene and represents a unique biological model to study various aspects of genomic copy number variations. This study investigated the age-related copy number variability of the bovine TSPY gene, a new and unstudied aspect of the biology of TSPY that has been shown to vary among cattle breeds, individual bulls and somatic tissues. The subjects of this prospective 30-month long study were 25 Holstein bulls, sampled every six months. Real-time quantitative PCR was used to determine the relative TSPY copy number (rTSPY CN) and telomere length in the DNA samples extracted from blood. Twenty bulls showed an altered rTSPY CN after 30 months, although only 9 bulls showed a significant change (4 significant increase while 5 significant decrease, P<0.01). The sequential sampling provided the flow of rTSPY CN over six observations in 30 months and wide-spread variation of rTSPY CN was detected. Although a clear trend of the direction of change was not identifiable, the highly dynamic changes of individual rTSPY CN in aging bulls were observed here for the first time. In summary we have observed a highly variable rTSPY CN in bulls over a short period of time. Our results suggest the importance of further long term studies of the dynamics of rTSPY CN variablility.

Patterns of Inter-Chromosomal Gene Conversion on the Male-Specific Region of the Human Y Chromosome

The male-specific region of the human Y chromosome (MSY) is characterized by the lack of meiotic recombination and it has long been considered an evolutionary independent region of the human genome. In recent years, however, the idea that human MSY did not have an independent evolutionary history begun to emerge with the discovery that inter-chromosomal gene conversion (ICGC) can modulate the genetic diversity of some portions of this genomic region. Despite the study of the dynamics of this molecular mechanism in humans is still in its infancy, some peculiar features and consequences of it can be summarized. The main effect of ICGC is to increase the allelic diversity of MSY by generating a significant excess of clustered single nucleotide polymorphisms (SNPs) (defined as groups of two or more SNPs occurring in close proximity and on the same branch of the Y phylogeny). On the human MSY, 13 inter-chromosomal gene conversion hotspots (GCHs) have been identified so far, involving donor sequences mainly from the X-chromosome and, to a lesser extent, from autosomes. Most of the GCHs are evolutionary conserved and overlap with regions involved in aberrant X–Y crossing-over. This review mainly focuses on the dynamics and the current knowledge concerning the recombinational landscape of the human MSY in the form of ICGC, on how this molecular mechanism may influence the evolution of the MSY, and on how it could affect the information enclosed within a genomic region which, until recently, appeared to be an evolutionary independent unit.

Y chromosome palindromes and gene conversion

The presence of large and near-identical inverted repeat sequences (called palindromes) is a common feature of the constitutively haploid sex chromosomes of different species. Despite the fact palindromes originated in a non-recombining context, they have evolved a strong recombinational activity in the form of abundant arm-to-arm gene conversion. Their independent appearance in different species suggests they can have a profound biological significance that has yet to be fully clarified. It has been theorized that natural selection may have favored palindromic organization of male-specific genes and that the establishment of intra-palindrome gene conversion has strong adaptive significance. Arm-to-arm gene conversion allows the efficient removal of deleterious mutations, increases the fixation rate of beneficial mutations and has played an important role in modulating the equilibrium between gene loss and acquisition during Y chromosome evolution. Additionally, a palindromic organization of duplicates could favor the formation of unusual chromatin structures and could optimize the use of gene conversion as a mechanism to maintain the structural integrity of male-specific genes. In this review, we describe the structural features of palindromes on mammalian sex chromosomes and summarize different hypotheses regarding palindrome evolution and the functional benefits of arm-to-arm gene conversion on the unique haploid portion of the nuclear genome.

The Y chromosomes of the great apes

The great apes (orangutans, gorillas, chimpanzees, bonobos and humans) descended from a common ancestor around 13 million years ago, and since then their sex chromosomes have followed very different evolutionary paths. While great-ape X chromosomes are highly conserved, their Y chromosomes, reflecting the general lability and degeneration of this male-specific part of the genome since its early mammalian origin, have evolved rapidly both between and within species. Understanding great-ape Y chromosome structure, gene content and diversity would provide a valuable evolutionary context for the human Y, and would also illuminate sex-biased behaviours, and the effects of the evolutionary pressures exerted by different mating strategies on this male-specific part of the genome. High-quality Y-chromosome sequences are available for human and chimpanzee (and low-quality for gorilla). The chromosomes differ in size, sequence organisation and content, and while retaining a relatively stable set of ancestral single-copy genes, show considerable variation in content and copy number of ampliconic multi-copy genes. Studies of Y-chromosome diversity in other great apes are relatively undeveloped compared to those in humans, but have nevertheless provided insights into speciation, dispersal, and mating patterns. Future studies, including data from larger sample sizes of wild-born and geographically well-defined individuals, and full Y-chromosome sequences from bonobos, gorillas and orangutans, promise to further our understanding of population histories, male-biased behaviours, mutation processes, and the functions of Y-chromosomal genes.

Avian W and mammalian Y chromosomes convergently retained dosage-sensitive regulators

After birds diverged from mammals, different ancestral autosomes evolved into sex chromosomes in each lineage. In birds, females are ZW and males ZZ, but in mammals females are XX and males XY. We sequenced the chicken W chromosome, compared its gene content with our reconstruction of the ancestral autosomes, and followed the evolutionary trajectory of ancestral W-linked genes across birds. Avian W chromosomes evolved in parallel with mammalian Y chromosomes, preserving ancestral genes through selection to maintain the dosage of broadly-expressed regulators of key cellular processes. We propose that, like the human Y chromosome, the chicken W chromosome is essential for embryonic viability of the heterogametic sex. Unlike other sequenced sex chromosomes, the chicken W did not acquire and amplify genes specifically expressed in reproductive tissues. We speculate that the pressures that drive the acquisition of reproduction related genes on sex chromosomes may be specific to the male germ line.

Chromosome Aberrations and Fertility Disorders in Domestic Animals

The association between chromosomal abnormalities and reduced fertility in domestic animals is well recorded and has been studied for decades. Chromosome aberrations directly affect meiosis, gametogenesis, and the viability of zygotes and embryos. In some instances, balanced structural rearrangements can be transmitted, causing fertility problems in subsequent generations. Here, we aim to give a comprehensive overview of the current status and future prospects of clinical cytogenetics of animal reproduction by focusing on the advances in molecular cytogenetics during the genomics era. We describe how advancing knowledge about animal genomes has improved our understanding of connections between gross structural or molecular chromosome variations and reproductive disorders. Further, we expand on a key area of reproduction genetics: cytogenetics of animal gametes and embryos. Finally, we describe how traditional cytogenetics is interfacing with advanced genomics approaches, such as array technologies and next-generation sequencing, and speculate about the future prospects.

Genome sequence, population history, and pelage genetics of the endangered African wild dog (Lycaon pictus)

BACKGROUND

The African wild dog (Lycaon pictus) is an endangered African canid threatened by severe habitat fragmentation, human-wildlife conflict, and infectious disease. A highly specialized carnivore, it is distinguished by its social structure, dental morphology, absence of dewclaws, and colorful pelage.

RESULTS

We sequenced the genomes of two individuals from populations representing two distinct ecological histories (Laikipia County, Kenya and KwaZulu-Natal Province, South Africa). We reconstructed population demographic histories for the two individuals and scanned the genomes for evidence of selection.

CONCLUSIONS

We show that the African wild dog has undergone at least two effective population size reductions in the last 1,000,000 years. We found evidence of Lycaon individual-specific regions of low diversity, suggestive of inbreeding or population-specific selection. Further research is needed to clarify whether these population reductions and low diversity regions are characteristic of the species as a whole. We documented positive selection on the Lycaon mitochondrial genome. Finally, we identified several candidate genes (ASIP, MITF, MLPH, PMEL) that may play a role in the characteristic Lycaon pelage.

Expression and epigenomic landscape of the sex chromosomes in mouse post-meiotic male germ cells

BACKGROUND

During meiosis, the X and Y chromosomes are transcriptionally silenced. The persistence of repressive chromatin marks on the sex chromatin after meiosis initially led to the assumption that XY gene silencing persists to some extent in spermatids. Considering the many reports of XY-linked genes expressed and needed in the post-meiotic phase of mouse spermatogenesis, it is still unclear whether or not the mouse sex chromatin is a repressive or permissive environment, after meiosis.

RESULTS

To determine the transcriptional and chromatin state of the sex chromosomes after meiosis, we re-analyzed ten ChIP-Seq datasets performed on mouse round spermatids and four RNA-seq datasets from male germ cells purified at different stages of spermatogenesis. For this, we used the last version of the genome (mm10/GRCm38) and included reads that map to several genomic locations in order to properly interpret the high proportion of sex chromosome-encoded multicopy genes. Our study shows that coverage of active epigenetic marks H3K4me3 and Kcr is similar on the sex chromosomes and on autosomes. The post-meiotic sex chromatin nevertheless differs from autosomal chromatin in its enrichment in H3K9me3 and its depletion in H3K27me3 and H4 acetylation. We also identified a posttranslational modification, H3K27ac, which specifically accumulates on the Y chromosome. In parallel, we found that the X and Y chromosomes are enriched in genes expressed post-meiotically and display a higher proportion of spermatid-specific genes compared to autosomes. Finally, we observed that portions of chromosome 14 and of the sex chromosomes share specific features, such as enrichment in H3K9me3 and the presence of multicopy genes that are specifically expressed in round spermatids, suggesting that parts of chromosome 14 are under the same evolutionary constraints than the sex chromosomes.

CONCLUSIONS

Based on our expression and epigenomic studies, we conclude that, after meiosis, the mouse sex chromosomes are no longer silenced but are nevertheless regulated differently than autosomes and accumulate different chromatin marks. We propose that post-meiotic selective constraints are at the basis of the enrichment of spermatid-specific genes and of the peculiar chromatin composition of the sex chromosomes and of parts of chromosome 14.

A High-Resolution SNP Array-Based Linkage Map Anchors a New Domestic Cat Draft Genome Assembly and Provides Detailed Patterns of Recombination

High-resolution genetic and physical maps are invaluable tools for building accurate genome assemblies, and interpreting results of genome-wide association studies (GWAS). Previous genetic and physical maps anchored good quality draft assemblies of the domestic cat genome, enabling the discovery of numerous genes underlying hereditary disease and phenotypes of interest to the biomedical science and breeding communities. However, these maps lacked sufficient marker density to order thousands of shorter scaffolds in earlier assemblies, which instead relied heavily on comparative mapping with related species. A high-resolution map would aid in validating and ordering chromosome scaffolds from existing and new genome assemblies. Here, we describe a high-resolution genetic linkage map of the domestic cat genome based on genotyping 453 domestic cats from several multi-generational pedigrees on the Illumina 63K SNP array. The final maps include 58,055 SNP markers placed relative to 6637 markers with unique positions, distributed across all autosomes and the X chromosome. Our final sex-averaged maps span a total autosomal length of 4464 cM, the longest described linkage map for any mammal, confirming length estimates from a previous microsatellite-based map. The linkage map was used to order and orient the scaffolds from a substantially more contiguous domestic cat genome assembly (Felis catus v8.0), which incorporated ∼20 × coverage of Illumina fragment reads. The new genome assembly shows substantial improvements in contiguity, with a nearly fourfold increase in N50 scaffold size to 18 Mb. We use this map to report probable structural errors in previous maps and assemblies, and to describe features of the recombination landscape, including a massive (∼50 Mb) recombination desert (of virtually zero recombination) on the X chromosome that parallels a similar desert on the porcine X chromosome in both size and physical location.

A time- and cost-effective strategy to sequence mammalian Y Chromosomes: an application to the de novo assembly of gorilla Y

The mammalian Y Chromosome sequence, critical for studying male fertility and dispersal, is enriched in repeats and palindromes, and thus, is the most difficult component of the genome to assemble. Previously, expensive and labor-intensive BAC-based techniques were used to sequence the Y for a handful of mammalian species. Here, we present a much faster and more affordable strategy for sequencing and assembling mammalian Y Chromosomes of sufficient quality for most comparative genomics analyses and for conservation genetics applications. The strategy combines flow sorting, short- and long-read genome and transcriptome sequencing, and droplet digital PCR with novel and existing computational methods. It can be used to reconstruct sex chromosomes in a heterogametic sex of any species. We applied our strategy to produce a draft of the gorilla Y sequence. The resulting assembly allowed us to refine gene content, evaluate copy number of ampliconic gene families, locate species-specific palindromes, examine the repetitive element content, and produce sequence alignments with human and chimpanzee Y Chromosomes. Our results inform the evolution of the hominine (human, chimpanzee, and gorilla) Y Chromosomes. Surprisingly, we found the gorilla Y Chromosome to be similar to the human Y Chromosome, but not to the chimpanzee Y Chromosome. Moreover, we have utilized the assembled gorilla Y Chromosome sequence to design genetic markers for studying the male-specific dispersal of this endangered species.

The cryptic Y-autosome translocation in the small Indian mongoose, Herpestes auropunctatus, revealed by molecular cytogenetic approaches

In initial studies of the eutherian small Indian mongoose (Herpestes auropunctatus), the Y chromosome could not be identified in somatic cells. The male chromosome number is uniquely odd, 2n = 35, whereas that of females is 2n = 36. Previous reports indicated that this unique karyotype resulted from a translocation of the ancestral Y chromosome to an autosome. However, it has been difficult to identify the chromosomes that harbor the translocated Y chromosomal segment because it is an extremely small euchromatic region. Using a Southern blot analysis, we detected four conserved Y-linked genes, SRY, EIF2S3Y, KDM5D, and ZFY, in the male genome. We cloned homologues of these genes and determined their sequences, which showed high homology to genes in two carnivore species, cat and dog. To unambiguously identify the Y-bearing autosome, we performed immunostaining of pachytene spermatocytes using antibodies against SYCP3, γH2AX, and the centromere. We observed trivalent chromosomes, and the associations between the distal ends of the chromosomes were consistent with those of Y and X1 chromosomes. The centromere of the Y chromosome was located on the ancestral Y chromosomal segment. We mapped the complementary DNA (cDNA) clones of these genes to the male chromosomes using fluorescence in situ hybridization (FISH), and the linear localization of all genes was confirmed by two-colored FISH. These Y-linked genes were localized to the proximal region of the long arm of a single telomeric chromosome, and we successfully identified the chromosome harboring the ancestral Y chromosomal segment.

The Eutherian Pseudoautosomal Region

The pseudoautosomal region (PAR) is a unique segment of sequence homology between differentiated sex chromosomes where recombination occurs during meiosis. Molecular and functional properties of the PAR are distinctive from the autosomes and the remaining regions of the sex chromosomes. These include a higher rate of recombination than genome average, bias towards GC-substitutions and increased interindividual nucleotide divergence and mutations. As yet, the PAR has been physically demarcated in only 28 eutherian species representing 6 mammalian orders. Murid rodents have the smallest, gene-poorest and most diverged PARs. Other eutherian PARs are largely homologous but differ in size and gene content, being the smallest in equids and human/simian primates and much larger in other eutherians. Because pseudoautosomal genes escape X inactivation, their dosage changes with sex chromosome aneuploidies, whereas phenotypic effects of the latter depend on the size and gene content of the PAR. Thus, X monosomy is more viable in mice, humans and horses than in species with larger PARs. Presently, little is known about the functions of PAR genes in individual species, though human studies suggest their involvement in early embryonic development. The PAR is, thus, of evolutionary, genetic and biomedical significance and a 'research hotspot' in eutherian genomes.

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.

Mechanisms Underlying Mammalian Hybrid Sterility in Two Feline Interspecies Models

The phenomenon of male sterility in interspecies hybrids has been observed for over a century, however, few genes influencing this recurrent phenotype have been identified. Genetic investigations have been primarily limited to a small number of model organisms, thus limiting our understanding of the underlying molecular basis of this well-documented "rule of speciation." We utilized two interspecies hybrid cat breeds in a genome-wide association study employing the Illumina 63 K single-nucleotide polymorphism array. Collectively, we identified eight autosomal genes/gene regions underlying associations with hybrid male sterility (HMS) involved in the function of the blood-testis barrier, gamete structural development, and transcriptional regulation. We also identified several candidate hybrid sterility regions on the X chromosome, with most residing in close proximity to complex duplicated regions. Differential gene expression analyses revealed significant chromosome-wide upregulation of X chromosome transcripts in testes of sterile hybrids, which were enriched for genes involved in chromatin regulation of gene expression. Our expression results parallel those reported in Mus hybrids, supporting the "Large X-Effect" in mammalian HMS and the potential epigenetic basis for this phenomenon. These results support the value of the interspecies feline model as a powerful tool for comparison to rodent models of HMS, demonstrating unique aspects and potential commonalities that underpin mammalian reproductive isolation.

Estimating tempo and mode of Y chromosome turnover: explaining Y chromosome loss with the fragile Y hypothesis

Chromosomal sex determination is phylogenetically widespread, having arisen independently in many lineages. Decades of theoretical work provide predictions about sex chromosome differentiation that are well supported by observations in both XY and ZW systems. However, the phylogenetic scope of previous work gives us a limited understanding of the pace of sex chromosome gain and loss and why Y or W chromosomes are more often lost in some lineages than others, creating XO or ZO systems. To gain phylogenetic breadth we therefore assembled a database of 4724 beetle species' karyotypes and found substantial variation in sex chromosome systems. We used the data to estimate rates of Y chromosome gain and loss across a phylogeny of 1126 taxa estimated from seven genes. Contrary to our initial expectations, we find that highly degenerated Y chromosomes of many members of the suborder Polyphaga are rarely lost, and that cases of Y chromosome loss are strongly associated with chiasmatic segregation during male meiosis. We propose the "fragile Y" hypothesis, that recurrent selection to reduce recombination between the X and Y chromosome leads to the evolution of a small pseudoautosomal region (PAR), which, in taxa that require XY chiasmata for proper segregation during meiosis, increases the probability of aneuploid gamete production, with Y chromosome loss. This hypothesis predicts that taxa that evolve achiasmatic segregation during male meiosis will rarely lose the Y chromosome. We discuss data from mammals, which are consistent with our prediction.

An improved model for whole genome phylogenetic analysis by Fourier transform

DNA sequence similarity comparison is one of the major steps in computational phylogenetic studies. The sequence comparison of closely related DNA sequences and genomes is usually performed by multiple sequence alignments (MSA). While the MSA method is accurate for some types of sequences, it may produce incorrect results when DNA sequences undergone rearrangements as in many bacterial and viral genomes. It is also limited by its computational complexity for comparing large volumes of data. Previously, we proposed an alignment-free method that exploits the full information contents of DNA sequences by Discrete Fourier Transform (DFT), but still with some limitations. Here, we present a significantly improved method for the similarity comparison of DNA sequences by DFT. In this method, we map DNA sequences into 2-dimensional (2D) numerical sequences and then apply DFT to transform the 2D numerical sequences into frequency domain. In the 2D mapping, the nucleotide composition of a DNA sequence is a determinant factor and the 2D mapping reduces the nucleotide composition bias in distance measure, and thus improving the similarity measure of DNA sequences. To compare the DFT power spectra of DNA sequences with different lengths, we propose an improved even scaling algorithm to extend shorter DFT power spectra to the longest length of the underlying sequences. After the DFT power spectra are evenly scaled, the spectra are in the same dimensionality of the Fourier frequency space, then the Euclidean distances of full Fourier power spectra of the DNA sequences are used as the dissimilarity metrics. The improved DFT method, with increased computational performance by 2D numerical representation, can be applicable to any DNA sequences of different length ranges. We assess the accuracy of the improved DFT similarity measure in hierarchical clustering of different DNA sequences including simulated and real datasets. The method yields accurate and reliable phylogenetic trees and demonstrates that the improved DFT dissimilarity measure is an efficient and effective similarity measure of DNA sequences. Due to its high efficiency and accuracy, the proposed DFT similarity measure is successfully applied on phylogenetic analysis for individual genes and large whole bacterial genomes.

The Porcine TSPY Gene Is Tricopy but Not a Copy Number Variant

The testis-specific protein Y-encoded (TSPY) gene is situated on the mammalian Y-chromosome and exhibits some remarkable biological characteristics. It has the highest known copy number (CN) of all protein coding genes in the human and bovine genomes (up to 74 and 200, respectively) and also shows high individual variability. Although the biological function of TSPY has not yet been elucidated, its specific expression in the testis and several identified binding domains within the protein suggests roles in male reproduction. Here we describe the porcine TSPY, as a multicopy gene with three copies located on the short arm of the Y-chromosome with no variation at three exon loci among 20 animals of normal reproductive health from four breeds of domestic pigs (Piétrain, Landrace, Duroc and Yorkshire). To further investigate the speculation that porcine TSPY is not a copy number variant, we have included five Low-fertility boars and five boars with exceptional High-fertility records. Interestingly, there was no difference between the High- and Low-fertile groups, but we detected slightly lower TSPY CN at all three exons (2.56-2.85) in both groups, as compared to normal animals, which could be attributed to technical variability or somatic mosaicism. The results are based on both relative quantitative real-time PCR (qPCR) and droplet digital PCR (ddPCR). Chromosomal localization of the porcine TSPY was done using fluorescence in situ hybridization (FISH) with gene specific PCR probes.

Genome-Wide Search Identifies 1.9 Mb from the Polar Bear Y Chromosome for Evolutionary Analyses

The male-inherited Y chromosome is the major haploid fraction of the mammalian genome, rendering Y-linked sequences an indispensable resource for evolutionary research. However, despite recent large-scale genome sequencing approaches, only a handful of Y chromosome sequences have been characterized to date, mainly in model organisms. Using polar bear (Ursus maritimus) genomes, we compare two different in silico approaches to identify Y-linked sequences: 1) Similarity to known Y-linked genes and 2) difference in the average read depth of autosomal versus sex chromosomal scaffolds. Specifically, we mapped available genomic sequencing short reads from a male and a female polar bear against the reference genome and identify 112 Y-chromosomal scaffolds with a combined length of 1.9 Mb. We verified the in silico findings for the longer polar bear scaffolds by male-specific in vitro amplification, demonstrating the reliability of the average read depth approach. The obtained Y chromosome sequences contain protein-coding sequences, single nucleotide polymorphisms, microsatellites, and transposable elements that are useful for evolutionary studies. A high-resolution phylogeny of the polar bear patriline shows two highly divergent Y chromosome lineages, obtained from analysis of the identified Y scaffolds in 12 previously published male polar bear genomes. Moreover, we find evidence of gene conversion among ZFX and ZFY sequences in the giant panda lineage and in the ancestor of ursine and tremarctine bears. Thus, the identification of Y-linked scaffold sequences from unordered genome sequences yields valuable data to infer phylogenomic and population-genomic patterns in bears.

Incomplete meiotic sex chromosome inactivation in the domestic dog

BACKGROUND

In mammalian meiotic prophase, homologous chromosome recognition is aided by formation and repair of programmed DNA double-strand breaks (DSBs). Subsequently, stable associations form through homologous chromosome synapsis. In male mouse meiosis, the largely heterologous X and Y chromosomes synapse only in their short pseudoautosomal regions (PARs), and DSBs persist along the unsynapsed non-homologous arms of these sex chromosomes. Asynapsis of these arms and the persistent DSBs then trigger transcriptional silencing through meiotic sex chromosome inactivation (MSCI), resulting in formation of the XY body. This inactive state is partially maintained in post-meiotic haploid spermatids (postmeiotic sex chromatin repression, PSCR). For the human, establishment of MSCI and PSCR have also been reported, but X-linked gene silencing appears to be more variable compared to mouse. To gain more insight into the regulation and significance of MSCI and PSCR among different eutherian species, we have performed a global analysis of XY pairing dynamics, DSB repair, MSCI and PSCR in the domestic dog (Canis lupus familiaris), for which the complete genome sequence has recently become available, allowing a thorough comparative analyses.

RESULTS

In addition to PAR synapsis between X and Y, we observed extensive self-synapsis of part of the dog X chromosome, and rapid loss of known markers of DSB repair from that part of the X. Sequencing of RNA from purified spermatocytes and spermatids revealed establishment of MSCI. However, the self-synapsing region of the X displayed higher X-linked gene expression compared to the unsynapsed area in spermatocytes, and was post-meiotically reactivated in spermatids. In contrast, genes in the PAR, which are expected to escape MSCI, were expressed at very low levels in both spermatocytes and spermatids. Our comparative analysis was then used to identify two X-linked genes that may escape MSCI in spermatocytes, and 21 that are specifically re-activated in spermatids of human, mouse and dog.

CONCLUSIONS

Our data indicate that MSCI is incomplete in the dog. This may be partially explained by extensive, but transient, self-synapsis of the X chromosome, in association with rapid completion of meiotic DSB repair. In addition, our comparative analysis identifies novel candidate male fertility genes.

Recent Y chromosome divergence despite ancient origin of dioecy in poplars (Populus)

All species of the genus Populus (poplar, aspen) are dioecious, suggesting an ancient origin of this trait. Despite some empirical counter examples, theory suggests that nonrecombining sex-linked regions should quickly spread, eventually becoming heteromorphic chromosomes. In contrast, we show using whole-genome scans that the sex-associated region in Populus trichocarpa is small and much younger than the age of the genus. This indicates that sex determination is highly labile in poplar, consistent with recent evidence of 'turnover' of sex-determination regions in animals. We performed whole-genome resequencing of 52 P. trichocarpa (black cottonwood) and 34 Populus balsamifera (balsam poplar) individuals of known sex. Genomewide association studies in these unstructured populations identified 650 SNPs significantly associated with sex. We estimate the size of the sex-linked region to be ~100 kbp. All SNPs significantly associated with sex were in strong linkage disequilibrium despite the fact that they were mapped to six different chromosomes (plus 3 unmapped scaffolds) in version 2.2 of the reference genome. We show that this is likely due to genome misassembly. The segregation pattern of sex-associated SNPs revealed this to be an XY sex-determining system. Estimated divergence times of X and Y haplotype sequences (6-7 Ma) are much more recent than the divergence of P. trichocarpa (poplar) and Populus tremuloides (aspen). Consistent with this, in P. tremuloides, we found no XY haplotype divergence within the P. trichocarpa sex-determining region. These two species therefore have a different genomic architecture of sex, suggestive of at least one turnover event in the recent past.

The first case of 38,XX (SRY-positive) disorder of sex development in a cat

Background

SRY-positive XX testicular disorder of sex development (DSD) caused by X;Y translocations was not yet reported in domestic animals. In humans it is rarely diagnosed and a majority of clinical features resemble those which are typical for Klinefelter syndrome (KS). Here we describe the first case of SRY-positive XX DSD in a tortoiseshell cat with a rudimentary penis and a lack of scrotum.

Results

Molecular analysis showed the presence of two Y-linked genes (SRY and ZFY) and a normal sequence of the SRY gene. Application of classical cytogenetic techniques revealed two X chromosomes (38,XX), but further FISH studies with the use of the whole X chromosome painting probe and BAC probes specific to the Yp chromosome facilitated identification of Xp;Yp translocation. The SRY gene was localised at a distal position of Xp. The karyotype of the studied case was described as: 38,XX.ish der(X)t(X;Y)(p22;p12)(SRY+). Moreover, the X inactivation status assessed by a sequential R-banding and FISH with the SRY-specific probe showed a random inactivation of the derivative X^SRY chromosome.

Conclusions

Our study showed that among DSD tortoiseshell cats, apart from XXY trisomy and XX/XY chimerism, also SRY-positive XX cases may occur. It is hypothesized that the extremely rare occurrence of this abnormality in domestic animals, when compared with humans, may be associated with a different organisation of the Yp arm in these species.

DNA mutations of the cat: the good, the bad and the ugly

PRACTICAL RELEVANCE

The health of the cat is a complex interaction between its environment (nurture) and its genetics (nature). Over 70 genetic mutations (variants) have been defined in the cat, many involving diseases, structural abnormalities and clinically relevant health concerns. As more of the cat's genome is deciphered, less commonly will the term 'idiopathic' be used regarding the diagnosis of diseases and unique health conditions. State-of-the-art health care will include DNA profiling of the individual cat, and perhaps its tumor, to establish the best treatment approaches. Genetic testing and eventually whole genome sequencing should become routine diagnostics for feline health care.

GLOBAL IMPORTANCE

Cat breeds have disseminated around the world. Thus, practitioners should be aware of the breeds common to their region and the mutations found in those regional populations. Specific random-bred populations can also have defined genetic characteristics and mutations.

AUDIENCE

This review of 'the good, the bad and the ugly' DNA variants provides the current state of knowledge for genetic testing and genetic health management for cats. It is aimed at feline and general practitioners wanting to update and review the basics of genetics, what tests are available for cats and sources for genetic testing. The tables are intended to be used as references in the clinic. Practitioners with a high proportion of cat breeder clientele will especially benefit from the review.

EVIDENCE BASE

The data presented is extracted from peer-reviewed publications pertaining to mutation identification, and relevant articles concerning the heritable trait and/or disease. The author also draws upon personal experience and expertise in feline genetics.

Complex evolutionary trajectories of sex chromosomes across bird taxa

Sex-specific chromosomes, like the W of most female birds and the Y of male mammals, usually have lost most genes owing to a lack of recombination. We analyze newly available genomes of 17 bird species representing the avian phylogenetic range, and find that more than half of them do not have as fully degenerated W chromosomes as that of chicken. We show that avian sex chromosomes harbor tremendous diversity among species in their composition of pseudoautosomal regions and degree of Z/W differentiation. Punctuated events of shared or lineage-specific recombination suppression have produced a gradient of "evolutionary strata" along the Z chromosome, which initiates from the putative avian sex-determining gene DMRT1 and ends at the pseudoautosomal region. W-linked genes are subject to ongoing functional decay after recombination was suppressed, and the tempo of degeneration slows down in older strata. Overall, we unveil a complex history of avian sex chromosome evolution.

Supergenes and complex phenotypes

Understanding the molecular underpinnings of evolutionary adaptations is a central focus of modern evolutionary biology. Recent studies have uncovered a panoply of complex phenotypes, including locally adapted ecotypes and cryptic morphs, divergent social behaviours in birds and insects, as well as alternative metabolic pathways in plants and fungi, that are regulated by clusters of tightly linked loci. These 'supergenes' segregate as stable polymorphisms within or between natural populations and influence ecologically relevant traits. Some supergenes may span entire chromosomes, because selection for reduced recombination between a supergene and a nearby locus providing additional benefits can lead to locus expansions with dynamics similar to those known for sex chromosomes. In addition to allowing for the co-segregation of adaptive variation within species, supergenes may facilitate the spread of complex phenotypes across species boundaries. Application of new genomic methods is likely to lead to the discovery of many additional supergenes in a broad range of organisms and reveal similar genetic architectures for convergently evolved phenotypes.

Sequencing the mouse Y chromosome reveals convergent gene acquisition and amplification on both sex chromosomes

We sequenced the MSY (male-specific region of the Y chromosome) of the C57BL/6J strain of the laboratory mouse Mus musculus. In contrast to theories that Y chromosomes are heterochromatic and gene poor, the mouse MSY is 99.9% euchromatic and contains about 700 protein-coding genes. Only 2% of the MSY derives from the ancestral autosomes that gave rise to the mammalian sex chromosomes. Instead, all but 45 of the MSY's genes belong to three acquired, massively amplified gene families that have no homologs on primate MSYs but do have acquired, amplified homologs on the mouse X chromosome. The complete mouse MSY sequence brings to light dramatic forces in sex chromosome evolution: lineage-specific convergent acquisition and amplification of X-Y gene families, possibly fueled by antagonism between acquired X-Y homologs. The mouse MSY sequence presents opportunities for experimental studies of a sex-specific chromosome in its entirety, in a genetically tractable model organism.

Copy number variation in the horse genome

We constructed a 400K WG tiling oligoarray for the horse and applied it for the discovery of copy number variations (CNVs) in 38 normal horses of 16 diverse breeds, and the Przewalski horse. Probes on the array represented 18,763 autosomal and X-linked genes, and intergenic, sub-telomeric and chrY sequences. We identified 258 CNV regions (CNVRs) across all autosomes, chrX and chrUn, but not in chrY. CNVs comprised 1.3% of the horse genome with chr12 being most enriched. American Miniature horses had the highest and American Quarter Horses the lowest number of CNVs in relation to Thoroughbred reference. The Przewalski horse was similar to native ponies and draft breeds. The majority of CNVRs involved genes, while 20% were located in intergenic regions. Similar to previous studies in horses and other mammals, molecular functions of CNV-associated genes were predominantly in sensory perception, immunity and reproduction. The findings were integrated with previous studies to generate a composite genome-wide dataset of 1476 CNVRs. Of these, 301 CNVRs were shared between studies, while 1174 were novel and require further validation. Integrated data revealed that to date, 41 out of over 400 breeds of the domestic horse have been analyzed for CNVs, of which 11 new breeds were added in this study. Finally, the composite CNV dataset was applied in a pilot study for the discovery of CNVs in 6 horses with XY disorders of sexual development. A homozygous deletion involving AKR1C gene cluster in chr29 in two affected horses was considered possibly causative because of the known role of AKR1C genes in testicular androgen synthesis and sexual development. While the findings improve and integrate the knowledge of CNVs in horses, they also show that for effective discovery of variants of biomedical importance, more breeds and individuals need to be analyzed using comparable methodological approaches.

Genomes of individuals in a species vary in many ways, one of which is DNA copy number variation (CNV). This includes deletions, duplications, and complex rearrangements typically larger than 50 base-pairs. CNVs are part of normal genetic variation contributing to phenotypic diversity but can also be pathogenic and associated with diseases and disorders. In order to distinguish between the two, detailed knowledge about CNVs in the species of interest is needed. Here we studied the genomes of 38 normal horses of 16 diverse breeds, and identified 258 CNV regions. We integrated our findings with previously published horse CNVs and generated a composite dataset of ∼1400 CNVRs. Despite this large number, our analysis shows that CNV research in horses needs further improvement because the current data are based on 10% of horse breeds and that most CNVRs are study-specific and require validation. Finally, we analyzed CNVs in horses with disorders of sexual development and found in two male pseudo-hermaphrodites a large deletion disrupting a group of genes involved in sex hormone metabolism and sexual differentiation. The findings underline the possible role of CNVs in complex disorders such as development and reproduction.

Repetitive Sequence and Sex Chromosome Evolution in Vertebrates

Sex chromosomes are the most dynamic entity in any genome having unique morphology, gene content, and evolution. They have evolved multiple times and independently throughout vertebrate evolution. One of the major genomic changes that pertain to sex chromosomes involves the amplification of common repeats. It is hypothesized that such amplification of repeats facilitates the suppression of recombination, leading to the evolution of heteromorphic sex chromosomes through genetic degradation of Y or W chromosomes. Although contrasting evidence is available, it is clear that amplification of simple repetitive sequences played a major role in the evolution of Y and W chromosomes in vertebrates. In this review, we present a brief overview of the repetitive DNA classes that accumulated during sex chromosome evolution, mainly focusing on vertebrates, and discuss their possible role and potential function in this process.

Independent stratum formation on the avian sex chromosomes reveals inter-chromosomal gene conversion and predominance of purifying selection on the W chromosome

We used a comparative approach spanning three species and 90 million years to study the evolutionary history of the avian sex chromosomes. Using whole transcriptomes, we assembled the largest cross-species dataset of W-linked coding content to date. Our results show that recombination suppression in large portions of the avian sex chromosomes has evolved independently, and that long-term sex chromosome divergence is consistent with repeated and independent inversions spreading progressively to restrict recombination. In contrast, over short-term periods we observe heterogeneous and locus-specific divergence. We also uncover four instances of gene conversion between both highly diverged and recently evolved gametologs, suggesting a complex mosaic of recombination suppression across the sex chromosomes. Lastly, evidence from 16 gametologs reveal that the W chromosome is evolving with a significant contribution of purifying selection, consistent with previous findings that W-linked genes play an important role in encoding sex-specific fitness.

Switching on sex: transcriptional regulation of the testis-determining gene Sry

Mammalian sex determination hinges on the development of ovaries or testes, with testis fate being triggered by the expression of the transcription factor sex-determining region Y (Sry). Reduced or delayed Sry expression impairs testis development, highlighting the importance of its accurate spatiotemporal regulation and implying a potential role for SRY dysregulation in human intersex disorders. Several epigenetic modifiers, transcription factors and kinases are implicated in regulating Sry transcription, but it remains unclear whether or how this farrago of factors acts co-ordinately. Here we review our current understanding of Sry regulation and provide a model that assembles all known regulators into three modules, each converging on a single transcription factor that binds to the Sry promoter. We also discuss potential future avenues for discovering the cis-elements and trans-factors required for Sry regulation.

Transmissible [corrected] dog cancer genome reveals the origin and history of an ancient cell lineage

Canine transmissible venereal tumor (CTVT) is the oldest known somatic cell lineage. It is a transmissible cancer that propagates naturally in dogs. We sequenced the genomes of two CTVT tumors and found that CTVT has acquired 1.9 million somatic substitution mutations and bears evidence of exposure to ultraviolet light. CTVT is remarkably stable and lacks subclonal heterogeneity despite thousands of rearrangements, copy-number changes, and retrotransposon insertions. More than 10,000 genes carry nonsynonymous variants, and 646 genes have been lost. CTVT first arose in a dog with low genomic heterozygosity that may have lived about 11,000 years ago. The cancer spawned by this individual dispersed across continents about 500 years ago. Our results provide a genetic identikit of an ancient dog and demonstrate the robustness of mammalian somatic cells to survive for millennia despite a massive mutation burden.