Comparative gene mapping in the domestic cat (Felis catus)

Exploiting genomic synteny in Felidae: cross-species genome alignments and SNV discovery can aid conservation management

Background

While recent advances in genomics has enabled vast improvements in the quantification of genome-wide diversity and the identification of adaptive and deleterious alleles in model species, wildlife and non-model species have largely not reaped the same benefits. This has been attributed to the resources and infrastructure required to develop essential genomic datasets such as reference genomes. In the absence of a high-quality reference genome, cross-species alignments can provide reliable, cost-effective methods for single nucleotide variant (SNV) discovery. Here, we demonstrated the utility of cross-species genome alignment methods in gaining insights into population structure and functional genomic features in cheetah (Acinonyx jubatas), snow leopard (Panthera uncia) and Sumatran tiger (Panthera tigris sumatrae), relative to the domestic cat (Felis catus).

Results

Alignment of big cats to the domestic cat reference assembly yielded nearly complete sequence coverage of the reference genome. From this, 38,839,061 variants in cheetah, 15,504,143 in snow leopard and 13,414,953 in Sumatran tiger were discovered and annotated. This method was able to delineate population structure but limited in its ability to adequately detect rare variants. Enrichment analysis of fixed and species-specific SNVs revealed insights into adaptive traits, evolutionary history and the pathogenesis of heritable diseases.

Conclusions

The high degree of synteny among felid genomes enabled the successful application of the domestic cat reference in high-quality SNV detection. The datasets presented here provide a useful resource for future studies into population dynamics, evolutionary history and genetic and disease management of big cats. This cross-species method of variant discovery provides genomic context for identifying annotated gene regions essential to understanding adaptive and deleterious variants that can improve conservation outcomes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07899-2.

Applications and efficiencies of the first cat 63K DNA array

The development of high throughput SNP genotyping technologies has improved the genetic dissection of simple and complex traits in many species including cats. The properties of feline 62,897 SNPs Illumina Infinium iSelect DNA array are described using a dataset of over 2,000 feline samples, the most extensive to date, representing 41 cat breeds, a random bred population, and four wild felid species. Accuracy and efficiency of the array’s genotypes and its utility in performing population-based analyses were evaluated. Average marker distance across the array was 37,741 Kb, and across the dataset, only 1% (625) of the markers exhibited poor genotyping and only 0.35% (221) showed Mendelian errors. Marker polymorphism varied across cat breeds and the average minor allele frequency (MAF) of all markers across domestic cats was 0.21. Population structure analysis confirmed a Western to Eastern structural continuum of cat breeds. Genome-wide linkage disequilibrium ranged from 50–1,500 Kb for domestic cats and 750 Kb for European wildcats (Felis silvestris silvestris). Array use in trait association mapping was investigated under different modes of inheritance, selection and population sizes. The efficient array design and cat genotype dataset continues to advance the understanding of cat breeds and will support monogenic health studies across feline breeds and populations.

Aristaless-Like Homeobox protein 1 (ALX1) variant associated with craniofacial structure and frontonasal dysplasia in Burmese cats

Frontonasal dysplasia (FND) can have severe presentations that are medically and socially debilitating. Several genes are implicated in FND conditions, including Aristaless-Like Homeobox 1 (ALX1), which is associated with FND3. Breeds of cats are selected and bred for extremes in craniofacial morphologies. In particular, a lineage of Burmese cats with severe brachycephyla is extremely popular and is termed Contemporary Burmese. Genetic studies demonstrated that the brachycephyla of the Contemporary Burmese is a simple co-dominant trait, however, the homozygous cats have a severe craniofacial defect that is incompatible with life. The craniofacial defect of the Burmese was genetically analyzed over a 20 year period, using various genetic analysis techniques. Family-based linkage analysis localized the trait to cat chromosome B4. Genome-wide association studies and other genetic analyses of SNP data refined a critical region. Sequence analysis identified a 12bp in frame deletion in ALX1, c.496delCTCTCAGGACTG, which is 100% concordant with the craniofacial defect and not found in cats not related to the Contemporary Burmese.

Investigation of inherited diseases in cats

Practical relevance:

