Effective population size of an indigenous Swiss cattle breed estimated from linkage disequilibrium

Effective population size is an important parameter for the assessment of genetic diversity within a livestock population and its development over time. If pedigree information is not available, linkage disequilibrium (LD) analysis might offer an alternative perspective for the estimation of effective population size. In this study, 128 individuals of the Swiss Eringer breed were genotyped using the Illumina BovineSNP50 beadchip. We set bin size at 50 kb for LD analysis, assuming that LD for proximal single nucleotide polymorphism (SNP)-pairs reflects distant breeding history while LD from distal SNP-pairs would reflect near history. Recombination rates varied among different regions of the genome. The use of physical distances as an approximation of genetic distances (e.g. setting 1 Mb = 0.01 Morgan) led to an upward bias in LD-based estimates of effective population size for generations beyond 50, while estimates for recent history were unaffected. Correction for restricted sample size did not substantially affect these results. LD-based actual effective population size was estimated in the range of 87-149, whereas pedigree-based effective population size resulted in 321 individuals. For conservation purposes, requiring knowledge of recent history (<50 generations), approximation assuming constant recombination rate seemed adequate.

Pseudomyotonia, a muscle function disorder associated with an inherited ATP2A1 (SERCA1) defect in a Dutch Improved Red and White cross-breed calf

A Dutch Improved Red and White cross-breed heifer calf was evaluated for a muscular disorder resulting in exercise induced muscle stiffness. Clinical findings included generalized exercise-induced muscle spasms with normal response to muscle percussion. Electromyography showed no myotonic discharges, thus ruling out myotonia. Whereas histological examination of muscle tissue was unremarkable, Ca(2+)-ATPase activity of sarcoplasmatic reticulum membranes (SERCA1) was markedly decreased compared to control animals. Mutation analysis revealed the presence of a missense mutation in the ATP2A1 gene encoding the SERCA1 protein (p.Arg559Cys). The present case presents similarities to human Brody's disease, but also to pseudomyotonia and congenital muscular dystonia previously described in different cattle breeds.

Hexose, sialic acid and sialyllactose concentrations in the milk of dairy and non-dairy breed cows

It has been suggested that the ratio of lactose to milk oligosaccharides in mammalian milk/colostrum is based on the ratio of expression of a-lactalbumin and glycosyltransferases in the mammary epithelial cells. It has also been suggested that the high secretion of milk in dairy breed cows has been acquired by a high expression of a-lactalbumin expression. As there is a large difference of milk secretion level between dairy and non dairy breed cows, there may be a difference in the ratio of lactose to milk oligosaccharides in milks between dairy and non dairy breed cows. In this study, the concentrations of hexose, sialic acid as well as sialyllactoses, which are representative bovine milk oligosaccharides, were determined in the milks of dairy and non dairy breed cows. The concentration of hexose was significantly higher in the milks of non dairy breed cows than that of dairy breed cows, but there were no significant differences with respect to sialic acid and sialyllactose. The significant difference of the ratio of the concentrations of 3'- and 6'-sialyllactose to total hexose in milk was not observed between dairy and non dairy cows.

Analysis of pedigree and conformation data to explain genetic variability of the horse breed Franches-Montagnes

Franches-Montagnes is the only native horse breed in Switzerland, therefore special efforts should be made for ensuring its survival. The objectives of this study were to characterize the structure of this population as well as genetic variability with pedigree data, conformation traits and molecular markers. Studies were focused to clarify if this population is composed of a heavy- and a light-type subpopulation. Extended pedigree records of 3-year-old stallions (n = 68) and mares (n = 108) were available. Evaluations of body conformation traits as well as pedigree data and molecular markers did not support the two-subpopulation hypothesis. The generation interval ranged from 7.8 to 9.3 years. The complete generation equivalent was high (>12). The number of effective ancestors varied between 18.9 and 20.1, whereof 50% of the genetic variability was attributed to seven of them. Genetic contribution of Warmblood horses ranged from 36% to 42% and that of Coldblood horses from 4% to 6%. The average inbreeding coefficient reached 6%. Inbreeding effective population size was 114.5 when the average increase of the inbreeding coefficient per year since 1910 was taken. Our results suggest that bottleneck situations occurred because of selection of a small number of sire lines. Promotion of planned matings between parents that are less related is recommended in order to avoid a reduction of the genetic diversity.

Genetic diversity among horse populations with a special focus on the Franches-Montagnes breed

