Intestinal Tritrichomonas foetus infection in cats in Switzerland detected by in vitro cultivation and PCR

Tritrichomonas foetus, a parasite well known for its significance as venereally transmitted pathogen in cattle, has recently been identified as a cause of chronic large-bowel diarrhea in domestic cats in the US, UK, and, more recently, also in Norway. In a period of 3 months (October to December 2007), 45 cats of Switzerland suffering from chronic diarrhea were investigated for intestinal infections, including a search for trichomonads. A commercially available in vitro culture system was used to screen for infection, complemented with a PCR and subsequent amplicon sequencing to support speciation. The PCR is based upon amplification of a sequence derived from the internal transcribed spacer region 1 (ITS1) on the ribosomal RNA gene (rRNA) using primers designed to detect a broad range of genera and species belonging to the family of Trichomonadidae. The method was furthermore adapted to the uracil DNA glycosylase (UDG) system in order to prevent carry-over contamination and it included a recombinant internal control to track for inhibitory reactions. Eleven out of the 45 cats were culture-positive, as revealed by microscopic identification of trichomonadid organisms. One of the isolates was subjected to scanning electron microscopy and findings revealed the presence of three flagella, thus placing the isolate into the gender Tritrichomonas sp. PCR and subsequent amplicon sequencing were carried out with ten of the 11 isolates. A total homology with published T. foetus sequences was confirmed in all of the cases. T. foetus therefore appears to range among those organisms that can cause chronic diarrhea in cats in Switzerland.

Multiplex strategy for multilocus sequence typing, fla typing, and genetic determination of antimicrobial resistance of Campylobacter jejuni and Campylobacter coli isolates collected in Switzerland

We present an optimized multilocus sequence typing (MLST) scheme with universal primer sets for amplifying and sequencing the seven target genes of Campylobacter jejuni and Campylobacter coli. Typing was expanded by sequence determination of the genes flaA and flaB using optimized primer sets. This approach is compatible with the MLST and flaA schemes used in the PubMLST database and results in an additional typing method using the flaB gene sequence. An identification module based on the 16S rRNA and rpoB genes was included, as well as the genetic determination of macrolide and quinolone resistances based on mutations in the 23S rRNA and gyrA genes. Experimental procedures were simplified by multiplex PCR of the 13 target genes. This comprehensive approach was evaluated with C. jejuni and C. coli isolates collected in Switzerland. MLST of 329 strains resulted in 72 sequence types (STs) among the 186 C. jejuni strains and 39 STs for the 143 C. coli isolates. Fourteen (19%) of the C. jejuni and 20 (51%) of the C. coli STs had not been found previously. In total, 35% of the C. coli strains collected in Switzerland contained mutations conferring antibiotic resistance only to quinolone, 15% contained mutations conferring resistance only to macrolides, and 6% contained mutations conferring resistance to both classes of antibiotics. In C. jejuni, these values were 31% and 0% for quinolone and macrolide resistance, respectively. The rpoB sequence allowed phylogenetic differentiation between C. coli and C. jejuni, which was not possible by 16S rRNA gene analysis. An online Integrated Database Network System (SmartGene, Zug, Switzerland)-based platform for MLST data analysis specific to Campylobacter was implemented. This Web-based platform allowed automated allele and ST designation, as well as epidemiological analysis of data, thus streamlining and facilitating the analysis workflow. Data networking facilitates the exchange of information between collaborating centers. The described approach simplifies and improves the genotyping of Campylobacter, allowing cost- and time-efficient routine monitoring.

Occurrence, quantification, and genotyping of Mycoplasma conjunctivae in wild Caprinae with and without infectious keratoconjunctivitis

