Patterns of molecular genetic variation among cat breeds

Genetic variation in cat breeds was assessed utilizing a panel of short tandem repeat (STR) loci genotyped in 38 cat breeds and 284 single-nucleotide polymorphisms (SNPs) genotyped in 24 breeds. Population structure in cat breeds generally reflects their recent ancestry and absence of strong breed barriers between some breeds. There is a wide range in the robustness of population definition, from breeds demonstrating high definition to breeds with as little as a third of their genetic variation partitioning into a single population. Utilizing the STRUCTURE algorithm, there was no clear demarcation of the number of population subdivisions; 16 breeds could not be resolved into independent populations, the consequence of outcrossing in established breeds to recently developed breeds with common ancestry. These 16 breeds were divided into 6 populations. Ninety-six percent of cats in a sample set of 1040 were correctly assigned to their classified breed or breed group/population. Average breed STR heterozygosities ranged from moderate (0.53; Havana, Korat) to high (0.85; Norwegian Forest Cat, Manx). Most of the variation in cat breeds was observed within a breed population (83.7%), versus 16.3% of the variation observed between populations. The hierarchical relationships of cat breeds is poorly defined as demonstrated by phylogenetic trees generated from both STR and SNP data, though phylogeographic grouping of breeds derived completely or in part from Southeast Asian ancestors was apparent.     

Assessing SNP markers for assigning individuals to cattle populations

The effectiveness of single nucleotide polymorphisms (SNPs) for the assignment of cattle to their source breeds was investigated by analysing a panel of 90 SNPs assayed on 24 European breeds. Breed assignment was performed by comparing the Bayesian and frequentist methods implemented in the structure 2.2 and geneclass 2 software programs. The use of SNPs for the reallocation of known individuals to their breeds of origin and the assignment of unknown individuals was tested. In the reallocation tests, the methods implemented in structure 2.2 performed better than those in geneclass 2 , with 96% vs. 85% correct assignments respectively. In contrast, the methods implemented in geneclass 2 showed a greater correct assignment rate in allocating animals treated as unknowns to a reference dataset (62% vs. 51% and 80% vs. 65% in field tests 1 and 2 respectively). These results demonstrate that SNPs are suitable for the assignment of individuals to reference breeds. The results also indicate that structure 2.2 and geneclass 2 can be complementary tools to assess breed integrity and assignment. Our findings also stress the importance of a high-quality reference dataset in allocation studies.     

Incidence of pyometra in Swedish insured cats

Pyometra is a clinically relevant problem in intact female cats and dogs. The etiology is similar in both animal species, with the disease caused by bacterial infection of a progesterone-sensitized uterus. Here, we studied pyometra in cats with the aim to describe the incidence and probability of developing pyometra based on age and breed. The data used were reimbursed claims for veterinary care insurance or life insurance claims or both in cats insured in a Swedish insurance database from 1999 to 2006. The mean incidence rate (IR) for pyometra was about 17 cats per 10,000 cat years at risk (CYAR). Cats with pyometra were diagnosed at a median age of 4 years and a significant breed effect was observed. The breed with the highest IR (433 cats per 10,000 CYAR) was the Sphynx, and other breeds with IR over 60 cats per 10,000 CYAR were Siberian cat, Ocicat, Korat, Siamese, Ragdoll, Maine coon, and Bengal. Pyometra was more commonly diagnosed with increasing age, with a marked increase in cats older than 7 years. The mean case fatality rate in all cats was 5.7%, which is slightly higher than corresponding reports in dogs of 3% to 4%. Geographical location (urban or rural) did not affect the risk of developing the disease. The present study provides information of incidence and probability of developing pyometra based on age, breed, and urban or rural geographical location. These data may be useful for designing cat breeding programs in high-risk breeds and for future studies of the genetic background of the disease.     

Not another type of potato: MC1R and the russet coloration of Burmese cats