Genetic characterization helps to assure breed integrity and to assign individuals to defined populations. The objective of this study was to characterize genetic diversity in six horse breeds and to analyse the population structure of the Franches-Montagnes breed, especially with regard to the degree of introgression with Warmblood. A total of 402 alleles from 50 microsatellite loci were used. The average number of alleles per locus was significantly lower in Thoroughbreds and Arabians. Average heterozygosities between breeds ranged from 0.61 to 0.72. The overall average of the coefficient of gene differentiation because of breed differences was 0.100, with a range of 0.036-0.263. No significant correlation was found between this parameter and the number of alleles per locus. An increase in the number of homozygous loci with increasing inbreeding could not be shown for the Franches-Montagnes horses. The proportion of shared alleles, combined with the neighbour-joining method, defined clusters for Icelandic Horse, Comtois, Arabians and Franches-Montagnes. A more disparate clustering could be seen for European Warmbloods and Thoroughbreds, presumably from frequent grading-up of Warmbloods with Thoroughbreds. Grading-up effects were also observed when Bayesian and Monte Carlo resampling approaches were used for individual assignment to a given population. Individual breed assignments to defined reference populations will be very difficult when introgression has occurred. The Bayesian approach within the Franches-Montagnes breed differentiated individuals with varied proportions of Warmblood.

Do common genetic variants in endotoxin signaling pathway contribute to predisposition to alcoholic liver cirrhosis?

BACKGROUND: Tumor necrosis factor-alpha (TNF-alpha) and interleukin-1beta (IL-1beta), produced by endotoxin-activated Kupffer cells, play a key role in the pathogenesis of alcoholic liver cirrhosis (ALC). Alleles TNFA -238A, IL1B -31T and variant IL1RN*2 of repeat polymorphism in the gene encoding the IL-1 receptor antagonist increase production of TNF-alpha and IL-1beta, respectively. Alleles CD14 -159T, TLR4 c.896G and TLR4 c.1196T modify activation of Kupffer cells by endotoxin. We confirmed the published associations between these common variants and genetic predisposition to ALC by means of a large case-control association study conducted on two Central European populations. METHODS: The study population comprised a Czech sample of 198 ALC patients and 370 controls (MONICA project), and a German sample of 173 ALC patients and 331 controls (KORA-Augsburg), and 109 heavy drinkers without liver disease. RESULTS: Single locus analysis revealed no significant difference between patients and controls in all tested loci. Diplotype [IL1RN 2/ 2; IL1B -31T+] was associated with increased risk of ALC in the pilot study, but not in the validation samples. CONCLUSIONS: Although cytokine mediated immune reactions play a role in the pathogenesis of ALC, hereditary susceptibility caused by variants in the corresponding genes is low in Central European populations.

CT scan based determination of optimal bone corridor for atlantoaxial ventral screw fixation in miniature breed dogs

OBJECTIVE: To describe the most reliable insertion angle, corridor length and width to place a ventral transarticular atlantoaxial screw in miniature breed dogs. STUDY DESIGN: Retrospective CT imaging study. SAMPLE POPULATION: Cervical CT scans of toy breed dogs (n = 21). METHODS: Dogs were divided into 2 groups--group 1: no atlantoaxial abnormalities; group 2: atlantoaxial instability. Insertion angle in medial to lateral and ventral to dorsal direction was measured in group 1. Corridor length and width were measured in groups 1 and 2. Corridor width was measured at 3 points of the corridor. Each variable was measured 3 times and the mean used for statistical analysis. RESULTS: Mean +/- SD optimal transarticular atlantoaxial insertion angle was determined to be 40 +/- 1 degrees in medial to lateral direction from the midline and 20 +/- 1 degrees in ventral to dorsal direction from the floor of the neural canal of C2. Mean corridor length was 7 mm (range, 4.5-8.0 mm). Significant correlation was found between corridor length, body weight, and age. Mean bone corridor width ranged from 3 to 5 mm. Statistically significant differences were found between individuals, gender and measured side. CONCLUSIONS: Optimal placement of a transarticular screw for atlantoaxial joint stabilization is very demanding because the screw path corridor is very narrow.

Genetic predisposition to prostate cancer