The health of the cat mirrors a complex interaction between its environment (nurture) and its genetics (nature). To date, over 70 genetic mutations (variants) have been defined in the cat; many involve diseases, structural anomalies, coat color and texture, including numerous that are clinically relevant. This trend will continue as more of the feline genome is deciphered. Genetic testing, and eventually whole-genome sequencing, should become routine diagnostic tools in feline healthcare within the foreseeable future.

Global importance:

Cat breeds have dispersed around the world. Thus, feline medicine clinicians should be aware of breeds common to their region and common mutations found within those regional populations. Random-bred populations of domestic cats can also have defined genetic characteristics and mutations, which are equally worthy of understanding by feline medicine clinicians.

Outline:

This article reviews the chronology and evolution of genetic and genomic tools pertinent to feline medicine. Possible strategies for mapping genetic traits and defects, and how these impact on feline health, are also discussed. The focus is on three historical periods: (1) research conducted before the availability of the cat genome; (2) research performed immediately after the availability of sequences of the cat genome; and (3) current research that goes beyond one cat genome and utilizes the genome sequences of many cats.

Evidence base:

The data presented are extracted from peer-reviewed publications pertaining to mutation identification, and relevant articles concerning heritable traits and/or diseases. The authors draw upon their personal experience and expertise in feline genetics.

Molecular cloning of ion channels in Felis catus that are related to periodic paralyses in man: a contribution to the understanding of the genetic susceptibility to feline neck ventroflexion and paralysis

Neck ventroflexion in cats has different causes; however, the most common is the hypokalemia associated with flaccid paralysis secondary to chronic renal failure. In humans, the most common causes of acute flaccid paralysis are hypokalemia precipitated by thyrotoxicosis and familial forms linked to mutations in sodium, potassium, and calcium channel genes. Here, we describe the sequencing and analysis of skeletal muscle ion channels in Felis catus that could be related to periodic paralyses in humans, contributing to the understanding of the genetic susceptibility to feline neck ventroflexion and paralysis. We studied genomic DNA from eleven cats, including five animals that were hyperthyroid with hypokalemia, although only one presented with muscle weakness, and six healthy control domestic cats. We identified the ion channel ortholog genes KCNJ2, KCNJ12, KCNJ14, CACNA1S and SCN4A in the Felis catus genome, together with several polymorphic variants. Upon comparative alignment with other genomes, we found that Felis catus provides evidence for a high genomic conservation of ion channel sequences. Although we hypothesized that neck ventroflexion in cats could be associated with a thyrotoxic or familial periodic paralysis channel mutation, we did not identify any previously detected human channel mutation in the hyperthyroid cat presenting hypokalemia. However, based on the small number of affected cats in this study, we cannot yet rule out this molecular mechanism. Notwithstanding, hyperthyroidism should still be considered as a differential diagnosis in hypokalemic feline paralysis.

First WNK4-hypokalemia animal model identified by genome-wide association in Burmese cats

Burmese is an old and popular cat breed, however, several health concerns, such as hypokalemia and a craniofacial defect, are prevalent, endangering the general health of the breed. Hypokalemia, a subnormal serum potassium ion concentration ([K⁺]), most often occurs as a secondary problem but can occur as a primary problem, such as hypokalaemic periodic paralysis in humans, and as feline hypokalaemic periodic polymyopathy primarily in Burmese. The most characteristic clinical sign of hypokalemia in Burmese is a skeletal muscle weakness that is frequently episodic in nature, either generalized, or sometimes localized to the cervical and thoracic limb girdle muscles. Burmese hypokalemia is suspected to be a single locus autosomal recessive trait. A genome wide case-control study using the illumina Infinium Feline 63K iSelect DNA array was performed using 35 cases and 25 controls from the Burmese breed that identified a locus on chromosome E1 associated with hypokalemia. Within approximately 1.2 Mb of the highest associated SNP, two candidate genes were identified, KCNH4 and WNK4. Direct sequencing of the genes revealed a nonsense mutation, producing a premature stop codon within WNK4 (c.2899C>T), leading to a truncated protein that lacks the C-terminal coiled-coil domain and the highly conserved Akt1/SGK phosphorylation site. All cases were homozygous for the mutation. Although the exact mechanism causing hypokalemia has not been determined, extrapolation from the homologous human and mouse genes suggests the mechanism may involve a potassium-losing nephropathy. A genetic test to screen for the genetic defect within the active breeding population has been developed, which should lead to eradication of the mutation and improved general health within the breed. Moreover, the identified mutation may help clarify the role of the protein in K⁺ regulation and the cat represents the first animal model for WNK4-associated hypokalemia.