Mycoplasma conjunctivae, the causative agent of infectious keratoconjunctivitis (IKC), was recently detected in asymptomatic Alpine ibex (Capra ibex ibex). This suggested that an external source of infection may not be required for an IKC outbreak in wildlife but might be initiated by healthy carriers, which contradicted previous serologic investigations in chamois. Our aims were to 1) assess the prevalence of M. conjunctivae among asymptomatic ibex and Alpine chamois (Rupicapra rupicapra rupicapra) and its frequency in IKC-affected animals, 2) determine mycoplasma loads in different disease stages, and 3) characterize the M. conjunctivae strains involved. Eye swabs from 654 asymptomatic and 204 symptomatic animals were collected in diverse Swiss regions between 2008 and 2010, and tested by TaqMan real-time PCR. Data analysis was performed considering various patterns of IKC occurrence in the respective sampling regions. Strains from 24 animals were compared by cluster analysis. Prevalence of M. conjunctivae was 5.6% (95% confidence interval [CI]: 3.7-8.1%) in asymptomatic ibex and 5.8% (CI: 3.0-9.9%) in asymptomatic chamois, with significant differences between years and regions in both species. Detection frequency in symptomatic animals was significantly higher during IKC outbreaks than in nonepidemic situations (i.e., regular but low incidence or sporadic occurrence). Mycoplasma load was significantly lower in eyes from healthy carriers and animals with mild signs than from animals with moderate and severe signs. Although some strains were found in both asymptomatic and diseased animals of the same species, others apparently differed in their pathogenic potential depending on the infected species. Overall, we found a widespread occurrence of M. conjunctivae in wild Caprinae with and without IKC signs. Our results confirm the central role of M. conjunctivae in outbreaks but suggest that other infectious agents may be involved in IKC cases in nonepidemic situations. Additionally, presence and severity of signs are related to the quantity of M. conjunctivae in the eyes rather than to the strain. We propose that individual or environmental factors influence the clinical expression of the disease and that persistence of M. conjunctivae in populations of wild Caprinae cannot be excluded.

Unique mixed phenotype and unexpected functional effect revealed by novel compound heterozygosity mutations involving SCN5A.

BACKGROUND Functional characterization of mutations involving the SCN5A-encoded cardiac sodium channel has established the pathogenic mechanisms for type 3 long QT syndrome and type 1 Brugada syndrome and has provided key insights into the physiological importance of essential structure-function domains. OBJECTIVE This study sought to present the clinical and biophysical phenotypes discerned from compound heterozygosity mutations in SCN5A on different alleles in a toddler diagnosed with QT prolongation and fever-induced ventricular arrhythmias. METHODS A 22-month-old boy presented emergently with fever and refractory ventricular tachycardia. Despite restoration of sinus rhythm, the infant sustained profound neurological injury and died. Using polymerase chain reaction, denaturing high-performance liquid chromatography, and direct DNA sequencing, comprehensive open-reading frame/splice mutational analysis of the 12 known long QT syndrome susceptibility genes was performed. RESULTS The infant had 2 SCN5A mutations: a maternally inherited N-terminal frame shift/deletion (R34fs/60) and a paternally inherited missense mutation, R1195H. The mutations were engineered by site-directed mutagenesis and heterologously expressed transiently in HEK293 cells. As expected, the frame-shifted and prematurely truncated peptide, SCN5A-R34fs/60, showed no current. SCN5A-R1195H had normal peak and late current but abnormal voltage-dependent gating parameters. Surprisingly, co-expression of SCN5A-R34fs/60 with SCN5A-R1195H elicited a significant increase in late sodium current, whereas co-expression of SCN5A-WT with SCN5A-R34fs/60 did not. CONCLUSIONS A severe clinical phenotype characterized by fever-induced monomorphic ventricular tachycardia and QT interval prolongation emerged in a toddler with compound heterozygosity involving SCN5A: R34fs/60, and R1195H. Unexpectedly, the 94-amino-acid fusion peptide derived from the R34fs/60 mutation accentuated the late sodium current of R1195H-containing Na(V)1.5 channels in vitro.

Genome-wide analysis in German shepherd dogs reveals association of a locus on CFA 27 with atopic dermatitis