The Burmese is a breed of domestic cat that originated in Southeast Asia and was further developed in the United States. Variants in melanocortin 1 receptor (MC1R) causes common coat colour phenotypes in a variety of mammalian species but only limited colour variation in the domestic cat. Known as the extension (E) locus, melanocortin 1 receptor (MC1R) interacts with the agouti locus to produce the eumelanin and pheomelanin pigments. Recently, a novel reddish coloration, which is termed russet, was identified in the Burmese cat breed. Because this russet Burmese coloration changes with aging, MC1R was suggested as candidate gene. The similar colouration in specific lineages of Norwegian Forest cat known as amber (e) (c.250G>A; p.Asp84Asn) was excluded for this Burmese phenotype. The complete 954‐bp coding region of MC1R was directly sequenced in russet Burmese and suspected carriers. A 3‐bp deletion (c.439_441del) associated with the deletion of a phenyalanine (p.Phe146del) in the protein sequence was identified. All russet coloured cats were homozygous for the variant, and all obligate carriers were heterozygous, confirming that the deletion segregated concordantly with colouring in Burmese cats from the New Zealand foundation lineage. The variant was not identified in 442 cats from 26 different breeds and random‐bred cats. Twenty‐six Burmese from the USA did not have the variant. This MC1R variant defines a unique coloration and the second breed‐specific MC1R variant in cats. The interactions of the two recessive feline MC1R alleles (E  >  e, e^r) is unknown.     

Investigation of the role of opsin gene polymorphism in generalized progressive retinal atrophies in dogs

The generalized progressive retinal atrophies (gPRAs) form a group of retinal degenerations of pedigree dogs and cats, which have a variety of genetic origins (mostly unknown). We have examined the opsin gene for polymorphisms in several breeds of pedigree dog suffering from distinct forms of gPRA, by methods including single-strand conformation polymorphism analysis, microsatellite analysis and direct sequencing. The breeds examined included the Tibetan terrier, the miniature schnauzer, the Irish setter, the miniature poodle, the Labrador retriever and the English cocker spaniel, as well as individuals from breeds in which PRA has not been described and of mixed breed. Individuals from each of the named breeds suffering from PRA were compared with clinically normal dogs. Two polymorphisms were found. One, segregating within the Tibetan terrier population, but not seen in other breeds, was a synonymous transition at nucleotide position 780 in exon 3. Inheritance of this polymorphism suggests that opsin is unlikely to contain mutations causative of gPRA in this breed. The other polymorphism occurred between all miniature schnauzers examined and dogs of other breeds. It consisted of a single base insertion in intron 2. No polymorphisms in the opsin sequence were detected in any other breed. DNA sequencing allowed rigorous exclusion of mutations in opsin as a cause of gPRA in miniature poodles, English cocker spaniels or Labrador retrievers.     

Genomic diversity and population structure of Mexican and Spanish bovine Lidia breed

The Lidia bovine breed is distinguished for its low genetic exchangeability given its selection on aggressive behavior, its management uniqueness and its subdivided structure. In this study, we present a comprehensive genome‐wide analysis of genetic diversity, population structure and admixture of 468 animals from Mexican and Spanish Lidia breed populations and 64 samples belonging to 10 Spanish native and American‐creole breeds using 37 148 single nucleotide polymorphisms. We found similar average inbreeding values in the Lidia breed, with different distributions within groups; variability of inbreeding values among Spanish lineages was significant and no differences were found among the Mexican sub‐populations. Together, the high F_IS of the lineages and the behavior of the runs of homozygosity are consequences of the lineage's small effective population sizes, contributing to their inbreeding increase. Population admixture analysis discarded any influence on the genetic structure of the Lidia populations from the Spanish native and American‐creole breeds. In addition, both Lidia populations depicted different genetic origins, with the exception of some Mexican individuals whose origins traced back to recent Spanish importations.     

Genetic epidemiology of blood type,disease and trait variants,and genome-wide genetic diversity in over 11,000 domestic cats