State of cat genomics

Our knowledge of cat family biology was recently expanded to include a genomics perspective with the completion of a draft whole genome sequence of an Abyssinian cat. The utility of the new genome information has been demonstrated by applications ranging from disease gene discovery and comparative genomics to species conservation. Patterns of genomic organization among cats and inbred domestic cat breeds have illuminated our view of domestication, revealing linkage disequilibrium tracks consequent of breed formation, defining chromosome exchanges that punctuated major lineages of mammals and suggesting ancestral continental migration events that led to 37 modern species of Felidae. We review these recent advances here. As the genome resources develop, the cat is poised to make a major contribution to many areas in genetics and biology.

Development of goldfish macrophages in vitro

Over 100 years after the first description of macrophages by Metchnikoff, there are still questions as to the mechanisms leading to the heterogeneity of their lineage. Current views are based on the mononuclear phagocyte system (MPS) theory, where all mammalian macrophages are derived from circulating blood monocytes and ultimately from hematopoietic stem cells in the bone marrow. Our studies on the regulation of fish macrophage development, suggested that teleosts have alternate pathways of monopoiesis, which undoubtedly contribute to macrophage heterogeneity in the goldfish. Macrophage heterogeneity has been attributed to a network of positive and negative regulators of macrophage development, including soluble mediators known as colony-stimulating factors of which two (M-CSF and GM-CSF) promote formation and growth of mature macrophages. In contrast to our knowledge of CSFs and their receptors in mammals, there is no published information about fish macrophage CSFs. Since fish macrophages generate their own growth factors, it is reasonable to assume that pathways of fish macrophage development and hematopoiesis may be distinct from those of mammalian macrophages. More importantly, the presence of fish progenitor/stem cells and developing macrophages in long-term cultures, allowed us to address pathways of macrophage differentiation, which could not be addressed in mammalian macrophage culture systems. Characterization of primary kidney macrophage (PKM) cultures from goldfish hematopoietic tissues (kidney) indicated that three distinct subpopulations developed in response to endogenous macrophage growth factors. These macrophage subpopulations expressed several differentiation markers, including the hematopoietic stem cell antigen AC133, c-kit, granulin, CD63, macrosialin, c/EBPbeta, legumain, and the colony-stimulating factor receptor-1 (CSF-1R). In the goldfish, there appeared to be a stringent control between those early progenitors that self-renewed, and those that were recruited into the maturation pathways. We report that upon commitment, goldfish macrophages developed through two distinct differentiation pathways: one consistent with the "classical" pathway (MPS) of macrophage development (progenitors-->monocytes-->mature macrophages), and an "alternate" pathway (AP-macrophages) where mature macrophages appeared to rapidly develop from early progenitors in the absence of an intermediate monocyte stage. AP-macrophages represent a unique subset of spontaneously growing cells. Their self-renewal was promoted by endogenous macrophage growth factors (MGF), and effectively controlled by a novel soluble form of the CSF-1R (sCSF-1R). The discovery of sCSF-1R in the goldfish highlights the inherent complexity in the hematopoietic regulatory machinery of teleosts.

Microscopy and image analysis

This unit provides an overview of light microscopy, including objectives, light sources, filters, film, and color photography for fluorescence microscopy and fluorescence in situ hybridization (FISH). Computerized image-analysis systems currently used in clinical cytogenetics are also discussed.

Tyrosinase and tyrosinase related protein 1 alleles specify domestic cat coat color phenotypes of the albino and brown loci