Humans and dogs are both affected by the allergic skin disease atopic dermatitis (AD), caused by an interaction between genetic and environmental factors. The German shepherd dog (GSD) is a high-risk breed for canine AD (CAD). In this study, we used a Swedish cohort of GSDs as a model for human AD. Serum IgA levels are known to be lower in GSDs compared to other breeds. We detected significantly lower IgA levels in the CAD cases compared to controls (p = 1.1 × 10(-5)) in our study population. We also detected a separation within the GSD cohort, where dogs could be grouped into two different subpopulations. Disease prevalence differed significantly between the subpopulations contributing to population stratification (λ = 1.3), which was successfully corrected for using a mixed model approach. A genome-wide association analysis of CAD was performed (n cases = 91, n controls = 88). IgA levels were included in the model, due to the high correlation between CAD and low IgA levels. In addition, we detected a correlation between IgA levels and the age at the time of sampling (corr = 0.42, p = 3.0 × 10(-9)), thus age was included in the model. A genome-wide significant association was detected on chromosome 27 (praw = 3.1 × 10(-7), pgenome = 0.03). The total associated region was defined as a ~1.5-Mb-long haplotype including eight genes. Through targeted re-sequencing and additional genotyping of a subset of identified SNPs, we defined 11 smaller haplotype blocks within the associated region. Two blocks showed the strongest association to CAD. The ~209-kb region, defined by the two blocks, harbors only the PKP2 gene, encoding Plakophilin 2 expressed in the desmosomes and important for skin structure. Our results may yield further insight into the genetics behind both canine and human AD.

A multiplex PCR for the simultaneous detection and genotyping of the Echinococcus granulosus complex.

Echinococcus granulosus is characterized by high intra-specific variability (genotypes G1-G10) and according to the new molecular phylogeny of the genus Echinococcus, the E. granulosus complex has been divided into E. granulosus sensu stricto (G1-G3), E. equinus (G4), E. ortleppi (G5), and E. canadensis (G6-G10). The molecular characterization of E. granulosus isolates is fundamental to understand the spatio-temporal epidemiology of this complex in many endemic areas with the simultaneous occurrence of different Echinococcus species and genotypes. To simplify the genotyping of the E. granulosus complex we developed a single-tube multiplex PCR (mPCR) allowing three levels of discrimination: (i) Echinococcus genus, (ii) E. granulosus complex in common, and (iii) the specific genotype within the E. granulosus complex. The methodology was established with known DNA samples of the different strains/genotypes, confirmed on 42 already genotyped samples (Spain: 22 and Bulgaria: 20) and then successfully applied on 153 unknown samples (Tunisia: 114, Algeria: 26 and Argentina: 13). The sensitivity threshold of the mPCR was found to be 5 ng Echinoccoccus DNA in a mixture of up to 1 µg of foreign DNA and the specificity was 100% when template DNA from closely related members of the genus Taenia was used. Additionally to DNA samples, the mPCR can be carried out directly on boiled hydatid fluid or on alkaline-lysed frozen or fixed protoscoleces, thus avoiding classical DNA extractions. However, when using Echinococcus eggs obtained from fecal samples of infected dogs, the sensitivity of the mPCR was low (<40%). Thus, except for copro analysis, the mPCR described here has a high potential for a worldwide application in large-scale molecular epidemiological studies on the Echinococcus genus.

Learner preferences regarding integrating, sequencing and aligning virtual patients with other activities in the undergraduate medical curriculum: A focus group study

CONTEXT: E-learning resources, such as virtual patients (VPs), can be more effective when they are integrated in the curriculum. To gain insights that can inform guidelines for the curricular integration of VPs, we explored students' perceptions of scenarios with integrated and non-integrated VPs aimed at promoting clinical reasoning skills. METHODS: During their paediatric clerkship, 116 fifth-year medical students were given at least ten VPs embedded in eight integrated scenarios and as non-integrated add-ons. The scenarios differed in the sequencing and alignment of VPs and related educational activities, tutor involvement, number of VPs, relevance to assessment and involvement of real patients. We sought students' perceptions on the VP scenarios in focus group interviews with eight groups of 4-7 randomly selected students (n = 39). The interviews were recorded, transcribed and analysed qualitatively. RESULTS: The analysis resulted in six themes reflecting students' perceptions of important features for effective curricular integration of VPs: (i) continuous and stable online access, (ii) increasing complexity, adapted to students' knowledge, (iii) VP-related workload offset by elimination of other activities, (iv) optimal sequencing (e.g.: lecture--1 to 2 VP(s)--tutor-led small group discussion--real patient) and (V) optimal alignment of VPs and educational activities, (vi) inclusion of VP topics in assessment. CONCLUSIONS: The themes appear to offer starting points for the development of a framework to guide the curricular integration of VPs. Their impact needs to be confirmed by studies using quantitative controlled designs.