Prostate cancer (PC) is a significant economic and health burden in the U.S. and Europe but its causes are largely unknown. The most significant risk factors (after gender) are age and family history of the disease. A gene with high penetrance but low frequency on chromosome 1q, HPC 1, has been suggested to cause a proportion of the familial aggregation of PC but other more common genes, conferring less risk, are also thought to contribute to disease predisposition. We have pursued a strategy to study both types of genetic risk in PC. To identify high penetrance genes, affected men from thirteen families have been genotyped for genetic linkage analysis at six microsatellite markers spanning 45 cM of 1q24-25. Both LOD score and non-parametric statistics provide no significant support for HPC1 in this genomic region, although 3 of the families did combine to produce a LOD score of 0.9. These families will be included in a genome wide search for other PC predisposition genes as part of a multinational collaboration.^ For study of common genetic factors in PC development, leukocyte DNA samples from an unselected series of 55 patients and 67 controls have been examined for genetic differences in two other candidate genes, the androgen receptor gene, hAR, at Xq11-12, and the vitamin D receptor gene, hVDR, at 12q12-14. hAR was typed for two trinucleotide repeat length polymorphisms, (CAG)$\rm\sb{n}$ and (GGC)$\rm\sb{n},$ encoding polyglutamine and polyglycine tracts, respectively, which have been implicated in PC susceptibility. These data, combined with similarly processed patients and controls from the U.K. show no consistent association of allele length with PC risk. A novel finding, however, has been a significant association between the number of GGC repeats and the length of time between diagnosis and relapse in stage T1-T4 Caucasian patients irrespective of therapy and age of the patient. Of 49 patients who relapsed out of 108 entering the study, those with 16 or fewer GGC repeats had an average relapse-free-period of 101 (+/$-$7.7) months while for those with more than 16 repeats the period averaged 48 (+/$-$2.9) months, a difference of 2.1 fold or 4.4 years.^ The second gene, hVDR, was genotyped at two polymorphisms, a synonymous C/T substitution in exon 9 identified by differential TaqI enzymatic digestion and a variable length polyA tract in the 3$\sp\prime$ UTR. Although these polymorphisms are in strong linkage disequilibrium only the polyA region showed a possible association with PC risk. Men homozygous for alleles with fewer than 18 A's had an increased risk (OR = 3.0, p = 0.0578) compared to controls. This result is opposite to the findings of others and may either indicate off-setting random errors which together balance out to no significant overall effect or reflect more complex genetic and/or environmental associations.^ Overall, this research suggests that single gene familial predisposition may be less prominent in PC than in other cancers and that the characteristics of PC pathology may be useful in identifying the effects of common genetic factors. ^

NGS Nominated CELA1, HSPG2, and KCNK5 as Candidate Genes for Predisposition to Balkan Endemic Nephropathy

Balkan endemic nephropathy (BEN) is a familial chronic tubulointerstitial disease with insidious onset and slow progression leading to terminal renal failure. The results of molecular biological investigations propose that BEN is a multifactorial disease with genetic predisposition to environmental risk agents. Exome sequencing of 22 000 genes with Illumina Nextera Exome Enrichment Kit was performed on 22 DNA samples (11 Bulgarian patients and 11 Serbian patients). Software analysis was performed via NextGene, Provean, and PolyPhen. The frequency of all annotated genetic variants with deleterious/damaging effect was compared with those of European populations. Then we focused on nonannotated variants (with no data available about them and not found in healthy Bulgarian controls). There is no statistically significant difference between annotated variants in BEN patients and European populations. From nonannotated variants with more than 40% frequency in both patients' groups, we nominated 3 genes with possible deleterious/damaging variants-CELA1, HSPG2, and KCNK5. Mutant genes (CELA1, HSPG2, and KCNK5) in BEN patients encode proteins involved in basement membrane/extracellular matrix and vascular tone, tightly connected to process of angiogenesis. We suggest that an abnormal process of angiogenesis plays a key role in the molecular pathogenesis of BEN.

Doppler sonography of the medial arterial blood supply to the coxofemoral joints of 36 medium to large breed dogs and its relationship with radiographic signs of joint disease.

The medial arterial supply to 68 of the 72 coxofemoral joints of 36 medium to large breed dogs was examined ultrasonographically. The medial circumflex femoral artery and three branches were identified; the artery and its transverse branch were identified in all 68 joints, and the deep branch was identified in 61 joints, and the ascending branch was identified in 63. However, the acetabular and obturator branches were not identified. The pulsatility index, the mean velocity and the peak systolic velocity of the medial circumflex femoral artery were determined and associated with a radiographic score of degenerative coxofemoral joint disease and a lath distraction index (LDI). In joints with a LDI greater than 0.35, the pulsatility index was significantly lower (P=0.023) and its mean velocity was higher (P=0.005). However, no significant associations were observed in individual dogs when the measurements in both joints were taken into account.