In the largest DNA-based study of domestic cats to date, 11,036 individuals (10,419 pedigreed cats and 617 non-pedigreed cats) were genotyped via commercial panel testing elucidating the distribution and frequency of known disease, blood type, and physical trait associated genetic variants across cat breeds. This study provides allele frequencies for many disease-associated variants for the first time and provides updates on previously reported information with evidence suggesting that DNA testing has been effectively used to reduce disease associated variants within certain pedigreed cat populations over time. We identified 13 disease-associated variants in 47 breeds or breed types in which the variant had not previously been documented, highlighting the relevance of comprehensive genetic screening across breeds. Three disease-associated variants were discovered in non-pedigreed cats only. To investigate the causality of nine disease-associated variants in cats of different breed backgrounds our veterinarians conducted owner interviews, reviewed clinical records, and invited cats to have follow-up clinical examinations. Additionally, genetic variants determining blood types A, B and AB, which are relevant clinically and in cat breeding, were genotyped. Appearance-associated genetic variation in all cats is also discussed. Lastly, genome-wide SNP heterozygosity levels were calculated to obtain a comparable measure of the genetic diversity in different cat breeds. This study represents the first comprehensive exploration of informative Mendelian variants in felines by screening over 10,000 pedigreed cats. The results qualitatively contribute to the understanding of feline variant heritage and genetic diversity and demonstrate the clinical utility and importance of such information in supporting breeding programs and the research community. The work also highlights the crucial commitment of pedigreed cat breeders and registries in supporting the establishment of large genomic databases, that when combined with phenotype information can advance scientific understanding and provide insights that can be applied to improve the health and welfare of cats.     

Novel Y‐chromosome short tandem repeats in Sus scrofa and their variation in European wild boar and domestic pig populations

Y‐chromosome markers are important tools for studying male‐specific gene flow within and between populations, hybridization patterns and kinship. However, their use in non‐human mammals is often hampered by the lack of Y‐specific polymorphic markers. We identified new male‐specific short tandem repeats (STRs) in Sus scrofa using the available genome sequence. We selected four polymorphic loci (5–10 alleles per locus), falling in one duplicated and two single‐copy regions. A total of 32 haplotypes were found by screening 211 individuals from eight wild boar populations across Europe and five domestic pig populations. European wild boar were characterized by significantly higher levels of haplotype diversity compared to European domestic pigs (H_D = 0.904 ± 0.011 and H_D = 0.491 ± 0.077 respectively). Relationships among STR haplotypes were investigated by combining them with single nucleotide polymorphisms at two linked genes (AMELY and UTY) in a network analysis. A differentiation between wild and domestic populations was observed (F_ST = 0.229), with commercial breeds sharing no Y haplotype with the sampled wild boar. Similarly, a certain degree of geographic differentiation was observed across Europe, with a number of local private haplotypes and high diversity in northern populations. The described Y‐chromosome markers can be useful to track male inheritance and gene flow in wild and domestic populations, promising to provide insights into evolutionary and population genetics in Sus scrofa.     

Identification of high utility SNPs for population assignment and traceability purposes in the pig using high-throughput sequencing

The objectives of this study were to develop breed-specific single nucleotide polymorphisms (SNPs) in five pig breeds sequenced with Illumina's Genome Analyzer and to investigate their usefulness for breed assignment purposes. DNA pools were prepared for Duroc, Landrace, Large White, Pietrain and Wild Boar. The total number of animals used for sequencing was 153. SNP discovery was performed by aligning the filtered reads against Build 7 of the pig genome. A total of 313,964 high confidence SNPs were identified and analysed for the presence of breed-specific SNPs (defined in this context as SNPs for which one of the alleles was detected in only one breed). There were 29,146 putative breed-specific SNPs identified, of which 4441 were included in the PorcineSNP60 beadchip. Upon re-examining the genotypes obtained using the beadchip, 193 SNPs were confirmed as being breed specific. These 193 SNPs were subsequently used to assign an additional 490 individuals from the same breeds, using the sequenced individuals as reference populations. In total, four breed assignment tests were performed. Results showed that for all methods tested 99% of the animals were correctly assigned, with an average probability of assignment of at least 99.2%, indicating the high utility of breed-specific markers for breed assignment and traceability. This study provides a blueprint for the way next-generation sequencing technologies can be used for the identification of breed-specific SNPs, as well as evidence that these SNPs may be a powerful tool for breed assignment and traceability of animal products to their breeds of origin.     