The genes encoding enzymes of the tyrosinase family are strong candidates for coat color variation in mammals. To investigate their influence in domestic cat coat color, we determined the complete nucleotide coding sequence of the domestic cat genes tyrosinase (TYR)--a plausible candidate gene for the albino (C) locus, and tyrosinase related protein 1 (TYRP1)--a candidate gene for the brown (B) locus. Sequence variants between individuals exhibiting variation in pigmentation were submitted to association studies. In TYR, two nonsynonymous substitutions encoding TYR-G301R and TYR-G227W were associated with the siamese and burmese phenotypes of the albino locus, respectively. TYRP1 was mapped on chromosome D4 within 5 cM of a highly polymorphic microsatellite, previously found to be fixed in a cat breed selected for the chocolate (b) allele of the B locus, which reinforced TYRP1 as a candidate gene for the B locus in the domestic cat. Two DNA polymorphisms, one leading to a TYRP1-A3G substitution in the signal peptide and another to an in-frame insertion TYRP1-421ins17/18 caused by a donor splice site mutation in intron 6, were associated with the chocolate (b) allele. A premature UAG stop codon at position 100 of TYRP1 was associated with a second allele of the B locus, cinnamon (b(l)). The results provide very strong evidence that the specific nucleotide variants of feline TYR (chromosome D1) are causative of the siamese (c(s)) and burmese (c(b)) alleles of the albino locus, as well as nucleotide variants of TYRP1 (chromosome D4) as specifying the chocolate (b) and cinnamon (b(l)) alleles of the B locus.

Development and Characterization of CATS Markers for Genetic Linkage Mapping in the House Musk Shrew, Suncus murinus

To serve as an initial step in developing an ideal genetic marker map for the house musk shrew, Suncus murinus, 318 comparative anchor tagged sequence (CATS) primer pairs were assessed for polymorphism ascertainment and linkage mapping. Of the 112 (35.2%) CATS primer pairs that were successfully amplified by PCR in the shrew, 18 (16.1%) showed polymorphism between two mutant strains, BAN-kc, oeb and WZ. Linkage analysis of the polymorphic CATS markers and three visible mutant genes, kc, oeb and wz, genotyped in a 77 F2 mapping panel from a cross of the two mutant strains, assigned wz and five CATS markers into three linkage groups. Sequence analysis revealed that two (ADA and TXN) out of nine CATS amplified sequences had a total of six deletions of varying sizes and 17 single nucleotide polymorphisms (SNPs). BLAST search identified three CATS (ADA, CYP1A2, and TXN) products matching the genes from which they were originally designed, while the remaining six markers could not be identified. Together with the use of the detected SNPs as genetic markers, the five CATS markers linkage mapped in this species will serve as anchors in establishing the first framework map for locating loci affecting all heritable qualitative and quantitative traits in the musk shrew.

The feline major histocompatibility complex is rearranged by an inversion with a breakpoint in the distal class I region

In order to determine the genomic organization of the major histocompatibility complex (MHC) of the domestic cat (Felis catus), DNA probes for 61 markers were designed from human MHC reference sequences and used to construct feline MHC BAC contig map spanning ARE1 in the class II region to the olfactory receptor complex in the extended class I region. Selected BAC clones were then used to identify feline-specific probes for the three regions of the mammalian MHC (class II-class III-class I) for radiation hybrid mapping and fluorescent in situ hybridization to refine the organization of the domestic cat MHC. The results not only confirmed that the p-arm of domestic cat B2 is inverted relative to human Chromosome 6, but also demonstrated that one inversion breakpoint localized to the distal segment of the MHC class I between TRIM39 and TRIM26. The inversion thus disjoined the approximately 2.85 Mb of MHC containing class II-class III-class I (proximal region) from the approximately 0.50 Mb of MHC class I/extended class I region, such that TRIM39 is adjacent to the Chromosome B2 centromere and TRIM26 is adjacent to the B2 telomere in the domestic cat.

The Feline Genome Project

The compilation of a dense gene map and eventually a whole genome sequence (WGS) of the domestic cat holds considerable value for human genome annotation, for veterinary medicine, and for insight into the evolution of genome organization among mammals. Human association and veterinary studies of the cat, its domestic breeds, and its charismatic wild relatives of the family Felidae have rendered the species a powerful model for human hereditary diseases, for infectious disease agents, for adaptive evolutionary divergence, for conservation genetics, and for forensic applications. Here we review the advantages, rationale, and present strategy of a feline genome project, and we describe the disease models, comparative genomics, and biological applications posed by the full resolution of the cat's genome.

Evolution of mammalian genome organization inferred from comparative gene mapping