Next-generation sequencing of HIV-1 RNA genomes: determination of error rates and minimizing artificial recombination.

Next-generation sequencing (NGS) is a valuable tool for the detection and quantification of HIV-1 variants in vivo. However, these technologies require detailed characterization and control of artificially induced errors to be applicable for accurate haplotype reconstruction. To investigate the occurrence of substitutions, insertions, and deletions at the individual steps of RT-PCR and NGS, 454 pyrosequencing was performed on amplified and non-amplified HIV-1 genomes. Artificial recombination was explored by mixing five different HIV-1 clonal strains (5-virus-mix) and applying different RT-PCR conditions followed by 454 pyrosequencing. Error rates ranged from 0.04-0.66% and were similar in amplified and non-amplified samples. Discrepancies were observed between forward and reverse reads, indicating that most errors were introduced during the pyrosequencing step. Using the 5-virus-mix, non-optimized, standard RT-PCR conditions introduced artificial recombinants in a fraction of at least 30% of the reads that subsequently led to an underestimation of true haplotype frequencies. We minimized the fraction of recombinants down to 0.9-2.6% by optimized, artifact-reducing RT-PCR conditions. This approach enabled correct haplotype reconstruction and frequency estimations consistent with reference data obtained by single genome amplification. RT-PCR conditions are crucial for correct frequency estimation and analysis of haplotypes in heterogeneous virus populations. We developed an RT-PCR procedure to generate NGS data useful for reliable haplotype reconstruction and quantification.

Sequencing of TGF-beta pathway genes in familial cases of intracranial aneurysm.

BACKGROUND AND PURPOSE: Familial aggregation of intracranial aneurysms (IA) strongly suggests a genetic contribution to pathogenesis. However, genetic risk factors have yet to be defined. For families affected by aortic aneurysms, specific gene variants have been identified, many affecting the receptors to transforming growth factor-beta (TGF-beta). In recent work, we found that aortic and intracranial aneurysms may share a common genetic basis in some families. We hypothesized, therefore, that mutations in TGF-beta receptors might also play a role in IA pathogenesis. METHODS: To identify genetic variants in TGF-beta and its receptors, TGFB1, TGFBR1, TGFBR2, ACVR1, TGFBR3, and ENG were directly sequenced in 44 unrelated patients with familial IA. Novel variants were confirmed by restriction digestion analyses, and allele frequencies were analyzed in cases versus individuals without known intracranial disease. Similarly, allele frequencies of a subset of known SNPs in each gene were also analyzed for association with IA. RESULTS: No mutations were found in TGFB1, TGFBR1, TGFBR2, or ACVR1. Novel variants identified in ENG (p.A60E) and TGFBR3 (p.W112R) were not detected in at least 892 reference chromosomes. ENG p.A60E showed significant association with familial IA in case-control studies (P=0.0080). No association with IA could be found for any of the known polymorphisms tested. CONCLUSIONS: Mutations in TGF-beta receptor genes are not a major cause of IA. However, we identified rare variants in ENG and TGFBR3 that may be important for IA pathogenesis in a subset of families.

Characterization of osteopontin charge forms produced by osteoblasts

