Sequencing and comparative genomic analysis of 1227 Felis catus cDNA sequences enriched for developmental, clinical and nutritional phenotypes
© Irizarry et al; licensee BioMed Central Ltd. 2012
Received: 20 April 2011
Accepted: 18 January 2012
Published: 18 January 2012
The feline genome is valuable to the veterinary and model organism genomics communities because the cat is an obligate carnivore and a model for endangered felids. The initial public release of the Felis catus genome assembly provided a framework for investigating the genomic basis of feline biology. However, the entire set of protein coding genes has not been elucidated.
We identified and characterized 1227 protein coding feline sequences, of which 913 map to public sequences and 314 are novel. These sequences have been deposited into NCBI's genbank database and complement public genomic resources by providing additional protein coding sequences that fill in some of the gaps in the feline genome assembly. Through functional and comparative genomic analyses, we gained an understanding of the role of these sequences in feline development, nutrition and health. Specifically, we identified 104 orthologs of human genes associated with Mendelian disorders. We detected negative selection within sequences with gene ontology annotations associated with intracellular trafficking, cytoskeleton and muscle functions. We detected relatively less negative selection on protein sequences encoding extracellular networks, apoptotic pathways and mitochondrial gene ontology annotations. Additionally, we characterized feline cDNA sequences that have mouse orthologs associated with clinical, nutritional and developmental phenotypes. Together, this analysis provides an overview of the value of our cDNA sequences and enhances our understanding of how the feline genome is similar to, and different from other mammalian genomes.
The cDNA sequences reported here expand existing feline genomic resources by providing high-quality sequences annotated with comparative genomic information providing functional, clinical, nutritional and orthologous gene information.
KeywordsFeline bioinformatics comparative genomics cDNA annotation gene ontology OMIM ortholog nutrition phenotype
The domestic cat, Felis catus, is a member of the family Felidae and represents the Feliformia branch of the order Carnivora . The domestic cat is an important companion animal and veterinary species. There are roughly 82 million companion cats living in more than 35 million US households . The domestic cat also has substantial value as a model organism for comparative mammalian genomics because it is an obligate carnivore , unlike the dog which is an omnivore . Additionally, the domestic cat is an important model organism for Felidae because of its close phylogenetic relationship to the wildcat (Felis silvestris), the sand cat (Felis margarita), the black-footed cat (Felis nigripes) and the jungle cat (Felis chaus). It can also serve as a model for more distantly related felid species including pumas such as the Cheetah (Acinonyx jubatus), lynx species, ocelots [5, 6], and members of panthera including the lion (Panthera leo), the tiger (Panthera tigris), and snow leopard (Uncia uncia) . A major goal of feline genomics is to identify and decode both cat-specific biology as well as conserved mammalian biology. The identification of feline-specific biochemistry and physiology is required in order to better understand the unique nutritional and veterinary needs of cats and to enhance the wellness of domestic cats as well as the health and management of captive felid species.
A number of cat-specific biological adaptations have been described to date. Cats exhibit a variety of evolutionary adaptations thought to be associated with their predatory behaviour and obligate carnivore status. For example, domestic cats exhibit distinct distal forelimb anatomical adaptations associated with predation [8, 9], as well as sensory adaptations in both sound perception [10–12] and visual acuity [13, 14]. At a molecular level, cats exhibit differences in the regulation of sugar transporters  resulting in lower liver glucose transporter activity  and differences in carbohydrate metabolism compared to omnivores . Because the carnivore diet is relatively high in amino acid content, adult cats maintain blood glucose levels from gluconeogenesis of glucogenic amino acids, lactic acid and glycerol . Compared to omnivorous mammals, in which gluconeogenesis occurs in the post absorptive state, cats exhibit the greatest extent of gluconeogenesis right after a meal during the absorptive state .
Amino acid biosynthesis and deficiency has been relatively well studied in domestic cats. Cats have dietary requirements for the amino acids taurine , arginine , cysteine and, methionine . Arginine deficiency in cats has been associated with rapid onset of hyperammonemia characterized by severe signs of ammonia toxicity . The sulphur containing amino acids cysteine and methionine are normally present in high amounts in animal flesh and are required for normal feline development [24, 25]. The beta-amino sulfonic acid taurine is required in cats because, unlike many other species which can conjugate bile acids to either glycine or taurine for secretion of bile salts into bile, cats can only use taurine. Unlike dogs, cats have evolved limited capacity to synthesize taurine , subsequently, taurine deficiency in cats is associated with abnormal cardiac , immune , neurological , platelet , reproductive  and retinal  dysfunctions. The recent description of the taurine transporter knock out mouse underscores the biological roles of taurine in mammals .
Although many aspects of feline-specific biology have been elucidated to date, bioinformatics methods and comparative genomics approaches can provide a mechanism for producing a number of plausible and useful biological hypotheses from feline cDNA sequences.
The 2007 release of the feline genome  marked the beginning of the feline genomics era, which was followed by the identification of close to 1 million single nucleotide polymorphisms across cat breeds  which further extends the repertoire of genomic tools for investigating the genomic basis of feline phenotypes. In this paper, we describe the sequencing of additional feline cDNA sequences and demonstrate the utility of employing comparative genomics methods to investigate, not only the roles of these cDNA sequences, but the extent to which these feline sequences diverge from other mammalian orthologous sequences.
Our working hypothesis is that conservation among human, mouse, dog and cat orthologs underscores conserved mammalian biology while feline sequence divergence among mammalian orthologs provides potential insight into cat-specific biology. Specifically, we employ a computational comparative gene expression analysis to map the cDNA sequences to anatomical information, developmental timelines, cells and pathology terms. Additionally, we utilize the gene ontology annotation, in combination with measures of synonymous and non-synonymous differences in orthologous protein sequences, to better understand which of the cDNA sequences are likely to represent conserved mammalian biology and which are more likely to represent feline-specific biology. We organize these results into biological processes, cellular localization and molecular function in order to more easily interpret the results. Finally, we map these feline cDNA sequences to orthologs in other species in order to identify (1) phenotypes, (2) biochemical pathways and (3) human diseases in an attempt to better understand the roles of these cDNA sequences in feline development, nutrition and disease.
Sequencing and Orthologue Identification
For the set of 2081 cDNA sequences, the shortest and longest sequences were 353 and 4750 nucleotides respectively. The average nucleotide length was 1349 nucleotides with a standard deviation of 567 nucleotides. The 2081 protein sequence set exhibited a shortest and longest sequence of 41 and 1128 amino acids respectively. The average protein sequence length was 279 amino acids with a standard deviation of 149 amino acids.
This set of sequences was used to blast against the set of known human cDNA and protein sequences to identify the best human match (see Figure 1). Additionally, these 2081 cDNA sequences were blasted against known and ab initio feline cDNA and protein sequences from ensemble  to identify sequences for which public feline sequence data exists. Subsequently, these sequences were aligned using a global alignment algorithm to remove sequences for which the best blast hit represented only local homology. After manual review of all of the global nucleotide and protein alignments, a set of 1227 non-redundant feline sequences were selected as high confidence, high quality feline sequences. Within the set of 1227 sequences, 913 known sequences and 314 novel sequences were identified for which 914 were successfully mapped to their corresponding dog, human and mouse orthologs. Although additional non-redundant feline cDNA sequences we identified mapped to three or fewer orthologs across the four species, we limited our subsequent analysis to only those sequences for which all three non-feline species orthologs were confidently identified. This decision was made to ensure that our functional and comparative analysis would include only feline cDNA sequences for which dog, mouse and human orthologs were identified. Of the 914 orthologous sequence set, 844 sequences corresponded to known feline sequences and 70 corresponded to novel sequences (see Figure 1). Additional file 1, Table S1 contains the complete set of 1227 non-redundant nucleotide and protein sequences. The complete set of 914 orthologous sequences is listed in Additional file 2, Table S2 along with the designation of known or novel and the corresponding ensembl gene, transcript and protein identifiers for the dog, human and mouse orthologs.
Comparative Gene Expression Analysis
The sequences we report were obtained from extensive sequencing of 21 individual tissue cDNA libraries and 1 pooled cDNA library. It is well known that while some genes may exhibit rather narrow ranges of expression across tissues and cell types, many genes exhibit expression across numerous tissues and cell types [37, 38]. We chose to leverage the orthologous relationships among our sequences to infer gene expression patterns across a set of anatomical regions.
Inferred Anatomical Gene Expression Patterns
Number of Genes
% of Genes
head and neck
islets of Langerhans
Inferred Cell Type Gene Expression Patterns
Number of Genes
% of Genes
cardiac muscle cell
natural killer cell
retinal pigment epithelium
skeletal muscle cell
smooth muscle cell
Inferred Pathology Gene Expression Patterns
Number of Genes
% of Genes
carcinoma in situ
papillary serous carcinoma
Inferred Developmental Gene Expression Patterns
Number of Genes
% of Genes
Taken together these results suggest that the genes encoding the cDNA sequences we have identified exhibit considerably larger breadth of expression than would be suggested from the initial tissues that were sequenced. The broad extent of tissue, cell type, developmental and pathological expression annotation suggests that these sequences may include sequences underlying tissue and organ development as well as contributing to specific pathological conditions. In order to better understand the biological role of these genes we chose to combine the expression annotation with other functional and comparative annotation types.
Gene Ontology Annotation Analysis
Gene ontology (GO) annotation was performed on the feline sequences using the previously identified comparative genomics ortholog relationships. Gene ontology terms were mapped from human annotation files to feline orthologs. The initial gene ontology human molecular function annotation file contained 73,467 function annotation terms mapped to 21,956 human gene identifiers, corresponding to 3,085 unique gene ontology function terms. The cellular location gene ontology annotation file contained 975 unique terms mapped to 21,956 human genes resulting in 69,556 gene-term relationships. The biological process gene ontology annotation also contained 21,956 human gene identifiers consisting of 6518 unique gene ontology process annotation terms represented by 89,968 gene-to-GO entries.
The mapping of gene ontology functional annotation terms onto the non-redundant full length sequences resulted in 901 of our feline cDNA sequences becoming associated with 647 unique gene ontology molecular function annotation terms resulting in 3219 annotation-gene relationships. Repeating the procedure to map the cellular location annotation, we mapped 3423 gene-annotation relationships corresponding to 337 unique location annotation terms covering the set of 901 genes. Mapping the biological process annotation terms produced 4247 gene-to-GO annotations of which 1441 unique gene ontology process annotations mapped successfully to 901 genes.
Typically gene ontology annotation terms are filtered using an enrichment criterion that is calculated from a hypergeometric null model to describe the number of annotation terms one might expect to occur within a gene set of a given size and a GO annotation distribution of particular parameters. Although such an approach is necessary when attempting to determine the biological role of a gene set, such as up-regulated or down-regulated genes in a gene expression study, we did not calculate an enrichment of gene ontology terms, instead we combined the gene ontology annotation with measures of evolutionary selection using non-synonymous (dN) versus synonymous (dS) codon statistics as a means of exploring the evolutionary relationships that exist among the different gene ontology annotations across our cDNA sequences. A well accepted approach for identifying evidence of positive selection is to identify genes exhibiting significantly larger rates of non-synonymous substitutions per non-synonymous site than synonymous substitutions per synonymous site. Evidence of fixation exists when the ratio of non-synonymous substitution rate to synonymous substitution rate equals zero (dN/dS = 0).
Evidence of negative selection exists when dN/dS < 1 and evidence of positive selection exists when dN/dS > 1. We recognize that using the dN/dS value across an entire gene is an extremely conservative measure of selection, and that smaller regions within a gene may exhibit local signals of positive selection . However, we chose the conservative approach in order to minimize reporting false positives due to the possibility of sequencing errors.
Instead of considering all of the genes we identified as a single gene set, we chose to select gene subsets using SQL queries in MySQL to identify cDNA sequences sharing gene ontology annotation terms for which we calculated an average dN/dS value. From this analysis, we were able to identify annotation types exhibiting low dN/dS values, corresponding to greater levels of sequence conservation across species. We were also able to identify annotation terms that exhibited considerably higher dN/dS values indicating less negative selection in the act on some types of genes. Because we chose to employ a stringent criteria for positive selection, we did not identify genes exhibiting strong signals of positive selection, instead, we were able to identify genes and annotation types with different levels of selection pressure acting on them. Beginning with the gene ontology location annotation, an SQL query was performed such that the genes exhibiting the same location annotation terms were grouped together and the average dN/dS value was calculated for cat versus dog, cat versus human and cat versus mouse. Location annotations occurring within gene sets that exhibit extremely low dN/dS values and very low standard deviation of the dN/dS value for each species were selected as negatively selected location annotation gene sets.
Location annotation terms associated with genes exhibiting a larger dN/dS value included extracellular region (48 genes, dN/dS = 0.12), extracellular matrix (10 genes, dN/dS = 0.13), extracellular space (25 genes, dN/dS = 0.13) and integral to plasma membrane (26 genes, dN/dS = 0.14). Some of the location terms with the greatest dN/dS values were associated with the mitochondria, for example respiratory chain (8 genes, dN/dS = 0.1), mitochondrial membrane (5 genes, dN/dS = 0.12), mitochondrial matrix (6 genes, dN/dS = 0.16), mitochondrial large ribosomal subunit (2 genes, dN/dS = 0.19) and mitochondrial small subunit (4 genes, dN/dS = 0.31).
Examples of process terms relating to these themes include neuroprotection (3 genes, dN/dS = 0.13), positive regulation of apoptosis (2 genes, dN/dS = 0.13), response to estradiol stimulus (6 genes, dN/dS = 0.14), response to peptide hormone stimulus (3 genes, dN/dS = 0.14), induction of apoptosis (10 genes, dN/dS = 0.15), positive regulation of MAPKKK cascade (2 genes, dN/dS = 0.16), cell surface receptor linked signalling pathway (5 genes, dN/dS = 0.16), response to antibiotic (4 genes, dN/dS = 0.18), response to hyperoxia (2 genes, dN/dS = 0.19), activation of caspase activity (5 genes, dN/dS = 0.21), response to cAMP (3 genes, dN/dS = 0.22), positive regulation of neuron apoptosis (3 genes, dN/dS = 0.24) and blood coagulation (6 genes, dN/dS = 0.22).
Function annotation terms associated with higher dN/dS values include transcription regulator activity (11 genes, dN/dS = 0.11), G-protein coupled receptor activity (5 genes, dN/dS = 0.11), DNA directed RNA polymerase activity (5 genes, dN/dS = 0.13), signal transducer activity (13 genes, dN/dS = 0.14), promoter binding (4 genes, dN/dS = 0.13), guanyl-nucleotide exchange factor activity (6 genes, dN/dS = 0.17), phosphoprotein phosphatase inhibitor activity (3 genes, dN/dS = 0.18), GTPase activator activity (6 genes, dN/dS = 0.19), vitamin D binding (2 genes, dN/dS = 0.26) and transmembrane receptor activity (2 genes, dN/dS = 0.39).
The themes observed in this data provide insight into the inner workings of the cell and shed light on the evolutionary constraints that act on different components of the intracellular machinery. The fact that the these feline sequences include a distribution of gene products, some of which are strongly conserved across human/mouse/dog, suggests that these sequences include genes that play very important roles in critical cellular processes and correspond to conserved mammalian cellular biology. However, some genes map to protein products that have relatively less selective pressure acting on them. These gene products are also important because they represent the targets of adaptive evolution within the cell. While microtubule structure and function must be highly conserved, regulatory gene products are freer to evolve new interactions that may increase fitness of the cell. Figures 3 through 5 contain the three types of gene ontology annotation together with the average dN/dS values for genes exhibiting the same annotation types. Although this analysis of dN/dS values across our genes provided a gene level picture of our data, we wanted to investigate the large-scale pattern of dN/dS values across our cDNA sequences.
GeneGO Analysis of Orthologous Genes by dN/dS Value
In order to gain a more global view of how the feline cDNA sequences compared to other species, a set of 711 cDNA sequences having orthologs containing gene ontology annotation across dog, mouse and human were analysed to detect any non-random patterns across the genes, species and annotations. We sorted a list of 711 genes by dN/dS value and identified 3 groupings, corresponding to the top 25% of dN/dS values, the bottom 25% of dN/dS values and the middle 50% of dN/dS values. Each list was used to query the GeneGO annotation database for metabolic pathways.
Summary Statistics for GeneGO Annotation Analysis
dN/dS < 0.0149
0.0149 < dN/dS < 0.1229
dN/dS > 0.1229
Common networks between top 25% and bottom 25%
We examined the metabolic networks of these genes in GeneGO. We observed that the group with lower dN/dS values exhibited fewer numbers of amino acid type metabolic networks than the group with larger dN/dS values. Our examination of these metabolic network annotations across the groups of genes provides insight into an interesting pattern that was not apparent from the gene level gene ontology analysis described in the preceding section.
Conserved and Divergent Metabolic Pathways
(L)-leucine pathways and transport
Branched-chain amino acid metabolism
Estrone and Estradiol metabolism
N-acyl-sphingosine phosphate pathway
HETE, HPETE and Leukotriene4 metabolism
Glutamic acid pathway
Tryptophan, Phenylalanine, Methionine metabolism
L-glutamate pathways and transport
Tryptophan, Phenylalanine, Tyramine, Methionine metabolism and transport
Glycolysis, Glucogenesis and glucose transport
Glutamic acid pathways and transport
Identification of Metabolic and Biochemical Pathways
Based on the GeneGO findings, we wanted to gain further insight into the biochemical role of the feline cDNA sequences. We chose to further explore how our cDNA sequences map onto metabolic pathways by identifying a set of pathways for which at least one pathway member has been identified in the set of our orthologous cat cDNA sequences. This analysis identified ten distinct classes of biochemical pathways for which 112 feline cDNA sequences have been mapped to 75 different pathways.
The categories of pathways include amino acid metabolism, biosynthesis of secondary metabolites, carbohydrate metabolism, energy metabolism, lipid metabolism, nucleotide metabolism as well as glycan biosynthesis and metabolism, metabolism of cofactors and vitamins and xenobiotic biodegradation and metabolism.
Genes Mapped to KEGG Pathways
Number of Genes
Amino Acid Metabolism
Alanine and aspartate metabolism
Arginine and proline metabolism
Glycine, serine and threonine metabolism
Phenylalanine, tyrosine and tryptophan biosynthesis
Starch and sucrose metabolism
Urea cycle and metabolism of amino groups
Valine, leucine and isoleucine biosynthesis
Valine, leucine and isoleucine degradation
Biosynthesis of Secondary Metabolites
Alkaloid biosynthesis II
Limonene and pinene degradation
Ascorbate and aldarate metabolism
Citrate cycle (TCA cycle)
Fructose and mannose metabolism
Glyoxylate and dicarboxylate metabolism
Inositol phosphate metabolism
Nucleotide sugars metabolism
Pentose and glucuronate interconversions
Pentose phosphate pathway
Glycan Biosynthesis and Metabolism
Glycan structures - biosynthesis 1
Glycan structures - biosynthesis 2
Glycosphingolipid biosynthesis - lacto and neolacto series
Heparan sulfate biosynthesis
Keratan sulfate biosynthesis
alpha-Linolenic acid metabolism
Androgen and estrogen metabolism
Arachidonic acid metabolism
Bile acid biosynthesis
Biosynthesis of steroids
Biosynthesis of unsaturated fatty acids
Ether lipid metabolism
Fatty acid elongation in mitochondria
Fatty acid metabolism
Linoleic acid metabolism
Metabolism of Cofactors and Vitamins
Nicotinate and nicotinamide metabolism
Pantothenate and CoA biosynthesis
Porphyrin and chlorophyll metabolism
Vitamin B6 metabolism
Metabolism of Other Aminoacids
Xenobiotics Biodegradation and Metabolism
1- and 2-Methylnaphthalene degradation
Benzoate degradation via CoA ligation
Drug metabolism - cytochrome P450
Drug metabolism - other enzymes
Metabolism of xenobiotics by cytochrome P450
Comparative Phenotype Analysis
Phenotype Modules and Feline Disorders
-Ras association (RalGDS/AF-6) domain family member 1
-solute carrier family 22 (organic cation/carnitine transporter), member 5
(primary carnitine deficiency )
-cysteine-rich, angiogenic inducer, 61
(prostate cancer )
-transmembrane protein 38A
-eukaryotic translation initiation factor 2-alpha kinase 1
-snail homolog 1 (Drosophila)
(colon cancer )
-interleukin 1 receptor antagonist
-dilated dorsal aorta
-abnormal mitral valve morphology
-abnormal cardiac output
-abnormal myocardial fiber physiology
-abnormal outflow tract
-abnormal coronary artery morphology
-mitral valve dysplasia
-tricuspid valve dysplasia
-ventricular septal defects
-atrial septal defects
-dynamic subaortic stenosis
-TGFB-induced factor homeobox 2
(Alzheimer's disease )
-E binding protein 1
-potassium inwardly-rectifying channel, subfamily J, member 1
-retinol dehydrogenase 12 (all-trans/9-cis/11-cis)
(Leber's congenital amaurosis )
-arginine vasopressin receptor 1A
-peptidylprolyl cis/trans isomerase, NIMA-interacting 1
(Alzheimer's disease )
-abnormal mesoderm development
-abnormal proximal/distal developmental patterning
-abnormal rostral/caudal developmental patterning
-embryonic growth arrest
-abnormal trophoblast layer morphology
-abnormal white adipose tissue
-decreased renal glomerular filtration rate
-decreased cholesterol levels
-abnormal intestine morphology
-post natal growth retardation
-retinal neuronal layer morphology
-abnormal retinal apoptosis
-abnormal circulating corticosterone level
-abnormal adrenal gland morphology
-decreased primordial germ cell number
-seminiferous tubule degeneration
-chronic kidney disease
-polycystic kidney disease
Immune and Hematopoietic
-tumor protein D52-like 2
(childhood leukemia )
-fragile histidine triad gene
-beclin 1, autophagy related
(hyper-IgM syndrome )
-solute carrier family 35, member C1
(leukocyte adhesion deficiency II )
-linker for activation of T cell
-bridging integrator 3
-interleukin 1 receptor antagonist
-abnormal macrophage physiology
-abnormal T-cell proliferation
-abnormal B-cell proliferation
-increased susceptibility to infection
-decreased granulocyte number
-abnormal erythrocyte morphology
-decreased platelet number
-abnormal immune system biology
-increased spleen germinal cell number
-increased spleen germinal size
-abnormal class switch recombination
-abnormal somatic hypermutation frequency
-abnormal lymph node primary follicle
-abnormal leukocyte adhesion
-abnormal cytokine secretion
-abnormal interferon secretion
-increased IgE levels
-increased IgG1 levels
-increased IgM levels
- increased interferon gamma secretion
- increased interleukin 10 secretion
- increased interleukin 4 secretion
- enlarged spleen
- liver inflammation
- abnormal chemokine secretion
- abnormal macrophage recruitment
- increased susceptibility to endotoxin shock
- chronic inflammation,
- increased interleukin 1beta secretion
- increased interleukin 17 secretion
- increased interleukin 6 secretion
-pyruvate kinase deficiency
Energy, Nutrition and Homeostasis
-phosphatidylserine synthase 2
-glycerol kinase 2
-NADH dehydrogenase (ubiquinone) Fe-S protein 4, 18 kDa (NADH-coenzyme Q reductase)
(Leigh syndrome )
-NAD(P)H dehydrogenase, quinone 1
(childhood acute lymphoblastic leukemia )
-solute carrier family 22 (organic cation/carnitine transporter), member 5
(primary carnitine deficiency )-syntaxin 4
-abnormal phospholipid level
-increased circulating follicle stimulating hormone
-decreased circulating glucose level
-increased fatty acid level
-abnormal body weight
-decreased body temperature
-decreased oxygen consumption
-abnormal liver morphology
-decreased circulating glucose
-abnormal glucose homeostasis
-decreased circulating carnitine
-increased circulating ammonia level
-impaired glucose tolerance
-increased circulating insulin level
-type II diabetes melitus Congts
(hyper-IgM syndrome )
-caspase 9, apoptosis-related
-bridging integrator 3
-Ras association (RalGDS/AF-6) domain family member 1
-B-cell derived lymphoma
-increased sensitivity to oxidative stress
-decreased cellular sensitivity to gamma irradiation
-increased incidence of ionizing radiation induced tumors
-increased tumor incidence
-cutaneous mast cell tumors
-visceral mast cell tumors
-feline injection site sarcomas
-squamous cell carcinoma
-mammary gland tumors
-bridging integrator 3
-mal, T-cell differentiation protein
(protection against invasive pneumococcal disease, bacteremia, malaria and tuberculosis )
-NADH dehydrogenase (ubiquinone) Fe-S protein 4, 18 kDa (NADH-coenzyme Q reductase)
(Leigh syndrome )
-caspase 9, apoptosis-related cysteine peptidase
-snail homolog 1(Drosophila)
(colon cancer )
-abnormal lens fiber morphology
-abnormal optic nerve nerve morphology
-abnormal eye electrophysiology
-optic nerve atrophy
-decreased startle response
-abnormal olfactory epithelium morphology
- abnormal Meckel's cartilage morphology
Behavior, Neurological and Nervous System
-N-ethylmaleimide-sensitive factor attachment protein, alpha
-carbohydrate sulfotransferase 10
-transmembrane protein 176B
-solute carrier family 35, member C1
(leukocyte adhesion deficiency II )
-glycerol kinase 2
-NAD(P)H dehydrogenase, quinone 2
(Alzheimer's disease )
-NADH dehydrogenase (ubiquinone) Fe-S protein 4, 18 kDa (NADH-coenzyme Q reductase)
(Leigh syndrome )
-caspase 9, apoptosis-related cysteine peptidase
-RAB3B, member RAS oncogene family
-snail homolog 1(Drosophila)
(colon cancer )
-mal, T-cell differentiation protein
(protection against invasive pneumococcal disease, bacteremia, malaria and tuberculosis )
-abnormal motor coordination
-impaired righting response
-abnormal learning and memory
-abnormal spatial learning, ataxia
-abnormal maternal nurturing
-abnormal suckling behavior
-abnormal motor learning
-abnormal spatial learning
-abnormal spatial working memory
-abnormal emotion and affect
-abnormal motor control
-abnormal nest building
-abnormal stationary movement
-abnormal CNS synaptic transmission
-reduced long term potentiation
-abnormal excitatory post synaptic potential
-abnormal brain commissure morphology
-abnormal brain development
-abnormal embryonic neuroepithelium layer differentiation
-decreased neuron apoptosis
-decreased neurotransmitter release
-enhanced paired-pulse facilitation
-open neural tube
-abnormal cerebellar granuale layer
-abnormal Purkinje cell layer, small cerebellum
-abnormal brain ventricle morphology
-abnormal cerbral cortex morphology
- abnormal forebrain morphology
-abnormal neuron morphology
-abnormal neuron physiology
-astrocytosis, brain vacuoles
-lysosomal storage diseases -idiopathic vestibular disease
-congenital unilateral vestibular disease
The cardiac module consists of eight genes and is associated with the following eight phenotypes: cardiac hypertrophy, dilated dorsal aorta, abnormal mitral valve morphology, abnormal cardiac output, abnormal myocardial fiber physiology, enlarged heart, abnormal outflow tract and abnormal coronary artery morphology. This module contains genes that are of relevance to feline cardiac disease such as hypertrophic cardiomyopathy and developmental defects of the heart.
The developmental-patterning module consists of seven genes and is associated with phenotypes that include abnormal mesoderm development, abnormal proximal/distal developmental patterning and abnormal rostral/caudal developmental patterning. Within this module we identified genes associated with distinct cell differentiation and specification properties such as embryonic growth arrest, abnormal trophoblast layer morphology and abnormal white adipose tissue. Additional phenotypes within this module were associated with retinal formation, renal function, intestine morphology as well as cholesterol, triglyceride and corticosterone levels. The phenotypes within this module may be useful in dissecting the genetic mechanisms underlying inherited developmental abnormalities in both domestic and endangered felids.
The third module is an immune and hematopoietic module that contains nine genes and represents phenotypes associated with specific cell types and lineages including macrophage physiology, spleen germinal cell number, granulocyte number, platelet number, T-cell and B-cell proliferation and hematopoiesis. Furthermore, this module exhibited phenotypes associated with susceptibility and resistance to pathogens such as abnormal immune system biology, abnormal class switch recombination and altered rate of infection. Some of the phenotypes in this module, including abnormal somatic hypermutation frequency and lymphoid hyperplasia were related to cancer, perhaps representing the immune surveillance component to the control of tumorigenesis within the body. Finally, some of the phenotypes within this module were associated with the modulation of specific immunologically important molecules such as cytokine secretion, interferon secretion, IgE levels, IgG1levels, and IgM levels. Genes within this module may offer some insight into feline specific immunological and inflammatory disorders.
The fourth module, energy/nutrition and homeostasis consists of six genes and exhibits a number of phenotypes associated with energy production and regulation within cells. Some of these phenotypes include decreased circulating glucose level, decreased oxygen consumption, abnormal gluconeogenesis and increased glucagon. Other phenotypes include endocrine level regulation of the organism such as abnormal body weight and decreased body temperature. Additionally, there were phenotypes associated with diseases of energy metabolism such as diabetes, these phenotypes included insulin resistance, abnormal glucose homeostasis and increased circulating insulin level. These phenotypes provide a context for better understanding of the unique nutritional and energy requirements of the cat.
The fifth module has five genes and encodes a tumorigenesis module associated with the following phenotypes: B-cell derived lymphoma, increased sensitivity to oxidative stress, increased apoptosis, decreased cellular sensitivity to gamma irradiation, increased incidence of ionizing radiation induced tumors, increased tumor incidence, malignant tumors, adenocarcinoma. The genes in this module may provide a useful gene set for investigating the genetic basis of feline lymphoma and carcinoma.
Module six is a sensory systems module, containing five genes, and is associated with the following visual phenotypes: abnormal lens fiber morphology, cataracts, abnormal optic nerve morphology, abnormal eye electrophysiology, abnormal vision, blindness and optic nerve atrophy. Cats exhibit vision related abnormalities under certain nutritional deficiencies; the genes associated with these phenotypes may provide a better understanding of the observed link between feline nutrition and visual function. Other phenotypes within this module include both hyperekplexia and decreased startle response, which may underlie feline adaptations required for successful predation.
The seventh module is a behavioral/neurological and nervous system set that contains 11 genes. The behavioral phenotypes arising from this module span traits as diverse as motor coordination and balance through learning, memory and gait. Additional phenotypes in this module are associated with emotion and affect as well as vocalization and maternal behavior. Within this module, we identified a number of phenotypes underlying neuronal specific physiological mechanisms such as altered synaptic transmission, altered long term potentiation, abnormal excitatory post synaptic potentials and decreased neurotransmitter release. This module contains a variety of developmentally important nervous system phenotypes having anatomical or histological annotations. These include abnormal brain commissure morphology, abnormal brain development, abnormal embryonic neuroepithelium layer differentiation as well as open neural tube, abnormal cerebellar granule layer, abnormal Purkinje cell layer, small cerebellum, abnormal brain ventricle morphology, abnormal cerebral cortex morphology and abnormal forebrain and hindbrain morphology. Finally, we identified specific CNS phenotypes of clinical importance such as abnormal neuron morphology, abnormal neuron physiology, astrocytosis, brain stem haemorrhage, gliosis and inter cranial haemorrhage.
We chose to focus on a relatively small number of gene-phenotype relationships in order to explore a relatively high resolution picture of important feline phenotypes that may be representative of our cDNA sequences. Our goal was to determine if any of our cDNA sequences were associated with phenotypes that may be of value in understanding the genetic basis of feline specific biology. Our analysis demonstrates that some of our cDNA sequences are indeed associated, through comparative genomics sequence analysis using the mammalian phenotype browser database, with phenotypes that are extremely important in feline health and disease. These modules and related genes provide an important and extremely useful candidate gene set for domestic cat functional genomics.
Orthologous OMIM Diseases
List of OMIM Diseases
2-methyl-3-hydroxybutyryl-CoA dehydrogenase deficiency
Glycogen storage disease, type 0
Gonadal dysgenesis, 46XY, partial, with minifascicular neuropathy
Griscelli syndrome, type 2
3-methylglutaconic aciduria, type I
Acyl-CoA dehydrogenase, short-chain, deficiency of
Hemolytic anemia due to bisphosphoglycerate mutase deficiency
Adrenal cortical carcinoma
Homocysteine plasma level
Aldolase A deficiency
Amyotrophic lateral sclerosis 10
Hyperleucinemia-isoleucinemia or hypervalinemia
Arthrogryposis multiplex congenita, distal, type 1
Hypervalinemia or hyperleucine-isoleucinemia
Bartter syndrome, type 2
Hypokalemic periodic paralysis
Hypomagnesemia, renal, with ocular involvement
Hypotrichosis, localized, autosomal recessive
Immunodeficiency with hyper IgM, type 4
Breast cancer, sporadic
Brugada syndrome 2
Leukoencephalopathy with vanishing white matter
Lipoid adrenal hyperplasia
Cardiomyopathy, dilated, 1M,
Megakaryoblastic leukemia, acute
Cardiomyopathy, dilated, 1N
Mental retardation, X-linked syndromic
Cardiomyopathy, dilated, 1Z
Methemoglobinemia due to cytochrome b5 deficiency
Cardiomyopathy, familial hypertrophic
Methylmalonyl-CoA epimerase deficiency
Carnitine acetyltransferase deficiency
Microphthalmia, syndromic 6
Carnitine deficiency, systemic primary
Mitochondrial complex I deficiency
Cataract, posterior polar 2
Myopathy due to phosphoglycerate mutase deficiency
Cerebral dysgenesis, neuropathy, ichthyosis, and palmoplantar keratoderma syndrome
Myopathy with exercise intolerance, Swedish type
Ceroid lipofuscinosis, neuronal 8
Charcot-Marie-Tooth disease, axonal, type 2F
Oral-facial-digital syndrome 1
Charcot-Marie-Tooth neuropathy, X-linked dominant, 1
Colon cancer, advanced
Combined oxidative phosphorylation deficiency 2
Phenylketonuria - dihydropteridine reductase deficiency
Combined oxidative phosphorylation deficiency 5
Phosphoglycerate kinase 1 deficiency
Congenital disorder of glycosylation, type IIc
Porphyria cutanea tarda
Cutis laxa, autosomal dominant
Ribose 5-phosphate isomerase deficiency
Deafness, autosomal recessive 63
Spastic paraplegia 31
Spondylocostal dysostosis, autosomal recessive 3
Diamond-Blackfan anemia 6
Epilepsy, neonatal myoclonic, with suppression-burst pattern
Temperature-sensitive apoptosis, cellular
Esophageal carcinoma, somatic
Transcobalamin II deficiency
Tyrosinemia, type I
Generalized epilepsy with febrile seizures
Ventricular tachycardia, catecholaminergic polymorphic, 2
Glutamine deficiency, congenital
Within the set of OMIM diseases, we identified biochemical and metabolic diseases such as disorders of oxidative phosphorylation and glycosylation as well as D-2-hydroxyglutaric aciduria, glycogen storage disease, phenylketonuria due to dihydropteridine reductase deficiency and phosphoglycerate kinase 1 deficiency. Among the disease annotations associated with cancers, we found that our cDNA sequences were associated with specific types of OMIM annotations including breast cancer, colon cancer, esophageal carcinoma, lung cancer, pancreatic cancer and ovarian cancer to name a few. We also identified diseases of the sensory systems, such as cataracts and deafness. Finally, we discovered a variety of orthologs of human genes implicated in specific disorders, including Leigh syndrome, Hyper-IgD syndrome, immunodeficiency associated with hyper IgM, Griscelli syndrome, STAR syndrome and Charcot-Marie-Tooth disease, retinitis pigmentosa and generalized epilepsy. Together, these OMIM annotations provide a diverse picture of the genes across diseases and offer a unique context for understanding the role of these cDNA sequences in feline health and disease.
We identified 1227 feline cDNA sequences derived from tissues obtained from ten cats and performed extensive comparative genomics functional analysis to elucidate the computationally derived comparative gene expression analysis patterns, biochemical functions and phenotypes associated with these sequences. Our cDNA sequences and associated comparative and functional analysis provide an initial perspective on feline biology as viewed through our set of 1227 cDNA sequences. Although it is predicted that the number of feline protein coding genes encoded in the cat genome is in the order of 20,000 to 25,000, similar to most other mammalian genomes, the number of known published cat protein coding gene sequences is much lower at 2099 sequences (NCBI databse, 2011). These 1227 cDNA/gene sequences represent a rich set of potential targets for genetic association studies, biologically relevant diets and pharmacologically active compounds which can be developed to enhance the well-being of companion cats worldwide. Additionally, these sequences have value in similar applications for endangered felids.
Our strategy to identify a set of 1227 high quality and high confidence cDNA sequences from feline tissue samples expands the expressed sequence data for domestic cat. Although we initially obtained over 3000 cDNA sequences, we chose to filter our sequences so that the set we describe would be of the most value for the feline genomics community. Specifically, the conservative strategy outlined in Figure 1 resulted in a set of 913 known sequences and 314 novel sequences (1227 sequences in total) of which 914 orthologous clusters across feline, human, dog and mouse were identified (for which 844 were known cDNA sequences and 70 were novel cDNA sequences). The genes corresponding to these 914 orthologs were used as input sequences for a variety of bioinformatics and computational analyses aimed at providing an initial perspective on the physiological and pathological roles of these sequences in feline development, nutrition and health. Although we have identified a number of interesting results using computational and sequence comparison methods, our analysis only identifies the potential roles of these genes based on comparative analysis in other species. However, validating these results and proving the function of these genes will require molecular and biochemical experimental analysis. The results of our inferred expression analysis provide a set of gene expression patterns consistent with the source tissues used for cDNA production. Of the 21 source tissues used as starting material, inferred expression patterns from each anatomical region were detected with greater than 100 genes being associated in each case. It is interesting to note that each of these tissues exhibited relatively high gene expression numbers (i.e., numbers of genes associated with anatomical expression), which is what one would expect if the inferred expression patterns were an accurate representation of the true expression patterns of the source tissues. Tissues such as brain (725 genes), heart (629 genes), pancreas (568 genes) and testis (703 genes) exhibit inferred expression of more than 60% of the genes encoding our 1227 cDNA sequences. Inferred cellular expression patterns correlated with cell types expected in the source tissues including glial cells and neurons (432 genes and 124 genes respectively), retinal pigment epithelium cells (514 genes), and skeletal muscle cells (499 genes). Together, these results provide an expression framework for understanding the roles of these cDNA sequences in feline physiology and pathology. Because greater than 70% of our cDNA sequences were associated with embryological expression patterns we were not surprised to discover that a significant number of developmental phenotypes were associated with our set of cDNA sequences. Specifically, we identified genes associated with abnormal heart morphology and abnormal cardiac blood flow, abnormal mesoderm development, abnormal developmental patterning and abnormal retinal neuronal layer morphology. These phenotypes are consistent with the expression and role of genes identified in the source tissues selected for cDNA sequencing. The fact that the inferred expression patterns exhibit greater breadth of expression than the starting tissues is in line with the notion that genes tend to be expressed in complex spatial and temporal patterns. It may be the case that the inferred expression patterns include some anatomical, cellular and/or developmental expression patterns which may be false positives, however the overall picture of expression provided by this analysis greatly enhances the value of these cDNA sequences in genomic applications.
Interestingly, our analysis of gene ontology in the context of dN/dS values of individual orthologous cDNA sequences provides insight into how the domestic cat is both similar to and differs from other mammals. We detected evidence of negative selection acting on genes associated with microtubules and the actin cytoskeleton, suggesting that genes associated with these cellular structures are fairly well conserved among mammals [41, 42]. Additionally, we identified gene ontology annotation terms affiliated with the nucleus, the chromosomes and DNA replication exhibiting relatively low values of dN/dS along with orthologs associated with transcriptional regulation and translational elongation. Similar values were obtained for genes annotated as G-protein beta/gamma binding and trans-Golgi network trafficking, vesicle and endoplasmic reticulum compartment membrane and SNARE complex. This is not surprising given that the housekeeping functions of mammalian cells are relatively well conserved. All cells must transmit information from the genome into RNA and protein components in a manner that maintains the appropriate subcellular compartmentalization of molecular functions. Intracellular trafficking that diverges from cellular requirements is likely to exhibit relatively deleterious consequences leading to negative selection compared to cells that function appropriately. Microtubules are involved in cellular integrity, cell motility and cell division; all of these processes are critical for cell viability [41, 43].
In comparison to these highly conserved orthologs which mediate the core cellular processes, we detect evidence of considerably less negative selection acting on orthologs associated with transmembrane receptors, apoptotic signals, guanyl-nucleotide exchange factors and GPCR activity. Additionally, we identified evidence of less negative selection among orthologs associated with extracellular spaces, mitochondrial membrane affiliation and integral proteins of the plasma membrane. Unlike the highly conserved orthologs with intracellular functions, these orthologs form the basis of interactions across cells, through the extracellular space into the nucleus and organelles by a variety of signal transduction mechanisms for which multiple paralogous genes exist in each species. Such patterns of selection have been identified by others and represent evolutionary patterns of selection that may be associated with positive selection in different evolutionary lineages . Moreover, these cDNA sequences might encode proteins for which extracellular environment plays a selective role during evolution.
It is well documented that paralogs diverge at a greater rate than orthologs [45, 46]. Because our analysis did not include the entire set of genes from the cat, we cannot rule out the possibility that some of our orthologs are not true orthologs. It is worthwhile to point out that our analysis included only cat, dog, mouse and human genes which effectively limits the detection of evolutionary selection using the dN/dS ratio because some of these species diverged more than 100 million years ago. Nonetheless, it is interesting that others have observed similar patterns of divergence in protein networks operating at the cellular periphery and within the extracellular space [44, 47].
Our analysis identified orthologs associated with respiratory chain and mitochondria as exhibiting relatively lower levels of negative selection. It is possible that the predatory status of cats resulted in adaptive changes in energy production and oxidative phosphorylation that facilitate the high energy requirements of predation.
It is interesting that we detect evidence of divergence within apoptotic genes in the cat compared to other mammalian species. This may underlie species specific differences in adaptation, such as what might be expected to have happened as obligate carnivores diverged from a common ancestor of omnivores and herbivores. The high protein requirements coupled with enhanced predatory fitness may have co-evolved with differences in cellular response to stress and cellular apoptosis, both within and outside of the brain.
This hypothesis is supported by the metabolic network analysis in GeneGO where the top 25% dN/dS values were associated with metabolic pathways implicated in non-carbohydrate roles. The metabolic network analysis performed with GeneGO demonstrated that genes in the group with smaller dN/dS values are associated with metabolic networks most involved in carbohydrate metabolism, while the genes in the larger dN/dS value group are in metabolic networks most involved in amino acid metabolism. This suggests that depending on metabolic requirements, the evolution rate may not be the same across all metabolic networks, and obligate carnivores like cats, may exhibit relatively less negative selection acting on genes involved in amino acid metabolism and more neutral selection acting on carbohydrate associated genes. This result is in agreement with the observation that cats exhibit different dietary requirements for amino acids taurine , arginine , cysteine and, methionine . In contrast to dogs, cats are unable to synthesize taurine from cysteine , subsequently, taurine deficiency in cats is associated with a variety of clinically important conditions including cardiac  immune , neurological , platelet , reproductive  and retinal  dysfunctions. Additionally, cats exhibit rapid onset of ammonia toxicity resulting from arginine deficiency and, in severe cases, may die within 24 hours [23, 48].
Through the use of KEGG pathway annotation, we identified domestic cat genes involved in a variety of amino acid related pathways including the metabolism of alanine, aspartate, arginine, proline, glutamate, glycine, serine, threonine, histidine, lysine, methionine, phenylalanine, tyrosine and tryptophan. We identified specific pathways in amino acid metabolism, which tend to differ between obligate carnivores and omnivorous mammals . These include six genes involved in tryptophan metabolism which are of value for cats because they are unable to synthesize niacin from tryptophan, as compared to omnivores . Additionally we identified three genes involved in arginine metabolism, which is an essential amino acid in cats . We identified genes involved in glutamate metabolism, which may provide insight into the metabolic consequences of the low levels of ornithine produced from glutamate in cats .
We also identified genes associated with pathways underlying lipid metabolism, including genes participating in biochemical pathways of linoleic, alpha-linoleic acid and arachidonic acids, which is important and noteworthy because cats cannot use linoleic acid for the biosynthesis of arachidonic acid . Further analysis of these genes may provide clues about feline biochemistry associated with arachidonic acid which may be important in feline reproduction . Finally, we identified genes involved in the metabolism of retinol, which represent another very important gene set because cats are unable to synthesize retinol from beta-carotene .
The metabolism and biosynthesis of cofactors, vitamins and glycans is important in the nutrition and health of animals. Within these biochemical pathways, we identified three genes associated with folate metabolism, seven genes involved in glutathione metabolism and two genes associated with keratin sulfate biosynthesis, two genes associated with N-glycan biosynthesis and three genes associated with pathothenate and CoA biosynthesis. Some of these genes may provide value as important biological markers for monitoring oxidative stress, apoptosis and immune function in cats .
Collectively, many of these genes and their associated pathways are important for feline health and nutrition because they represent biochemical processes that cats have adapted to accommodate the narrow dietary range of an obligate carnivore in contrast to omnivorous mammals. The subsequent characterization of these genes and pathways may provide a genomic foundation for understanding how obligate carnivores differ from other animals in both health and disease.
Our functional and evolutionary analysis suggests that through divergent evolutionary trajectories, different species evolve slightly different biochemical processes of cells, tissues and organs that contribute to the manifestation of species specific adaptations and disorders. The domestic cat is known to suffer from a number of hereditary diseases, many of which have counterparts in other species like humans and dogs . As part of our investigation into the biological significance of our cDNA sequences, we employed a comparative genomics approach to discover the phenotypes associated with these sequences. Our approach leveraged the mammalian phenotype ontology that has been developed as part of the mouse genome database . We decided to select a relatively small number of genes for which a considerable number of important phenotypes may be associated.
Our phenotype data was obtained from previously published mouse phenotyping studies using transgenic or knockout mice. Subsequently, they should be considered as related to, rather than exactly, the true phenotypes that might arise in the cat. Because our method relies upon orthologous relationships between cat and mouse genes, it is worthwhile to point out that inaccurate mappings between orthologs may lead to inaccurate predictions of phenotypes. Furthermore, as we have described throughout this paper, the cat exhibits some strong similarities to general biological processes that are shared with mammals. The cat also has well documented differences when compared to omnivorous animals. Therefore, one must consider the phenotype analysis as a general thematic picture of the functional consequences of our cDNA sequences rather than as a one-to-one mapping of gene-phenotype associations within our cDNA sequences.
We identified seven phenotypic modules exhibiting 136 phenotypes arising from only 38 genes. Many of the genes we identified exhibit numerous phenotypes, both within and across modules. Such pleiotropic effects underlie the complexity of mammalian genomes and provide context for future genomic studies. We selected these gene-phenotype associations to provide a detailed, but yet tractable picture of how our cDNA sequences might map onto anatomical and physiological traits.
Within the cardiac module, we identified eight genes associated with phenotypes relating to cardiac disease in cats. Some of the genes within this module include tropomodulin 1, snail homolog 1 and an interleukin receptor antagonist. This module includes phenotypes of cardiac hypertrophy and mitral valve defects, both of which are known hereditary diseases in cats . These genes provide examples of the types of phenotypes that might arise from perturbations of cat genes underlying inherited feline cardiac diseases, such as aortic stenosis, atrial-septal defect, mitral valve displasia, tetralogy of Fallot and ventricular-septal defect [53, 54].
Our developmental module consists of seven genes and includes a TGFbeta induced homeobox transcription factor as well as the signaling molecule arginine-vasopressin. The phenotypes associated with this module include developmental patterning across both the proximal/distal axis and the rostral/caudal axis. The phenotypes also include cellular specification and patterning such as mesoderm development, trophoblast layer morphology and adipose tissue differentiation, to name a few. Domestic cats exhibit a variety of developmental defects, such as polydactyly, hip dysplasia, sacrococcygeal dysgenesis, portocaval shunt, open central fontanel, open lateral fontanel and thoracic hemivertebra [54–57]. The cDNA sequences we describe may include genes that are responsible for abnormal developmental conditions in domestic and endangered felids.
We identified a sensory module, which contains five genes such as NADH dehydrogenase (ubiuinone) Fe-S protein 4 and caspase 9 apoptosis-related cysteine pepidase. This module includes the phenotypes of cataracts, blindness and optic nerve atrophy. Examples of inherited sensory system disorders in the domestic cat include cataracts, corneal dystrophy (stromal and endothelial), progressive retinal atrophy and glaucoma . The overlap between retinal and ocular phenotypes and inherited feline diseases suggests that there are specific genomic regions, represented by our cDNA sequences, which may include aspects of the genetic mechanisms of these debilitating diseases in cats. It is interesting to note that our sensory module includes genes involved in energy production. This is not surprising as retinal tissue is known to exhibit relatively high energy requirements and depletion of energy in this tissue has been associated with blindness and other vision defects .
Within our energy and homeostasis module, we identified genes like glycerol kinase 2, NAD(P)H dehydrogenase quinone 1 and NADH dehydrogenase (ubiquinone) Fe-S protein 4. The phenotypes within this module are associated with traits of clinical and adaptive importance in the cat. For example, our comparative phenotype analysis identified phenotypes of insulin resistance, increased circulating insulin level and impaired glucose tolerance; traits associated with the feline hereditary disease of diabetes mellitus . This module also contains phenotypes such as abnormal gluconeogenesis, increased glucagon, abnormal glucose homeostasis and increased circulating ammonia level, which are important in felid nutrition as cats use gluconeogenesis as a predominant form of energy production and are susceptible to ammonia toxicity [17, 18]. The genes in this module are of value in exploring some of the fundamental metabolic and biochemical differences between obligate carnivores and omnivores. Moreover, these genes may provide a genomic basis for specific diets that can reduce the incidence of feline disorders associated with specific nutritional deficiencies.
Within other modules, we identified phenotypes associated with cancer, such as increased tumor incidence, malignant tumors and B-cell derived lymphoma which may provide clues to the genetic susceptibility cats have for hereditary lymphoma . Among the behavioral phenotypes within the nervous system module, we identified a number of traits that may represent predator specific adaptations of cats. For example, we identified cDNA sequences associated with spatial learning, balance, righting response, gate and motor coordination; traits that are almost synonymous with cats and of extreme adaptive value for an apex hyper predator.
The comparative genomics analysis of OMIM diseases within our cDNA sequence data set provides a final perspective on the importance of our reported sequences in the health of domestic cats. Many of the diseases identified in the OMIM mapping are also represented by phenotypes within the modules. This independent annotation demonstrates that our analysis converges even though OMIM analysis leverages human orthology relationships and the phenotype analysis leverages murine orthology relationships. It is worth noting the limitation of sequence based comparative genomics approaches. They can provide considerable insight into the functional role of our cDNA sequences, but must ultimately be proven through focused and carefully designed genomics studies in cats. Nonetheless, our cDNA sequences and associated analysis provide considerable value through the identification of many interesting clinically and nutritionally relevant feline genes.
The set of diseases and phenotypes provides a starting point for candidate gene approaches and for the selection of biomarkers for monitoring nutrition and health. By combining diverse types of annotation, we can better understand the function of a given gene in a breadth of tissues and organ systems and of the biological processes it is involved in the organismal level, as well as its role in disease. For example, we identified genes associated with expression in the heart, and with a number of cardiac phenotypes, including cardiac hypertrophy, abnormal outflow tract and abnormal mitral valve morphology, as well as the OMIM disease annotation of dilated cardiomyopathy. These are of direct relevance to feline disease, since hypertrophic cardiomyopathy is a common clinical concern in cats .
The recent development of a 70,000 SNP feline bead array by Hill's Pet Nutrition and the Morris Animal Foundation provides an important and powerful resource for conducting gene association studies in the domestic cat, and related endangered species. However, even in the absence of whole-genome genetic association approaches, our characterization of these 1227 cDNA sequences provides an extremely valuable resource for candidate gene approaches aimed at investigating the genetic basis of feline phenotypes. It will be interesting to see how our comparative and functional analysis of these 1227 cDNA sequences compares to the data produced from high throughput sequencing and future genetic studies within and across different breeds in the domestic cat. It is likely that some of our functional annotations may turn out not to hold, and it is equally likely that some of them will. Through collaborative efforts, it will be possible to begin unravelling the genetic mechanisms underlying feline health and disease.
We report the identification of 1227 feline cDNA sequences of which, 913 correspond to higher quality versions of public feline sequences and 314 correspond to novel feline sequences for which no known public sequence data exists. Our comprehensive functional analysis identified a number of physiologically important biochemical pathways that these sequences are involved in as well as of the developmental, clinical and nutritional relevant phenotypes they are associated with.
Construction of feline tissue specific cDNA libraries
The study protocol was reviewed and approved by the Institutional Animal Care and Use Committee. All cats were immunized against feline panleukopenia, calici, rhinotracheitis, and rabies. Cats were housed with 10 - 12 other cats and food was continuously available throughout the day until their daily caloric requirements were consumed. Cats were housed in spacious rooms with natural light that varies with seasonal changes. Cats experienced behavioral enrichment through interactions with each other, by daily interaction and play time with caretakers, large windows and sun porches to watch the natural landscape and access to toys. At the end of their natural life, cats were euthanized for humane purposes and tissues were stored at -80C.
Total RNA was purified from 21 feline tissues (brain, kidney medulla/cortex, spleen, heart, liver, lung, skeletal muscle, thyroid gland, lymph node, pancreas, adrenal gland, tongue, colon, mammary gland, neonatal thymus, brain and testes) collected from 10 domestic short-haired cats postmortem, three cell lines derived from kidney, brain, lung, and 1 tissue pool using standard procedures as described in . The purity and integrity of each RNA sample was assessed by spectrophotometry and gel electrophoresis. Forty normalized cDNA libraries were constructed by Agencourt Inc. (Beckman-Coulter Genomics), 22 with standard inserts (1.2 kb) and 18 with long inserts (> 4 kb). The first and second cDNA strands were synthesized using optimized methods, and cDNAs were size selected prior to cloning. The size-selected cDNAs were directionally cloned into the pAGEN vector by polishing and restriction digest, creating a 5' blunt end and a 3' overhang.
Each cDNA library was subsequently tested for specific quality control measures (average insert size, number of independent clones and percentage of recombinant clones), and normalized to reduce the proportion of highly abundant mRNAs. Normalization was performed by dividing each library into two populations, using the first for in vitro transcription of biotinylated RNA, and the second to generate single stranded phagemid DNA. The two populations were then mixed, and self-hybridized DNA-RNA molecules corresponding to over-represented mRNAs were removed. The remaining single stranded DNA molecules were primed for second strand synthesis and the resulting clones were transformed into bacteria, yielding the normalized libraries.
Sequencing of feline cDNA libraries
Plasmids were purified from each library using a large-scale automated protocol, the SprintPrep® Solid Phase Reversible Immobilization procedure. Sequencing reactions were performed in 384-well plates using BigDye® Version 3.1 direct cycle sequencing (Applied Biosystems, CA). Sequencing reactions were purified using the CleanSeq® dye-terminator removal kit (Agencourt, Inc.), and resolved by capillary electrophoresis using the ABI3730 Genetic Analyzer (Applied Biosystems, CA). Sequencing reads were processed using Phred and quality scores for each run were monitored using the Agencourt, Inc. Galaxy LIMS system. Sequencing of these cDNA libraries yielded a total of 919,676 EST reads.
Data Management and Analysis
The sequence data, annotation data and the data resulting from sequence analysis were loaded into the MySQL relational database version 5 to facilitate data management and analysis .
Sequence Filtering and Ortholog Detection
A set of 3035 full length feline cDNA sequences were obtained from the analysis of the sequencing data and used to identify a set of high confidence cDNA sequences. All cDNA sequences were translated in 6 reading frames and the longest protein coding sequence obtained was noted. These cDNA and protein sequences were clustered using blast to identify a set of non-redundant nucleotide and non-redundant protein sequences using a stringency of 95% or greater as criteria for identifying redundant sequences. For each cluster, the longest representative sequence was chosen as the non-redundant representative. The intersection of non-redundant nucleotide sequences and non-redundant protein sequences was used as the set of non-redundant sequences.
The BLAST programs, blastp and blastn [62, 63], were run with the non-redundant full length feline sequences as query and the target species sequences downloaded from ENSEMBL ftp://ftp.ensembl.org/pub/current/fasta/  as subject sequences. The subject sequences for each of the four species (dog, cat, human and mouse) were: Homo_sapiens.GRCh37.60.cdna.all.fa (containing 147,141 sequences), Homo_sapiens.GRCh37.60.pep.all.fa (containing 81,968 sequences), Mus_musculus.NCBIM37.60.cdna.all.fa (containing 82,508 sequences), Mus_musculus.NCBIM37.60.pep.all.fa (containing 50,959 sequences),
Canis_familiaris.BROADD2.60.cdna.all.fa (containing 27,301 sequences),
Canis_familiaris.BROADD2.60.pep.all.fa (containing 25,559 sequences),
Felis_catus.CAT.60.cdna.all.fa (containing16,332 sequences) and
Felis_catus.CAT.60.pep.all.fa (containing 15,048 sequences). Because the human sequence sets contain the greatest number of target sequences, 147,141 nucleotide sequences and 81,968 protein sequences, the set of non-redundant sequences were mapped to the human sequences. Additionally, the full length sequences were mapped to the set of known feline cDNA and protein sequences in order to classify the full length non-redundant feline sequences as either known or novel, where known indicates that the sequence is represented by a feline sequence in the public ensembl transcript/protein sequence data while novel indicates that the sequence does not have a representative transcript or protein sequence in the ensembl data set.
Because the public feline data does not contain all of the protein coding genes, it was not possible to perform an ortholog search using the standard reciprocal best hit approach. Instead, the blast results were filtered using an iterative heuristic process of selecting blast hits with specific match lengths, gaps, number mismatches and percent identity. In total, eight iterative steps were performed beginning with the most stringent and ending with the least stringent. Each step identified a set of qualifying non-redundant full length sequences. The first and most stringent step imposed the requirement that the blast match_length must be equal to the smallest of the two sequences (query or subject) and the number of mismatches = 0, number of gaps = 0, and the percent identity ≥ 99%. A second filter was used to add additional sequences to the results of the first step, and any sequences that had not been identified in the first step were added to the set of results. The second step used a blast match_length ratio of ≥ 0.99, number mismatches = 0, number gaps = 0 and percent identity ≥ 99%. A third step identified additional sequences that satisfied the third step criteria and for which the first two steps did select the non-redundant full length sequence. The third step criteria were blast match length ratio ≥ 0.87, number of mismatches ≤ 4, number of gaps = 0, and percent identity ≥ 99%. The iterative process continued for a total of eight steps with each subsequent step relaxing the filtering criteria in order to identify sequences that were not identified in the previous step. Fourth step criteria were blast match length ratio ≥ 0.725, number mismatches ≤ 5, number of gaps = 0, and percent identity ≥ 99%. Fifth step criteria were blast match length ≥ 0.69, number mismatches ≤ 4, number gaps ≤ 1 and percent identity ≥ 99%. Sixth step criteria included blast match length ≥ 0.625, number mismatches ≤ 8, number of gaps ≤ 1, and percent identity ≥ 98%. Seventh step criteria included blast match length ratio ≥ 0.575, number mismatches ≤ 13, number gaps ≤ 2 and percent identity ≥ 97%. The eighth step criteria were blast match length ≥ 0.52, number mismatches ≤ 12, number of gaps ≤ 2 and percent identity ≥ 97%.
The resulting set of non-redundant full length sequences were considered to represent the high quality feline cDNA and protein sequences. These high quality sequences which mapped to a known public feline sequence were used to generate global nucleotide and protein alignments using the partial order alignment software POA http://bioinfo.mbi.ucla.edu/poa2/POA_Online/Align.html [64, 65]. All alignments were manually inspected to ensure that each non-redundant full length feline sequence mapped to the correct public feline sequence.
Comparative Expression Analysis
In order to infer anatomical and cellular expression patterns of our sequences, four expression annotation files were downloaded from the public biomart http://www.biomart.org  web server. Because we mapped our sequences to their corresponding human orthologs, we downloaded the human biomart egenetics annotation data sets mapped on top of the ensembl gene 60 version human gene identifiers. The four annotation sets we obtained included human ensembl gene identifiers mapped to (1) a set of anatomical terms, (2) a set of cell types, (3) a list of pathological terms and (4) a list of developmental stages ranging from weeks to years.
Although our sequences represent a subset of gene products, we found value in identifying the spectrum of expression patterns these sequences may exhibit beyond the tissue libraries that we used. The mapping was accomplished by loading each of the four gene expression annotation files into the MySQL relational database and performing SQL queries that joined these expression tables to our orthologous gene set using the ensembl human gene identifier.
dN/dS Codon Substitution Rate Calculations
In order to better understand the evolutionary relationships between the feline cDNA sequences and the orthologous sequences in dog, human and mouse, we calculated dN/dS values for orthologous sequences across the different species. Phylogenetic Analysis by Maximum Likelihood (PAML version 4.4) software was used to run the codon stats using the "codeml" program. Codon stats were computed where it was possible (DNA and protein sequence availability for the orthologs) with basic model (NSSites = 0) ω = dN/dS, the ratio of nonsynonymous/synonymous substitution rates
The ω ratio is a measure of natural selection acting on the protein. Simplistically, values of ω < 1, = 1, and > 1 means negative purifying selection, neutral evolution, and positive selection respectively. PAL 2NAL [67–69] was used to create codon alignments between the cDNAs and the proteins to input to PAML program which computes the dN, dS and ω ratio. Codon substitution rate data was loaded into the MySQL relational database and used to assess the evolutionary pressure exerted on specific groups of genes. The gene groups were derived from other annotation types, such as gene ontology and phenotype annotation results.
Gene Ontology Annotation
Gene ontology annotation was added to the set of orthologous sequences via comparative sequence analysis. Because the orthologous sequences were already mapped to the human transcripts and proteins, we decided to download the gene ontology annotation files corresponding to biological process, molecular function and cellular localization in order to annotate the non-redundant full length feline sequences with gene ontology terms . The gene ontology annotation files linking the gene ontology terms to the ensembl human gene identifiers were obtained from biomart http://www.biomart.org. Each feline sequence we identified was annotated with all of the gene ontology terms associated with the orthologous human gene. In this manner, we were able to identify a larger set of gene ontology annotations per feline gene than we could have accomplished if we limited the annotation mapping to only the feline cDNA sequences we identified. Through this greedy algorithm, we were able to gain a more comprehensive understanding of the genes we identified. SQL queries in MySQL database were used to map the human gene ontology annotation terms to the orthologous feline genes encoding the cDNA sequences we identified.
GeneGO Metabolic Network Analysis
Metabolic networks of feline sequences (using ENSEMBL id to upload) were performed using the MetaCore software (GeneGO, St. Joseph, MI). MetaCore identifies networks based on a manually curated database containing known molecular interactions, functions, and disease interrelationships. The networks are identified by the probability that a random set of genes the same size as the input list would give rise to a particular mapping by chance. Therefore, an enrichment of biological relevant pathways or networks can be found.
KEGG Pathway Annotation
Pathway associations were identified using a comparative genomics approach. Because our orthologous sequences were mapped to human orthologs, it was possible to use the human pathway association information to map the pathways on the orthologous feline sequences. This was accomplished using SQL queries to join the KEGG http://www.genome.jp/kegg/pathway.html  and Biocarta http://www.biocarta.com/genes/index.asp pathway data  that has been associated with human ensembl gene identifiers with feline gene identifiers. Additional pathways were identified using the David Bioinformatics Database http://david.abcc.ncifcrf.gov/  through a gene set search using the ensembl gene identifiers for the set of human orthologs of the feline sequences we identified.
Comparative Phenotype Mapping
Gene specific phenotype annotation derived from mouse knockouts and/or transgenic strains is compiled and made publicly available at the Mouse Genome Database http://www.informatics.jax.org/ . The phenotype annotation is structured within the mammalian phenotype ontology which provides an acyclic graph of mammalian morphological and physiological phenotypes. Because the mouse phenotype data is associated with each mouse gene, it was possible to link the mammalian phenotype ontology to the feline non-redundant full length sequences through a two step process. First, the mammalian phenotype annotations linked to mouse gene identifiers were obtained and loaded into the MySQL database. Next, the appropriate SQL query was performed which created a table that joined the phenotype information with our feline sequence data. The resulting phenotype annotations on top of the feline orthologous gene set provide an additional mechanism for understanding the role of these cDNA sequences in cat development, health and disease.
OMIM Disease Mapping
A comparative genomics map of our feline sequences annotated with the OMIM http://www.ncbi.nlm.nih.gov/omim  disease information was generated using two different approaches. The first approach utilized MIM disease data that was produced from biomart and anchored to the human ensembl gene identifiers. The resulting annotation file was loaded into the relational database and an appropriate SQL query was used to connect the disease information to the feline sequences through the orthologous relationships that were previously determined. The resulting mapping provided formal associations between feline cDNA sequences and OMIM disease information http://www.ncbi.nlm.nih.gov/omim.
A second method of mapping the feline sequence data to the OMIM data was used to improve the set of OMIM annotated feline cDNA sequences. Specifically, the set of human ensembl gene identifiers corresponding to the orthologs for the feline cDNA sequences were used to query the David Bioinformatics database for OMIM disease information. The resulting file downloaded from the David Database contained human ensembl gene identifiers and OMIM disease identifiers. This file was loaded in the MySQL database and linked with the non-redundant feline cDNA sequences using an appropriate SQL query.
This work was supported by Hills Pet Nutrition, Inc. The authors would like to thank Dr. Steve Henriksen at Western University of Health Sciences for his help in providing the research infrastructure needed to complete this project. The authors would like to thank Dr. Teresa Morishita, Dr. Steve Waldhalm and Dr. Phil Nelson for freeing KJI from teaching responsibilities and providing additional research time needed to complete this project. The authors wish to thank Dr. Barry Watson for facilitating the creation of the Western University of Health Sciences comparative genomics partnership with Hills Pet Nutrition, Inc.
- Carnivora: Dogs, cats, bears, raccoons, weasels, mongooses, hyenas, seals, walruses, etc. [http://tolweb.org/Carnivora/15971/2000.01.01]
- Tarditi CR, Grahn RA, Evans JJ, Kurushima JD, Lyons LA: Mitochondrial DNA sequencing of cat hair: an informative forensic tool. J Forensic Sci. 2011, 56 (Suppl 1): S36-46.PubMed CentralView ArticlePubMedGoogle Scholar
- Bradshaw JW, Goodwin D, Legrand-Defretin V, Nott HM: Food selection by the domestic cat, an obligate carnivore. Comp Biochem Physiol A Physiol. 1996, 114 (3): 205-209.View ArticlePubMedGoogle Scholar
- Batchelor DJ, Al-Rammahi M, Moran AW, Brand JG, Li X, Haskins M, German AJ, Shirazi-Beechey SP: Sodium/glucose cotransporter-1, sweet receptor, and disaccharidase expression in the intestine of the domestic dog and cat: two species of different dietary habit. Am J Physiol Regul Integr Comp Physiol. 2011, 300 (1): R67-75.PubMed CentralView ArticlePubMedGoogle Scholar
- O'Brien SJ, Johnson WE: The evolution of cats. Genomic paw prints in the DNA of the world's wild cats have clarified the cat family tree and uncovered several remarkable migrations in their past. Scientific American. 2007, 297 (1): 68-75.View ArticlePubMedGoogle Scholar
- Johnson WE, Eizirik E, Pecon-Slattery J, Murphy WJ, Antunes A, Teeling E, O'Brien SJ: The late Miocene radiation of modern Felidae: a genetic assessment. Science. 2006, 311 (5757): 73-77.View ArticlePubMedGoogle Scholar
- Mota P, Amaral S, Martins L, de Lourdes Pereira M, Oliveira PJ, Ramalho-Santos J: Mitochondrial bioenergetics of testicular cells from the domestic cat (Felis catus)-a model for endangered species. Reprod Toxicol. 2009, 27 (2): 111-116.View ArticlePubMedGoogle Scholar
- Pettersson LG, Perfiliev S, Zotova E, Lundberg A: Role of claws and pads in taking and holding food in cats. Neurosci Res. 1998, 31 (4): 343-346.View ArticlePubMedGoogle Scholar
- Perfiliev S, Pettersson LG, Lundberg A: Control of claw movements in cats. Neurosci Res. 1998, 31 (4): 337-342.View ArticlePubMedGoogle Scholar
- Populin LC, Yin TC: Pinna movements of the cat during sound localization. J Neurosci. 1998, 18 (11): 4233-4243.PubMedGoogle Scholar
- Stecker GC, Mickey BJ, Macpherson EA, Middlebrooks JC: Spatial sensitivity in field PAF of cat auditory cortex. J Neurophysiol. 2003, 89 (6): 2889-2903.View ArticlePubMedGoogle Scholar
- Las L, Shapira AH, Nelken I: Functional gradients of auditory sensitivity along the anterior ectosylvian sulcus of the cat. J Neurosci. 2008, 28 (14): 3657-3667.View ArticlePubMedGoogle Scholar
- Jeffery N, Cox PG: Do agility and skull architecture influence the geometry of the mammalian vestibulo-ocular reflex?. J Anat. 2010, 216 (4): 496-509.PubMed CentralView ArticlePubMedGoogle Scholar
- Homman-Ludiye J, Manger PR, Bourne JA: Immunohistochemical parcellation of the ferret (Mustela putorius) visual cortex reveals substantial homology with the cat (Felis catus). J Comp Neurol. 2010, 518 (21): 4439-4462.View ArticlePubMedGoogle Scholar
- Buddington RK, Chen JW, Diamond JM: Dietary regulation of intestinal brush-border sugar and amino acid transport in carnivores. Am J Physiol. 1991, 261 (4 Pt 2): R793-801.PubMedGoogle Scholar
- Arai T, Washizu T, Sako T, Sasaki M, Motoyoshi S: D-glucose transport activities in erythrocytes and hepatocytes of dogs, cats and cattle. Comp Biochem Physiol Comp Physiol. 1992, 102 (2): 285-287.View ArticlePubMedGoogle Scholar
- Washizu T, Tanaka A, Sako T, Washizu M, Arai T: Comparison of the activities of enzymes related to glycolysis and gluconeogenesis in the liver of dogs and cats. Res Vet Sci. 1999, 67 (2): 205-206.View ArticlePubMedGoogle Scholar
- Silva SV, Mercer JR: Effect of protein intake on amino acid catabolism and gluconeogenesis by isolated hepatocytes from the cat (Felis domestica). Comp Biochem Physiol B. 1985, 80 (3): 603-607.PubMedGoogle Scholar
- Kettelhut IC, Foss MC, Migliorini RH: Glucose homeostasis in a carnivorous animal (cat) and in rats fed a high-protein diet. Am J Physiol. 1980, 239 (5): R437-444.PubMedGoogle Scholar
- Morris JH, Rogers QR: The metabolic basis for the taurine requirement of cats. Adv Exp Med Biol. 1992, 315: 33-44.View ArticlePubMedGoogle Scholar
- Morris JG: Nutritional and metabolic responses to arginine deficiency in carnivores. J Nutr. 1985, 115 (4): 524-531.PubMedGoogle Scholar
- Smalley KA, Rogers QR, Morris JG: Methionine requirement of kittens given amino acid diets containing adequate cystine. Br J Nutr. 1983, 49 (3): 411-417.View ArticlePubMedGoogle Scholar
- Stewart PM, Batshaw M, Valle D, Walser M: Effects of arginine-free meals on ureagenesis in cats. Am J Physiol. 1981, 241 (4): E310-315.PubMedGoogle Scholar
- Sturman JA, Gargano AD, Messing JM, Imaki H: Feline maternal taurine deficiency: effect on mother and offspring. J Nutr. 1986, 116 (4): 655-667.PubMedGoogle Scholar
- Chesney RW, Helms RA, Christensen M, Budreau AM, Han X, Sturman JA: The role of taurine in infant nutrition. Adv Exp Med Biol. 1998, 442: 463-476.View ArticlePubMedGoogle Scholar
- Morris JG: Idiosyncratic nutrient requirements of cats appear to be diet-induced evolutionary adaptations. Nutr Res Rev. 2002, 15 (1): 153-168.View ArticlePubMedGoogle Scholar
- Pion PD, Kittleson MD, Thomas WP, Skiles ML, Rogers QR: Clinical findings in cats with dilated cardiomyopathy and relationship of findings to taurine deficiency. J Am Vet Med Assoc. 1992, 201 (2): 267-274.PubMedGoogle Scholar
- Schuller-Levis G, Mehta PD, Rudelli R, Sturman J: Immunologic consequences of taurine deficiency in cats. J Leukoc Biol. 1990, 47 (4): 321-331.PubMedGoogle Scholar
- Saransaari P, Oja SS: Taurine in the developing cat: uptake and release in different brain areas. Neurochem Res. 1994, 19 (1): 77-82.View ArticlePubMedGoogle Scholar
- Welles EG, Boudreaux MK, Tyler JW: Platelet, antithrombin, and fibrinolytic activities in taurine-deficient and taurine-replete cats. Am J Vet Res. 1993, 54 (8): 1235-1243.PubMedGoogle Scholar
- Buff S, Donze A, Guerin P, Guillaud J, Fontbonne A, Menezo Y: Taurine and hypotaurine in spermatozoa and epididymal fluid of cats. J Reprod Fertil Suppl. 2001, 57: 93-95.PubMedGoogle Scholar
- Imaki H, Messing J, Sturman JA: Extensive taurine depletion and retinal degeneration in cats treated with beta-alanine for 40 weeks. Adv Exp Med Biol. 1998, 442: 449-460.View ArticlePubMedGoogle Scholar
- Warskulat U, Heller-Stilb B, Oermann E, Zilles K, Haas H, Lang F, Haussinger D: Phenotype of the taurine transporter knockout mouse. Methods Enzymol. 2007, 428: 439-458.View ArticlePubMedGoogle Scholar
- Pontius JU, Mullikin JC, Smith DR, Agencourt Sequencing Team, Lindblad-Toh K, Gnerre S, Clamp M, Chang J, Stephens R, Neelam B, Volfovsky N, Schäffer AA, Agarwala R, Narfström K, Murphy WJ, Giger U, Roca AL, Antunes A, Menotti-Raymond M, Yuhki N, Pecon-Slattery J, Johnson WE, Bourque G, Tesler G, NISC Comparative Sequencing Program, O'Brien SJ: Initial sequence and comparative analysis of the cat genome. Genome Res. 2007, 17 (11): 1675-1689.PubMed CentralView ArticlePubMedGoogle Scholar
- Mullikin JC, Hansen NF, Shen L, Ebling H, Donahue WF, Tao W, Saranga DJ, Brand A, Rubenfield MJ, Young AC, Cruz P, NISC Comparative Sequencing Program, Driscoll C, David V, Al-Murrani SW, Locniskar MF, Abrahamsen MS, O'Brien SJ, Smith DR, Brockman JA: Light whole genome sequence for SNP discovery across domestic cat breeds. BMC Genomics. 2010, 11: 406-PubMed CentralView ArticlePubMedGoogle Scholar
- Flicek P, Amode MR, Barrell D, Beal K, Brent S, Chen Y, Clapham P, Coates G, Fairley S, Fitzgerald S, Gordon L, Hendrix M, Hourlier T, Johnson N, Kähäri A, Keefe D, Keenan S, Kinsella R, Kokocinski F, Kulesha E, Larsson P, Longden I, McLaren W, Overduin B, Pritchard B, Riat HS, Rios D, Ritchie GR, Ruffier M, Schuster M: Ensembl 2011. Nucleic Acids Res. 2011, D800-806. 39 DatabaseGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5 (7): 621-628.View ArticlePubMedGoogle Scholar
- Ramskold D, Wang ET, Burge CB, Sandberg R: An abundance of ubiquitously expressed genes revealed by tissue transcriptome sequence data. PLoS Comput Biol. 2009, 5 (12): e1000598-PubMed CentralView ArticlePubMedGoogle Scholar
- Nielsen R: Statistical tests of selective neutrality in the age of genomics. Heredity. 2001, 86 (Pt 6): 641-647.View ArticlePubMedGoogle Scholar
- Blake JA, Bult CJ, Eppig JT, Kadin JA, Richardson JE: The Mouse Genome Database genotypes::phenotypes. Nucleic Acids Res. 2009, D712-719. 37 DatabaseGoogle Scholar
- Yamaguti GG, Lourenco GJ, Silveira VS, Tone LG, Lopes LF, Lima CS: Increased risk for acute lymphoblastic leukemia in children with cytochrome P450A1 (CYP1A1)- and NAD(P)H:quinone oxidoreductase 1 (NQO1)-inherited gene variants. Acta Haematol. 2010, 124 (3): 182-184.View ArticlePubMedGoogle Scholar
- Strater J, Herter I, Merkel G, Hinz U, Weitz J, Moller P: Expression and prognostic significance of APAF-1, caspase-8 and caspase-9 in stage II/III colon carcinoma: caspase-8 and caspase-9 is associated with poor prognosis. Int J Cancer. 2010, 127 (4): 873-880.PubMedGoogle Scholar
- Guo C, Ahmad T, Beckly J, Cummings JR, Hancock L, Geremia A, Cooney R, Pathan S, Jewell DP: Association of caspase-9 and RUNX3 with inflammatory bowel disease. Tissue Antigens. 2011, 77 (1): 23-29.View ArticlePubMedGoogle Scholar
- Aubier M, Levy J, Clerici C, Neukirch F, Cabrieres F, Herman D: Protective effect of theophylline on bronchial hyperresponsiveness in patients with allergic rhinitis. The American review of respiratory disease. 1991, 143 (2): 346-350.View ArticlePubMedGoogle Scholar
- Khor CC, Chapman SJ, Vannberg FO, Dunne A, Murphy C, Ling EY, Frodsham AJ, Walley AJ, Kyrieleis O, Khan A, Aucan C, Segal S, Moore CE, Knox K, Campbell SJ, Lienhardt C, Scott A, Aaby P, Sow OY, Grignani RT, Sillah J, Sirugo G, Peshu N, Williams TN, Maitland K, Davies RJ, Kwiatkowski DP, Day NP, Yala D, Crook DW, et al: A Mal functional variant is associated with protection against invasive pneumococcal disease, bacteremia, malaria and tuberculosis. Nat Genet. 2007, 39 (4): 523-528.PubMed CentralView ArticlePubMedGoogle Scholar
- Ginsberg JS, Caco CC, Brill-Edwards PA, Panju AA, Bona R, Demers CM, Tuters LM, Nugent P, McGinnis J, Grant BM: Venous thrombosis in patients who have undergone major hip or knee surgery: detection with compression US and impedance plethysmography. Radiology. 1991, 181 (3): 651-654.View ArticlePubMedGoogle Scholar
- Franci C, Gallen M, Alameda F, Baro T, Iglesias M, Virtanen I, Garcia de Herreros A: Snail1 protein in the stroma as a new putative prognosis marker for colon tumours. PLoS One. 2009, 4 (5): e5595-PubMed CentralView ArticlePubMedGoogle Scholar
- Peirson SN, Butler JN: RNA extraction from mammalian tissues. Methods Mol Biol. 2007, 362: 315-327.View ArticlePubMedGoogle Scholar
- Leshinsky-Silver E, Lebre AS, Minai L, Saada A, Steffann J, Cohen S, Rotig A, Munnich A, Lev D, Lerman-Sagie T: NDUFS4 mutations cause Leigh syndrome with predominant brainstem involvement. Mol Genet Metab. 2009, 97 (3): 185-189.View ArticlePubMedGoogle Scholar
- Gatalica Z, Lele SM, Rampy BA, Norris BA: The expression of Fhit protein is related inversely to disease progression in patients with breast carcinoma. Cancer. 2000, 88 (6): 1378-1383.View ArticlePubMedGoogle Scholar
- Jaeger PA, Wyss-Coray T: Beclin 1 complex in autophagy and Alzheimer disease. Archives of neurology. 2010, 67 (10): 1181-1184.View ArticlePubMedGoogle Scholar
- Quintana A, Kruse SE, Kapur RP, Sanz E, Palmiter RD: Complex I deficiency due to loss of Ndufs4 in the brain results in progressive encephalopathy resembling Leigh syndrome. Proc Natl Acad Sci USA. 2010, 107 (24): 10996-11001.PubMed CentralView ArticlePubMedGoogle Scholar
- Abbott JA: Feline hypertrophic cardiomyopathy: an update. Vet Clin North Am Small Anim Pract. 2010, 40 (4): 685-700.View ArticlePubMedGoogle Scholar
- Reddy S, Jia S, Geoffrey R, Lorier R, Suchi M, Broeckel U, Hessner MJ, Verbsky J: An autoinflammatory disease due to homozygous deletion of the IL1RN locus. The New England journal of medicine. 2009, 360 (23): 2438-2444.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang X, Bai HM, Chen L, Li B, Lu YC: Reduced expression of LC3B-II and Beclin 1 in glioblastoma multiforme indicates a down-regulated autophagic capacity that relates to the progression of astrocytic tumors. Journal of clinical neuroscience: official journal of the Neurosurgical Society of Australasia. 2010, 17 (12): 1515-1519.View ArticleGoogle Scholar
- Koukourakis MI, Giatromanolaki A, Sivridis E, Pitiakoudis M, Gatter KC, Harris AL: Beclin 1 over- and underexpression in colorectal cancer: distinct patterns relate to prognosis and tumour hypoxia. British journal of cancer. 2010, 103 (8): 1209-1214.PubMed CentralView ArticlePubMedGoogle Scholar
- Hagen L, Pena-Diaz J, Kavli B, Otterlei M, Slupphaug G, Krokan HE: Genomic uracil and human disease. Experimental cell research. 2006, 312 (14): 2666-2672.View ArticlePubMedGoogle Scholar
- Oceandy D, Cartwright EJ, Neyses L: Ras-association domain family member 1A (RASSF1A)-where the heart and cancer meet. Trends in cardiovascular medicine. 2009, 19 (8): 262-267.View ArticlePubMedGoogle Scholar
- Oceandy D, Pickard A, Prehar S, Zi M, Mohamed TM, Stanley PJ, Baudoin-Stanley F, Nadif R, Tommasi S, Pfeifer GP, Armesilla AL, Cartwright EJ, Neyses L: Tumor suppressor Ras-association domain family 1 isoform A is a novel regulator of cardiac hypertrophy. Circulation. 2009, 120 (7): 607-616.View ArticlePubMedGoogle Scholar
- Wilson HM: Feline alimentary lymphoma: demystifying the enigma. Top Companion Anim Med. 2008, 23 (4): 177-184.View ArticlePubMedGoogle Scholar
- MySQL Relational Database Server. [http://www.mysql.com/]
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402.PubMed CentralView ArticlePubMedGoogle Scholar
- Basic Local Alignment Search Tool. [http://blast.ncbi.nlm.nih.gov/Blast.cgi]
- Grasso C, Lee C: Combining partial order alignment and progressive multiple sequence alignment increases alignment speed and scalability to very large alignment problems. Bioinformatics. 2004, 20 (10): 1546-1556.View ArticlePubMedGoogle Scholar
- Lee C, Grasso C, Sharlow MF: Multiple sequence alignment using partial order graphs. Bioinformatics. 2002, 18 (3): 452-464.View ArticlePubMedGoogle Scholar
- Haider S, Ballester B, Smedley D, Zhang J, Rice P, Kasprzyk A: BioMart Central Portal--unified access to biological data. Nucleic Acids Res. 2009, W23-27. 37 Web ServerGoogle Scholar
- Yang Z: PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci. 1997, 13 (5): 555-556.PubMedGoogle Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24 (8): 1586-1591.View ArticlePubMedGoogle Scholar
- Suyama M, Torrents D, Bork P: PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 2006, W609-612. 34 Web ServerGoogle Scholar
- Barrell D, Dimmer E, Huntley RP, Binns D, O'Donovan C, Apweiler R: The GOA database in 2009--an integrated Gene Ontology Annotation resource. Nucleic Acids Res. 2009, D396-403. 37 DatabaseGoogle Scholar
- Ogata H, Goto S, Fujibuchi W, Kanehisa M: Computation with the KEGG pathway database. Biosystems. 1998, 47 (1-2): 119-128.View ArticlePubMedGoogle Scholar
- Mlecnik B, Scheideler M, Hackl H, Hartler J, Sanchez-Cabo F, Trajanoski Z: PathwayExplorer: web service for visualizing high-throughput expression data on biological pathways. Nucleic Acids Res. 2005, W633-637. 33 Web ServerGoogle Scholar
- Sherman BT, Huang da W, Tan Q, Guo Y, Bour S, Liu D, Stephens R, Baseler MW, Lane HC, Lempicki RA: DAVID Knowledgebase: a gene-centered database integrating heterogeneous gene annotation resources to facilitate high-throughput gene functional analysis. BMC Bioinformatics. 2007, 8: 426-PubMed CentralView ArticlePubMedGoogle Scholar
- Hamosh A, Scott AF, Amberger JS, Bocchini CA, McKusick VA: Online Mendelian Inheritance in Man (OMIM), a knowledgebase of human genes and genetic disorders. Nucleic Acids Res. 2005, D514-517. 33 DatabaseGoogle Scholar
- Wang Y, Kelly MA, Cowan TM, Longo N: A missense mutation in the OCTN2 gene associated with residual carnitine transport activity. Human mutation. 2000, 15 (3): 238-245.View ArticlePubMedGoogle Scholar
- D'Antonio KB, Toubaji A, Albadine R, Mondul AM, Platz EA, Netto GJ, Getzenberg RH: Extracellular matrix associated protein CYR61 is linked to prostate cancer development. The Journal of urology. 2010, 183 (4): 1604-1610.PubMed CentralView ArticlePubMedGoogle Scholar
- Arman A, Soylu O, Yildirim A, Furman A, Ercelen N, Aydogan H, Coker A, Tezel T: Interleukin-1 receptor antagonist gene VNTR polymorphism is associated with coronary artery disease. Arquivos brasileiros de cardiologia. 2008, 91 (5): 293-298.View ArticlePubMedGoogle Scholar
- Di Domenico F, Sultana R, Tiu GF, Scheff NN, Perluigi M, Cini C, Butterfield DA: Protein levels of heat shock proteins 27, 32, 60, 70, 90 and thioredoxin-1 in amnestic mild cognitive impairment: an investigation on the role of cellular stress response in the progression of Alzheimer disease. Brain research. 2010, 1333: 72-81.PubMed CentralView ArticlePubMedGoogle Scholar
- Lu CW, Lin JH, Rajawat YS, Jerng H, Rami TG, Sanchez X, DeFreitas G, Carabello B, DeMayo F, Kearney DL, Miller G, Li H, Pfaffinger PJ, Bowles NE, Khoury DS, Towbin JA, et al: Functional and clinical characterization of a mutation in KCNJ2 associated with Andersen-Tawil syndrome. Journal of medical genetics. 2006, 43 (8): 653-659.PubMed CentralView ArticlePubMedGoogle Scholar
- Priori SG, Pandit SV, Rivolta I, Berenfeld O, Ronchetti E, Dhamoon A, Napolitano C, Anumonwo J, di Barletta MR, Gudapakkam S, Bosi G, Stramba-Badiale M, Jalife J: A novel form of short QT syndrome (SQT3) is caused by a mutation in the KCNJ2 gene. Circulation research. 2005, 96 (7): 800-807.View ArticlePubMedGoogle Scholar
- Jacobson SG, Cideciyan AV, Aleman TS, Sumaroka A, Schwartz SB, Windsor EA, Roman AJ, Heon E, Stone EM, Thompson DA: RDH12 and RPE65, visual cycle genes causing leber congenital amaurosis, differ in disease expression. Investigative ophthalmology & visual science. 2007, 48 (1): 332-338.View ArticleGoogle Scholar
- Maruszak A, Safranow K, Gustaw K, Kijanowska-Haładyna B, Jakubowska K, Olszewska M, Styczyńska M, Berdyński M, Tysarowski A, Chlubek D, Siedlecki J, Barcikowska M, Zekanowski C: PIN1 gene variants in Alzheimer's disease. BMC medical genetics. 2009, 10: 115-PubMed CentralView ArticlePubMedGoogle Scholar
- Barbaric D, Byth K, Dalla-Pozza L, Byrne JA: Expression of tumor protein D52-like genes in childhood leukemia at diagnosis: clinical and sample considerations. Leukemia research. 2006, 30 (11): 1355-1363.View ArticlePubMedGoogle Scholar
- Wierzbicki PM, Adrych K, Kartanowicz D, Wypych J, Stanislawowski M, Zwolinska-Wcislo M, Celinski K, Skrodzka D, Godlewski J, Korybalski B, Smoczynski M, Kmiec Z: Overexpression of the fragile histidine triad (FHIT) gene in inflammatory bowel disease. Journal of physiology and pharmacology: an official journal of the Polish Physiological Society. 2009, 60 (Suppl 4): 57-62.Google Scholar
- Helmus Y, Denecke J, Yakubenia S, Robinson P, Luhn K, Watson DL, McGrogan PJ, Vestweber D, Marquardt T, Wild MK: Leukocyte adhesion deficiency II patients with a dual defect of the GDP-fucose transporter. Blood. 2006, 107 (10): 3959-3966.View ArticlePubMedGoogle Scholar
- Zhao X, Graves C, Ames SJ, Fisher DE, Spanjaard RA: Mechanism of regulation and suppression of melanoma invasiveness by novel retinoic acid receptor-gamma target gene carbohydrate sulfotransferase 10. Cancer research. 2009, 69 (12): 5218-5225.View ArticlePubMedGoogle Scholar
- Hashimoto T, Nakai M: Increased hippocampal quinone reductase 2 in Alzheimer's disease. Neurosci Lett. 2011, 502 (1): 10-12.View ArticlePubMedGoogle Scholar
- Cheng J, Zhu Y, He S, Lu Y, Chen J, Han B, Petrillo M, Wrzeszczynski KO, Yang S, Dai P, Zhai S, Han D, Zhang MQ, Li W, Liu X, Li H, Chen ZY, Yuan H: Functional mutation of SMAC/DIABLO, encoding a mitochondrial proapoptotic protein, causes human progressive hearing loss DFNA64. Am J Hum Genet. 2011, 89 (1): 56-66.PubMed CentralView ArticlePubMedGoogle Scholar
- Pluta P, Cebula-Obrzut B, Ehemann V, Pluta A, Wierzbowska A, Piekarski J, Bilski A, Nejc D, Kordek R, Robak T, Smolewski P, Jeziorski A: Correlation of Smac/DIABLO protein expression with the clinico-pathological features of breast cancer patients. Neoplasma. 2011, 58 (5): 430-435.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.