Comparative genome analyses, including chromosome painting in over 40 diverse mammalian species, ordered gene maps from several representatives of different mammalian and vertebrate orders, and large-scale sequencing of the human and mouse genomes are beginning to provide insight into the rates and patterns of chromosomal evolution on a whole-genome scale, as well as into the forces that have sculpted the genomes of extant mammalian species.

Chromosome painting in farm, pet and wild animal species

Among the advanced karyotype analysis approaches embraced by animal cytogenetics during the past decade, chromosome painting has had the greatest impact. Generation of chromosome specific paints is considered pivotal to his development. Additionally, ability to use these paints across species (referred to as Zoo-FISH or comparative painting) is undisputedly the most important breakthrough that has contributed to our ability to compare karyotypes of a wide range of evolutionarily highly diverged chromosome painting, and makes them aware of the tools/resources available to carry out this research in a variety of animal species. An overview of the current status of comparative chromosome painting results across closely as well as distantly related species is presented. Findings from different studies show how some segmental syntenies are more conserved as compared to others. The comparisons provide insight into the likely constitution of a vertebrate/mammalian ancestral karyotype and help understand some of the intricacies about karyotype evolution. Importance of comparative painting in setting the stage for rapid development of gene maps in a number of economically important species is elaborated.

A radiation hybrid map of the cat genome: implications for comparative mapping

Ordered gene maps of mammalian species are becoming increasingly valued in assigning gene variants to function in human and animal models, as well as recapitulating the natural history of genome organization. To extend this power to the domestic cat, a radiation hybrid (RH) map of the cat was constructed integrating 424 Type I-coding genes with 176 microsatellite markers, providing coverage over all 20 feline chromosomes. Alignment of parallel RH maps of human and cat reveal 100 conserved segments ordered (CSOs) between the species, nearly three times the number observed with reciprocal chromosome painting analyses. The observed number is equivalent to theoretical predictions of the number of conserved segments to be found between cat and human, implying that 300-400 Type I gene markers is sufficient to reveal nearly all conserved segments for species that exhibit the most frequently observed "slow" rate of genome reorganization. The cat-human RH map comparisons provide a new genomic tool for comparative gene mapping in the cat and related Felidae, and provide confirmation that the cat genome organization is remarkably conserved compared with human. These data demonstrate that ordered RH-based gene maps provide the most precise assessment of comparing genomes, short of contig construction or full-sequence determination.

Extensive conservation of sex chromosome organization between cat and human revealed by parallel radiation hybrid mapping

A radiation hybrid (RH)-derived physical map of 25 markers on the feline X chromosome (including 19 Type I coding loci and 6 Type II microsatellite markers) was compared to homologous marker order on the human and mouse X chromosome maps. Complete conservation of synteny and marker order was observed between feline and human X chromosomes, whereas the same markers identified a minimum of seven rearranged syntenic segments between mouse and cat/human X chromosome marker order. Within the blocks, the feline, human, and mouse marker order was strongly conserved. Similarly, Y chromosome locus order was remarkably conserved between cat and human Y chromosomes, with only one marker (SMCY) position rearranged between the species. Tight linkage and a conserved gene order for a segment encoding three genes, DFFRY-DBY-UTY in human, mouse, and cat Y chromosomes, coupled with demonstrated deletion effects of these genes on reproductive impairment in both human and mouse, implicates the region as critical for Y-mediated sperm production.

A comparative gene map of the horse (Equus caballus)

A comparative gene map of the horse genome composed of 127 loci was assembled based on the new assignment of 68 equine type I loci and on data published previously. PCR primers based on consensus gene sequences conserved across mammalian species were used to amplify markers for assigning 68 equine type I loci to 27 horse synteny groups established previously with a horse-mouse somatic cell hybrid panel (SCHP, UC Davis). This increased the number of coding genes mapped to the horse genome by over 2-fold and allowed refinements of the comparative mapping data available for this species. In conjunction with 57 previous assignments of type I loci to the horse genome map, these data have allowed us to confirm the assignment of 24 equine synteny groups to their respective chromosomes, to provisionally assign nine synteny groups to chromosomes, and to further refine the genetic composition established with Zoo-FISH of two horse chromosomes. The equine type I markers developed in this study provide an important resource for the future development of the horse linkage and physical genome maps.