Osteopontin (OPN) is a highly-phosphorylated extracellular matrix protein localized in bone, kidney, placenta, T-lymphocytes, macrophages, smooth muscle of the vascular system, milk, urine, and plasma. In ROS 17/2.8 osteoblast-like osteosarcoma cells, 1,25-dihydroxyvitamin D3 [1,25(OH)2D 3] regulates OPN at the transcriptional level resulting in increased steady state mRNA levels and increased production of OPN protein, maximal at 48 hours. Using ROS 17/2.8 cells as an osteoblast model, OPN was purified from culture medium after three hour treatments of either vehicle (ethanol) or 1,25(OH)2D3 via barium citrate precipitation followed by immunoaffinity chromatography. ^ Here, further evidence of regulation of OPN by 1,25(OH)2D 3 at the posttranslational level is presented. Prior to the up-regulation of OPN at the transcriptional level, 1,25(OH)2D3 induces a shift in OPN isoelectric point (pI) detected on two-dimensional gels from pI 4.6 to pI 5.1. Loading equal amounts of [32P]-labeled OPN recovered from ROS 17/2.8 cells exposed to 1,25(OH)2D3 or vehicle alone for three hours reveals that the shift from pI 4.6 to 5.1 is the result of reduced phosphorylation. Using structural analogs to 1,25(OH) 2D3, analog AT [25-(OH)-16-ene-23-yne-D3], which triggers Ca2+ influx through voltage sensitive Ca2+ channels but does not bind to the vitamin D receptor, mimicked the OPN pI shift while analog BT [1,25(OH)2-22-ene-24-cyclopropyl-D 3], which binds to the vitamin D receptor but does not allow Ca 2+ influx, did not. Inclusion of the Ca2+ channel blocker nifedipine also blocks the charge shift conversion of OPN. Further analysis of the signaling pathway initiated by 1,25(OH)2D3 reveals that inhibition of the cyclic 3′,5′ -adenosine monophosphate-dependent kinase, protein kinase A, or inhibition of the cyclic 3′,5′-guanine monophosphate-dependent kinase, protein kinase G, also prevents the charge shift conversion. ^ Isolation of OPN from rat femurs and tibiae provides evidence for the existence of these two OPN charge forms in vivo, evidenced by differential migration on isoelectric focusing gels and sodium dodecyl sulfate-polyacrylamide gels. Peptide sequencing of rat long bone fractions revealed the presence of a presumed dentin specific protein, dentin matrix protein-1 (DMP-1). Western blot analysis confirmed the existence of DMP-1 in these fractions. ^ Using the OPN charge forms in functional assays, it was determined that the charge forms have differential roles in both cell surface and mineralization functions. In cell attachment assays and Ca2+ influx assays using PC-3 prostate cancer cells, the pI 5.1 charge form of OPN was found to permit binding and increase intracellular Ca2+ concentrations of PC-3 cells. The increase in intracellular Ca2+ concentration was found to be integrin αvβ3-dependent. In mineralization assays, the pI 4.6 charge form of OPN promoted hydroxyapatite formation, while the pI 5.1 charge form had improved Ca2+ binding ability. ^ In conclusion, these findings suggest that 1,25(OH) 2D3 regulates OPN not only at the transcriptional level, but also plays a role in determination of the OPN phosphorylation state. The latter involves a short term (less than three hours) treatment and is associated with membrane-initiated Ca2+ influx. Functional assays utilizing the two OPN charge forms reveal the dependence of OPN post-translational state on its function. ^

Clinical outcomes and molecular genotyping of Staphylococcus aureus isolated from milk samples of dairy primiparous Mediterranean buffaloes (Bubalus bubalis)