Genome-Wide Analyses Suggest Mechanisms Involving Early B-Cell Development in Canine IgA Deficiency.

Immunoglobulin A deficiency (IgAD) is the most common primary immune deficiency disorder in both humans and dogs, characterized by recurrent mucosal tract infections and a predisposition for allergic and other immune mediated diseases. In several dog breeds, low IgA levels have been observed at a high frequency and with a clinical resemblance to human IgAD. In this study, we used genome-wide association studies (GWAS) to identify genomic regions associated with low IgA levels in dogs as a comparative model for human IgAD. We used a novel percentile groups-approach to establish breed-specific cut-offs and to perform analyses in a close to continuous manner. GWAS performed in four breeds prone to low IgA levels (German shepherd, Golden retriever, Labrador retriever and Shar-Pei) identified 35 genomic loci suggestively associated (p <0.0005) to IgA levels. In German shepherd, three genomic regions (candidate genes include KIRREL3 and SERPINA9) were genome-wide significantly associated (p <0.0002) with IgA levels. A ~20kb long haplotype on CFA28, significantly associated (p = 0.0005) to IgA levels in Shar-Pei, was positioned within the first intron of the gene SLIT1. Both KIRREL3 and SLIT1 are highly expressed in the central nervous system and in bone marrow and are potentially important during B-cell development. SERPINA9 expression is restricted to B-cells and peaks at the time-point when B-cells proliferate into antibody-producing plasma cells. The suggestively associated regions were enriched for genes in Gene Ontology gene sets involving inflammation and early immune cell development.

Genetic analysis of white facial and leg markings in the Swiss Franches-Montagnes Horse Breed

White markings and spotting patterns in animal species are thought to be a result of the domestication process. They often serve for the identification of individuals but sometimes are accompanied by complex pathological syndromes. In the Swiss Franches-Montagnes horse population, white markings increased vastly in size and occurrence during the past 30 years, although the breeding goal demands a horse with as little depigmented areas as possible. In order to improve selection and avoid more excessive depigmentation on the population level, we estimated population parameters and breeding values for white head and anterior and posterior leg markings. Heritabilities and genetic correlations for the traits were high (h(2) > 0.5). A strong positive correlation was found between the chestnut allele at the melanocortin-1-receptor gene locus and the extent of white markings. Segregation analysis revealed that our data fit best to a model including a polygenic effect and a biallelic locus with a dominant-recessive mode of inheritance. The recessive allele was found to be the white trait-increasing allele. Multilocus linkage disequilibrium analysis allowed the mapping of the putative major locus to a chromosomal region on ECA3q harboring the KIT gene.

Imputation of sequence level genotypes in the Franches-Montagnes horse breed

BACKGROUND A cost-effective strategy to increase the density of available markers within a population is to sequence a small proportion of the population and impute whole-genome sequence data for the remaining population. Increased densities of typed markers are advantageous for genome-wide association studies (GWAS) and genomic predictions. METHODS We obtained genotypes for 54 602 SNPs (single nucleotide polymorphisms) in 1077 Franches-Montagnes (FM) horses and Illumina paired-end whole-genome sequencing data for 30 FM horses and 14 Warmblood horses. After variant calling, the sequence-derived SNP genotypes (~13 million SNPs) were used for genotype imputation with the software programs Beagle, Impute2 and FImpute. RESULTS The mean imputation accuracy of FM horses using Impute2 was 92.0%. Imputation accuracy using Beagle and FImpute was 74.3% and 77.2%, respectively. In addition, for Impute2 we determined the imputation accuracy of all individual horses in the validation population, which ranged from 85.7% to 99.8%. The subsequent inclusion of Warmblood sequence data further increased the correlation between true and imputed genotypes for most horses, especially for horses with a high level of admixture. The final imputation accuracy of the horses ranged from 91.2% to 99.5%. CONCLUSIONS Using Impute2, the imputation accuracy was higher than 91% for all horses in the validation population, which indicates that direct imputation of 50k SNP-chip data to sequence level genotypes is feasible in the FM population. The individual imputation accuracy depended mainly on the applied software and the level of admixture.