Genetic diversity and population structure of Sicilian sheep breeds using microsatellite markers

Genetic diversity studies in domestic animals aim at evaluating genetic variation within and across breeds mainly for conservation purposes. In Sicily, dairy sheep production represents an important resource for hilly and mountain areas economy. Their milk is used for the production of traditional raw milk cheeses, sometimes protected designation of origin (PDO) cheeses. In some cases, the quality of these products is linked to a specific breed, i.e. mono-breed labelled cheeses and it is therefore important to be able to distinguish the milk of a breed from that of others, in order to guarantee both the consumer and the breed itself. In order to investigate the genetic structure and to perform an assignment test, a total of 331 individuals (Barbaresca, BAR n = 57, Comisana, COM n = 65, Pinzirita, PIN n = 75, Sarda, SAR n = 64, and Valle del Belice, VDB n = 70) were analysed using a panel of 20 microsatellite markers. A total of 259 alleles were observed with average polymorphic information content equal to 0.76, showing that the microsatellites panel used was highly informative. Estimates of observed heterozygosity ranged from 0.65 in the BAR breed to 0.75 in the COM breed. The low value of genetic differentiation among breeds (F_st = 0.049) may indicate that these breeds are little differentiated probably due to common history and breeding practices. The low F_is and F_it values indicated low level of inbreeding within and among breeds. The unrooted neighbor-joining dendrogram obtained from the Reynold's genetic distances, and factorial correspondence analysis revealed a separation between BAR and the other sheep breeds. Recent migration rates were estimated, showing that four out of the five breeds have not received a significant proportion of migrants. Only for the PIN breed a recent introgression rate from the VDB breed (7.2%) was observed. The Bayesian assignment test showed that BAR and SAR breeds had a more definite genetic structure (proportion of assignment of 92% and 86.6%, respectively), whereas the lowest assignment value was found in the PIN breed (67.1%). Our results indicated high genetic variability, low inbreeding and low genetic differentiation, except for BAR breed, and were in accordance with geographical location, history, and breeding practices. The low robustness of the assignment test makes it unfeasible for traceability purposes, due to the high level of admixture, in particular for COM, PIN and VDB.     

Two missense mutations in exon 9 of caprine PRLR gene were associated with litter size

Guanzhong (n = 321) and Boer (n = 191) goat breeds were used to detect single nucleotide polymorphisms (SNPs) in the coding regions of the prolactin receptor (PRLR) gene by DNA sequencing and PCR‐RFLP. Two SNPs (c.1457G>A and c.1645G>A) were identified that caused amino acid variations p.Ser485Asn and p.Val548Met respectively. Statistical results indicated that the c.1457G>A and c.1645G>A SNPs were significantly associated with litter size in Boer and Guanzhong goat breeds. Further analysis revealed that combined genotype C4 (GGGG) and haplotype G‐G were better than the others for litter size in both goat breeds. These results might contribute to goat genetic resources and breeding.     

Use of linked loci as individuals or haplotypes for marker-assisted breed assignment

The objective of this study was to use simulation to evaluate the benefits of considering haplotypes of loci when linked single nucleotide polymorphisms are used for breed assignment. Three breeds of 10 000 females each were simulated under eight scenarios that differed according to the number of generations separating the breeds, size of breed founder populations and recombination rate between linked loci. Molecular genotypes consisted of 20 groups of three linked loci each. Breed assignment was performed in the final generation and was based on the frequency method. Haplotypes were reconstructed using the expectation–maximization algorithm. Accuracy of breed assignment was based on the frequency of correct breed assignment. Assignment accuracy increased as more genotypes (loci or haplotypes) were considered and more animals were used to estimate genotypic frequencies within breed. For most scenarios, use of haplotypes yielded equal or greater accuracies than when loci were considered independent. The advantage of haplotypes tended to increase as linkage disequilibrium between adjacent loci increased. The greatest advantage for using haplotypes was observed when recombination rate was low (0.001), breeds were separated by few generations (100), and a relatively large number of founder animals (110) was used to form new breeds. In this situation, 90% accuracy of breed assignment was achieved using nine to 14 haplotypes (i.e. 27–42 loci) depending on breed, vs. 39–57 individual loci.     