Staphylococcus aureus is one of the most important pathogens causing mastitis in dairy cows and in Mediterranean buffaloes. Genotype B (GTB) is contagious in dairy cows and may occur in up to 87% of cows of a dairy herd. It was the aim of this study to evaluate genotypes present, clinical outcomes, and prevalence of Staph. aureus in milk samples of primiparous Mediterranean dairy buffaloes. Two hundred composite milk samples originating from 40 primiparous buffaloes were collected from May to June 2012, at d 10, 30, 60, 90, and 150 d in milk (DIM) to perform somatic cell counts and bacteriological cultures. Daily milk yields were recorded. Before parturition until 40 to 50 DIM, all primiparous animals were housed separated from the pluriparous animals. Milking was performed in the same milking parlor, but the primiparous animals were milked first. After 50 DIM, the primiparous were mixed with the pluriparous animals, including the milking procedure. Individual quarter samples were collected from each animal, and aliquots of 1 mL were mixed and used for molecular identification and genotyping of Staph. aureus. The identification of Staph. aureus was performed verifying the presence of nuc gene by nuc gene PCR. All the nuc-positive isolates were subjected to genotype analysis by means of PCR amplification of the 16S-23S rRNA intergenic spacer region and analyzed by a miniaturized electrophoresis system. Of all 200 composite samples, 41 (20.5%) were positive for Staph. aureus, and no genotype other than GTB was identified. The prevalence of samples positive for Staph. aureus was 0% at 10 DIM and increased to a maximum of 22/40 (55%) at 90 DIM. During the period of interest, 14 buffaloes tested positive for Staph. aureus once, 6 were positive twice, and 5 were positive 3 times, whereas 15 animals were negative at every sampling. At 90 and 150 DIM, 7 (17.5%) and 3 buffaloes (7.5%), respectively, showed clinical mastitis (CM), and only 1 (2.5%) showed CM at both samplings. At 60, 90, and 150 DIM, 1 buffalo was found with subclinical mastitis at each sampling. At 30, 60, 90, and 150 DIM, 2.5 (1/40), 22.5 (9/40), 35 (14/40), and 10% (4/40) were considered affected by intramammary infection, respectively. Buffaloes with CM caused by Staph. aureus had statistically significantly higher mean somatic cell count values (6.06 ± 0.29, Log10 cells/mL ± standard deviation) and statistically significantly lower mean daily milk yields (7.15 ± 1.49, liters/animal per day) than healthy animals (4.69 ± 0.23 and 13.87 ± 2.64, respectively), buffaloes with IMI (4.82 ± 0.23 and 11.16 ± 1.80, respectively), or with subclinical mastitis (5.47 ± 0.10 and 10.33 ± 0.68, respectively). Based on our knowledge, this is the first time that Staph. aureus GTB has been identified in milk samples of dairy Mediterranean buffaloes.

Elucidation of Nucleic Acid–Drug Interactions by Tandem Mass Spectrometry

In continuation of the long tradition of mass spectrometric research at the University of Bern, our group focuses on the characterization of nucleic acids as therapeutic agents and as drug targets. This article provides a short overview of our recent work on platinated single-stranded and higher-order nucleic acids. Nearly three decades ago the development of soft ionization techniques opened a whole new chapter in the mass spectrometric analysis of not only nucleic acids themselves, but also their interactions with potential drug candidates. In contrast to modern next generation sequencing approaches, though, the goal of the tandem mass spectrometric investigation of nucleic acids is by no means the complete sequencing of genetic DNA, but rather the characterization of short therapeutic and regulatory oligonucleotides and the elucidation of nucleic acid–drug interactions. The influence of cisplatin binding on the gas-phase dissociation of nucleic acids was studied by the means of electrospray ionization tandem mass spectrometry. Experiments on native and modified DNA and RNA oligomers confirmed guanine base pairs as the preferred platination site and laid the basis for the formulation of a gas-phase fragmentation mechanism of platinated oligonucleotides. The study was extended to double stranded DNA and DNA quadruplexes. While duplexes are believed to be the main target of cisplatin in vivo, the recently discovered DNA quadruplexes constitute another promising target for anti-tumor drugs owing to their regulatory functions in the cell cycle.

Highlights in Analytical Chemistry in Switzerland: OMA & OPA - A Software Tool for Mass Spectrometric Sequencing of Nucleic Acids