Construction and characterization of an eightfold redundant dog genomic bacterial artificial chromosome library

A large insert canine genomic bacterial artificial chromosome (BAC) library was built from a Doberman pinscher. Approximately 166,000 clones were gridded on nine high-density hybridization filters. Insert analysis of randomly selected clones indicated a mean insert size of 155 kb and predicted 8.1 coverage of the canine genome. Two percent of the clones were nonrecombinant. Chromosomal fluorescence in situ hybridization studies of 60 BAC clones indicated no chimerism. The library was hybridized with dog PCR products representing eight genes (ADA, TNFA, GCA, MYB, HOXA, GUSB, THY1, and TOP1). The resulting positive clones were characterized and shown to be compatible with an eightfold redundant library.

Development of a feline whole genome radiation hybrid panel and comparative mapping of human chromosome 12 and 22 loci

A 5000-rad whole genome radiation hybrid panel is described for the domestic cat, derived from irradiated male feline fibroblasts fused to a recipient hamster cell line. A panel of 93 cell lines has an estimated retention frequency of 0.39 (range 0.13-0.71) based upon PCR typing of 54 feline markers. To test the panel's utility, we determined the order of 16 Type I (coding gene) loci, 14 Type II (microsatellite) loci, and 1 endogenous retroviral element on feline chromosomes B4 and D3. Assessment of marker order derived from the RH panel was compared to assignments of the same loci using interspecies backcross mapping data, human homologue positions, and human-cat chromosome painting homologies. Assessment of concordant and discordant marker order for these loci provides improved resolution into the evolution of subchromosomal genome organizations and the methods to track them in these species.

A genetic linkage map of microsatellites in the domestic cat (Felis catus)

Of the nonprimate mammalian species with developing comparative gene maps, the feline gene map (Felis catus, Order Carnivora, 2N = 38) displays the highest level of syntenic conservation with humans, with as few as 10 translocation exchanges discriminating the human and feline genome organization. To extend this model, a genetic linkage map of microsatellite loci in the feline genome has been constructed including 246 autosomal and 7 X-linked loci. Two hundred thirty-five dinucleotide (dC. dA)n. (dG. dT)n and 18 tetranucleotide repeat loci were identified and genotyped in a two-family, 108-member multigeneration interspecies backcross pedigree between the domestic cat (F. catus) and the Asian leopard cat (Prionailurus bengalensis). Two hundred twenty-nine loci were linked to at least one other marker with a lod score >/=3.0, identifying 34 linkage groups. Representative markers from each linkage group were assigned to specific cat chromosomes by somatic cell hybrid analysis, resulting in chromosomal assignments to 16 of the 19 feline chromosomes. Genome coverage spans approximately 2900 cM, and we estimate a genetic length for the sex-averaged map as 3300 cM. The map has an average intragroup intermarker spacing of 11 cM and provides a valuable resource for mapping phenotypic variation in the species and relating it to gene maps of other mammals, including human.

Emerging patterns of comparative genome organization in some mammalian species as revealed by Zoo-FISH

Although gene maps for a variety of evolutionarily diverged mammalian species have expanded rapidly during the past few years, until recently it has been difficult to precisely define chromosomal segments that are homologous between species. A solution to this problem has come from the development of Zoo-FISH, also known as cross-species chromosome painting. The use of Zoo-FISH to identify regions of chromosomal homology has allowed the transfer of information from map-rich species such as human and mouse to a wide variety of other species. From a Zoo-FISH analysis spanning four mammalian orders (Primates, Artiodactyla, Carnivora, and Perissodactyla), and involving eight species (human, pig, cattle, Indian muntjac, cat, American mink, harbor seal, and horse), three distinct classes of synteny conservation have been designated: (1) conservation of whole chromosome synteny, (2) conservation of large chromosomal blocks, and (3) conservation of neighboring segment combinations. This analysis has also made it possible to identify a set of chromosome segments (based on human chromosome equivalents) that probably made up the karyotype of the common ancestor of the four orders. This approach provides a basis for developing a picture of the ancestral mammalian karyotype, but a full understanding will depend on studies encompassing more diverse combinations of mammalian orders.