Analysis of the genetic variation of vascular endothelial growth factor gene in three Chinese indigenous cattle breeds

PCR–SSCP and DNA sequencing methods were employed to screen the genetic variation of vascular endothelial growth factor (VEGF) gene in 675 individuals belonging to three Chinese indigenous cattle breeds including Qinchuan (QC), Jiaxian Red (JX) and Nanyang (NY) breed. Three new single nucleotide polymorphisms (SNPs) (g.6765T > C ss130456744, g.6860A > G ss130456745, g.6893T > C ss130456746) were found. One SNP (g.6765T > C) was detected in intron II of VEGF gene in all three breeds and the other two SNPs (g.6860A > G, g.6893T > C) were in exon III of VEGF gene only in NY breed. Among them, two synonymous mutations of exon III were identified: CCA (Pro) > CCG (Pro) at position 65th amino acid (aa) and TGT (Cys) > TGC (Cys) at position 76th aa of VEGF(190aa) in NY breed. Our study revealed that NY breed exhibited the most abundant genetic diversity in VEGF gene within the three cattle breeds. Furthermore, JX cattle breed was more similar to QC breed than to NY breed. Our genetic data in the present study supported the hypothesis that the distribution pattern of Chinese indigenous cattle breeds was closely related to the geographical and climatic background again.     

Penalized classification for optimal statistical selection of markers from high-throughput genotyping: application in sheep breeds

The identification of individuals’ breed of origin has several practical applications in livestock and is useful in different biological contexts such as conservation genetics, breeding and authentication of animal products. In this paper, penalized multinomial regression was applied to identify the minimum number of single nucleotide polymorphisms (SNPs) from high-throughput genotyping data for individual assignment to dairy sheep breeds reared in Sicily. The combined use of penalized multinomial regression and stability selection reduced the number of SNPs required to 48. A final validation step on an independent population was carried out obtaining 100% correctly classified individuals. The results using independent analysis, such as admixture, F_st, principal component analysis and random forest, confirmed the ability of these methods in selecting distinctive markers. The identified SNPs may constitute a starting point for the development of a SNP based identification test as a tool for breed assignment and traceability of animal products.     

Study of polymorphisms in the promoter region of ovine β-lactoglobulin gene and phylogenetic analysis among the Valle del Belice breed and other sheep breeds considered as ancestors

The aim of this work was to sequence the promoter region of β-lactoglobulin (BLG) gene in four sheep breeds, in order to identify polymorphisms, infer and analyze haplotypes, and phylogenetic relationship among the Valle del Belice breed and the other three breeds considered as ancestors. Sequencing analysis and alignment of the obtained sequences showed the presence of 36 single nucleotide polymorphisms (SNPs) and one deletion. A total of 22 haplotypes found in “best” reconstruction were inferred considering the 37 polymorphic sites identified. Haplotypes were used for the reconstruction of a phylogenetic tree using the Neighbor-Joining algorithm. The number of polymorphisms identified showed high variability within breeds. Analysis of genetic diversity indexes showed that the Sarda breed presented the lowest nucleotide diversity, whereas the Comisana breed presented the highest one. Comparing the nucleotide diversity among breeds, the highest value was obtained between Valle del Belice and Pinzirita breeds, whereas the lowest one was between Valle del Belice and Sarda breeds. Considering that polymorphisms in the promoter region of BLG gene could have a functional role associated with milk composition, the lowest value of nucleotide diversity between Valle del Belice and Sarda breeds may be related to a higher similarity of milk composition of these two breeds compared to the others. Further analyses will be conducted in order to evaluate the possible correlation between the genetic diversity indexes and the BLG content in milk of our breeds.     