Oligonucleotides comprising unnatural building blocks, which interfere with the translation machinery, have gained increased attention for the treatment of gene-related diseases (e.g. antisense, RNAi). Due to structural modifications, synthetic oligonucleotides exhibit increased biostability and bioavailability upon administration. Consequently, classical enzyme-based sequencing methods are not applicable to their sequence elucidation and verification. Tandem mass spectrometry is the method of choice for performing such tasks, since gas-phase dissociation is not restricted to natural nucleic acids. However, tandem mass spectrometric analysis can generate product ion spectra of tremendous complexity, as the number of possible fragments grows rapidly with increasing sequence length. The fact that structural modifications affect the dissociation pathways greatly increases the variety of analytically valuable fragment ions. The gas-phase dissociation of oligonucleotides is characterized by the cleavage of one of the four bonds along the phosphodiester chain, by the accompanying loss of nucleases, and by the generation of internal fragments due to secondary backbone cleavage. For example, an 18-mer oligonucleotide yields a total number of 272’920 theoretical fragment ions. In contrast to the processing of peptide product ion spectra, which nowadays is highly automated, there is a lack of tools assisting the interpretation of oligonucleotide data. The existing web-based and stand-alone software applications are primarily designed for the sequence analysis of natural nucleic acids, but do not account for chemical modifications and adducts. Consequently, we developed a software to support the interpretation of mass spectrometric data of natural and modified nucleic acids and their adducts with chemotherapeutic agents.

The brown coat colour of Coppernecked goats is associated with a non-synonymous variant at the TYRP1 locus on chromosome 8

The recent development of a goat SNP genotyping microarray enables genome-wide association studies in this important livestock species. We investigated the genetic basis of the black and brown coat colour in Valais Blacknecked and Coppernecked goats. A genome-wide association analysis using goat SNP50 BeadChip genotypes of 22 cases and 23 controls allowed us to map the locus for the brown coat colour to goat chromosome 8. The TYRP1 gene is located within the associated chromosomal region, and TYRP1 variants cause similar coat colour phenotypes in different species. We thus considered TYRP1 as a strong positional and functional candidate. We resequenced the caprine TYRP1 gene by Sanger and Illumina sequencing and identified two non-synonymous variants, p.Ile478Thr and p.Gly496Asp, that might have a functional impact on the TYRP1 protein. However, based on the obtained pedigree and genotype data, the brown coat colour in these goats is not due to a single recessive loss-of-function allele. Surprisingly, the genotype distribution and the pedigree data suggest that the (496) Asp allele might possibly act in a dominant manner. The (496) Asp allele was present in 77 of 81 investigated Coppernecked goats and did not occur in black goats. This strongly suggests heterogeneity underlying the brown coat colour in Coppernecked goats. Functional experiments or targeted matings will be required to verify the unexpected preliminary findings.

Independent polled mutations leading to complex gene expression differences in cattle

The molecular regulation of horn growth in ruminants is still poorly understood. To investigate this process, we collected 1019 hornless (polled) animals from different cattle breeds. High-density SNP genotyping confirmed the presence of two different polled associated haplotypes in Simmental and Holstein cattle co-localized on BTA 1. We refined the critical region of the Simmental polled mutation to 212 kb and identified an overlapping region of 932 kb containing the Holstein polled mutation. Subsequently, whole genome sequencing of polled Simmental and Holstein cows was used to determine polled associated genomic variants. By genotyping larger cohorts of animals with known horn status we found a single perfectly associated insertion/deletion variant in Simmental and other beef cattle confirming the recently published possible Celtic polled mutation. We identified a total of 182 sequence variants as candidate mutations for polledness in Holstein cattle, including an 80 kb genomic duplication and three SNPs reported before. For the first time we showed that hornless cattle with scurs are obligate heterozygous for one of the polled mutations. This is in contrast to published complex inheritance models for the bovine scurs phenotype. Studying differential expression of the annotated genes and loci within the mapped region on BTA 1 revealed a locus (LOC100848215), known in cow and buffalo only, which is higher expressed in fetal tissue of wildtype horn buds compared to tissue of polled fetuses. This implicates that the presence of this long noncoding RNA is a prerequisite for horn bud formation. In addition, both transcripts associated with polledness in goat and sheep (FOXL2 and RXFP2), show an overexpression in horn buds confirming their importance during horn development in cattle.