Evidence of selection signatures that shape the Persian cat breed

The Persian cat is mainly characterized by an extremely brachycephalic face as part of the standard body conformation. Despite the popularity, world-wide distribution, and economic importance of the Persian cat as a fancy breed, little is known about the genetics of their hallmark morphology, brachycephaly. Over 800 cats from different breeds including Persian, non-Persian breeds (Abyssinian, Cornish Rex, Bengal, La Perm, Norwegian Forest, Maine Coon, Manx, Oriental, and Siamese), and Persian-derived breeds (British Shorthair, Scottish Fold, Selkirk Rex) were genotyped with the Illumina 63 K feline DNA array. The experimental strategy was composed of three main steps: (i) the Persian dataset was screened for runs of homozygosity to find and select highly homozygous regions; (ii) selected Persian homozygous regions were evaluated for the difference of homozygosity between Persians and those considered non-Persian breeds, and, (iii) the Persian homozygous regions most divergent from the non-Persian breeds were investigated by haplotype analysis in the Persian-derived breeds. Four regions with high homozygosity (H > 0.7) were detected, each with an average length of 1 Mb. Three regions can be considered unique to the Persian breed, with a less conservative haplotype pattern in the Persian-derived breeds. Moreover, two genes, CHL1 and CNTN6 known to determine face shape modification in humans, reside in one of the identified regions and therefore are positional candidates for the brachycephalic face in Persians. In total, the homozygous regions contained several neuronal genes that could be involved in the Persian cat behavior and can provide new insights into cat domestication.     

Nucleotide variability and haplotype heterogeneity at the porcine fat mass and obesity‐associated (FTO) gene

A large number of studies have confirmed that variants within the fat mass and obesity‐associated (FTO) gene are associated with higher obesity risk in humans. We and others have shown that FTO polymorphisms are associated with fat deposition and related traits in several pig populations, thus confirming the role of this gene in fatness across species. However, some differences observed in different pig populations may be derived, at least in part, from genetic heterogeneity at this locus. Here, we characterise the nucleotide variability and haplotype diversity of the porcine FTO gene in breeds having different predispositions to fat deposition traits. We resequenced 4749 bp of coding and non‐coding regions of the porcine FTO gene in 44 pigs of eight different breeds and identified 27 single nucleotide polymorphisms (SNPs) and four insertions/deletions. A positive Tajima's D‐value (P < 0.10) obtained in Italian Duroc pigs may be compatible with putative balancing selection. From the sequenced pig panel, 20 haplotypes were inferred, some of which clustered according to the breed of origin (Meishan and Italian Duroc). Genetic heterogeneity at this locus could complicate the dissection of the effects of this gene on fat deposition and production traits in pigs. This situation resembles, to some extent, what has been reported in humans, thus making the study of the porcine FTO gene variability especially interesting, as it could be used as a model to understand the complex and elusive role of this gene in human obesity.     

Candidate gene analysis of osteochondrosis in Spanish Purebred horses

Equine osteochondrosis (OC) is a frequent developmental orthopaedic disease with high economic impact on the equine industry and may lead to premature retirement of the animal as a result of chronic pain and lameness. The genetic background of OC includes different genes affecting several locations; however, these genetic associations have been tested in only one or few populations, lacking the validation in others. The aim of this study was to identify the genetic determinants of OC in the Spanish Purebred horse breed. For that purpose, we used a candidate gene approach to study the association between loci previously implicated in the onset and development of OC in other breeds and different OC locations using radiographic data from 144 individuals belonging to the Spanish Purebred horse breed. Of the 48 polymorphisms analysed, three single nucleotide polymorphisms (SNPs) located in the FAF1, FCN3 and COL1A2 genes were found to be associated with different locations of OC lesions. These data contribute insights into the complex gene networks underlying the multifactorial disease OC, and the associated SNPs could be used in a marker‐assisted selection strategy to improve horse health, welfare and competitive lifespan.