Roan coat color in livestock.

Since domestication, a wide variety of phenotypes including coat color variation has developed in livestock. This variation is mostly based on selective breeding. During the beginning of selective breeding, potential negative consequences did not become immediately evident due to low frequencies of homozygous animals and have been occasionally neglected. However, numerous studies of coat color genetics have been carried out over more than a century and, meanwhile, pleiotropic effects for several coat color genes, including disorders of even lethal impact, were described. Similar coat color phenotypes can often be found across species, caused either by conserved genes or by different genes. Even in the same species, more than one gene could cause the same or similar coat color phenotype. The roan coat color in livestock species is characterized by a mixture of white and colored hair in cattle, pig, sheep, goat, alpaca, and horse. So far, the genetic background of this phenotype is not fully understood, but KIT and its ligand KITLG (MGF) are major candidate genes in livestock species. For some of these species, pleiotropic effects such as subfertility in homozygous roan cattle or homozygous embryonic lethality in certain horse breeds have been described. This review aims to point out the similarities and differences of the roan phenotype across the following livestock species: cattle, pig, sheep, goat, alpaca, and horse; and provides the current state of knowledge on genetic background and pleiotropic effects.

A 50-kb deletion disrupting the RSPO2 gene is associated with tetradysmelia in Holstein Friesian cattle.

BACKGROUND: Tetradysmelia is a rare genetic disorder that is characterized by an extremely severe reduction of all limb parts distal of the scapula and pelvic girdle. We studied a Holstein Friesian backcross family with 24 offspring, among which six calves displayed autosomal recessive tetradysmelia. In order to identify the genetic basis of the disorder, we genotyped three affected calves, five dams and nine unaffected siblings using a Bovine Illumina 50 k BeadChip and sequenced the whole genome of the sire. RESULTS: Pathological examination of four tetradysmelia cases revealed a uniform and severe dysmelia of all limbs. Applying a homozygosity mapping approach, we identified a homozygous region of 10.54 Mb on chromosome 14 (Bos taurus BTA14). Only calves that were diagnosed with tetradysmelia shared a distinct homozygous haplotype for this region. We sequenced the whole genome of the cases' sire and searched for heterozygous single nucleotide polymorphisms (SNPs) and small variants on BTA14 that were uniquely present in the sire and absent from 3102 control whole-genome sequences of the 1000 Bull Genomes Project, but none were identified in the 10.54-Mb candidate region on BTA14. Therefore, we subsequently performed a more comprehensive analysis by also considering structural variants and detected a 50-kb deletion in the targeted chromosomal region that was in the heterozygous state in the cases' sire. Using PCR, we confirmed that this detected deletion segregated perfectly within the family with tetradysmelia. The deletion spanned three exons of the bovine R-spondin 2 (RSPO2) gene, which encode three domains of the respective protein. R-spondin 2 is a secreted ligand of leucine-rich repeats containing G protein-coupled receptors that enhance Wnt signalling and is involved in a broad range of developmental processes during embryogenesis. CONCLUSIONS: We identified a 50-kb deletion on BTA14 that disrupts the coding sequence of the RSPO2 gene and is associated with bovine tetradysmelia. To our knowledge, this is the first reported candidate causal mutation for tetradysmelia in a large animal model. Since signalling pathways involved in limb development are conserved across species, the observed inherited defect may serve as a model to further elucidate fundamental pathways of limb development.

Optic nerve hypoplasia in miniature poodle dogs: A preliminary genetic and candidate gene association study.

OBJECTIVE: To conduct a genetic and candidate gene association study with samples from phenotype-ascertained dogs to identify putative disease-associated gene/mutation for optic nerve hypoplasia (ONH) in the miniature poodle. ANIMALS STUDIED: A total of 43 miniature poodles from the United States and Europe, nine affected bilaterally with ONH, were included in the study. Pedigree information was recorded. PROCEDURES: A pedigree including all animals studied was assembled. Twenty-one genes typically expressed in ganglion cells or that are associated with ocular malformations and have a critical function in eye and neural retina development were selected. Exons and exon-intron boundaries of eight genes were sequenced in four ONH cases and four controls. Furthermore, cases and controls were genotyped with the Illumina CanineHD BeadChip to obtain genotypes for 13 additional candidate genes for haplotype association. RESULTS: The assembled pedigree connected all ONH-affected dogs to a possible common founder. Identified variants and haplotypes of the tested candidate genes did not segregate with the phenotype using Identity by Descent approach assuming autosomal recessive inheritance with variable but yet unknown penetrance. CONCLUSIONS: Pedigree analysis did not reveal the inheritance pattern. There is no evidence of association of the evaluated candidate genes with ONH; therefore, the screened candidate genes can provisionally be ruled out as causally associated with the disease.

A missense change in the ATG4D gene links aberrant autophagy to a neurodegenerative vacuolar storage disease.

Inherited neurodegenerative disorders are debilitating diseases that occur across different species. We have performed clinical, pathological and genetic studies to characterize a novel canine neurodegenerative disease present in the Lagotto Romagnolo dog breed. Affected dogs suffer from progressive cerebellar ataxia, sometimes accompanied by episodic nystagmus and behavioral changes. Histological examination revealed unique pathological changes, including profound neuronal cytoplasmic vacuolization in the nervous system, as well as spheroid formation and cytoplasmic aggregation of vacuoles in secretory epithelial tissues and mesenchymal cells. Genetic analyses uncovered a missense change, c.1288G>A; p.A430T, in the autophagy-related ATG4D gene on canine chromosome 20 with a highly significant disease association (p = 3.8 x 10-136) in a cohort of more than 2300 Lagotto Romagnolo dogs. ATG4D encodes a poorly characterized cysteine protease belonging to the macroautophagy pathway. Accordingly, our histological analyses indicated altered autophagic flux in affected tissues. The knockdown of the zebrafish homologue atg4da resulted in a widespread developmental disturbance and neurodegeneration in the central nervous system. Our study describes a previously unknown canine neurological disease with particular pathological features and implicates the ATG4D protein as an important autophagy mediator in neuronal homeostasis. The canine phenotype serves as a model to delineate the disease-causing pathological mechanism(s) and ATG4D function, and can also be used to explore treatment options. Furthermore, our results reveal a novel candidate gene for human neurodegeneration and enable the development of a genetic test for veterinary diagnostic and breeding purposes.

A GJA9 frameshift variant is associated with polyneuropathy in Leonberger dogs.

BACKGROUND: Many inherited polyneuropathies (PN) observed in dogs have clinical similarities to the genetically heterogeneous group of Charcot-Marie-Tooth (CMT) peripheral neuropathies in humans. The canine disorders collectively show a variable expression of progressive clinical signs and ages of onset, and different breed prevalences. Previously in the Leonberger breed, a variant highly associated with a juvenile-onset PN was identified in the canine orthologue of a CMT-associated gene. As this deletion in ARHGEF10 (termed LPN1) does not explain all cases, PN in this breed may encompass variants in several genes with similar clinical and histopathological features. RESULTS: A genome-wide comparison of 173 k SNP genotypes of 176 cases, excluding dogs homozygous for the ARHGEF10 variant, and 138 controls, was carried out to detect further PN-associated variants. A single suggestive significant association signal on CFA15 was found. The genome of a PN-affected Leonberger homozygous for the associated haplotype was sequenced and variants in the 7.7 Mb sized critical interval were identified. These variants were filtered against a database of variants observed in 202 genomes of various dog breeds and 3 wolves, and 6 private variants in protein-coding genes, all in complete linkage disequilibrium, plus 92 non-coding variants were revealed. Five of the coding variants were predicted to have low or moderate effect on the encoded protein, whereas a 2 bp deletion in GJA9 results in a frameshift of high impact. GJA9 encodes connexin 59, a connexin gap junction family protein, and belongs to a group of CMT-associated genes that have emerged as important components of peripheral myelinated nerve fibers. The association between the GJA9 variant and PN was confirmed in an independent cohort of 296 cases and 312 controls. Population studies showed a dominant mode of inheritance, an average age of onset of approximately 6 years, and incomplete penetrance. CONCLUSIONS: This GJA9 variant represents a highly probable candidate variant for another form of PN in Leonberger dogs, which we have designated LPN2, and a new candidate gene for CMT disease. To date, approximately every third PN-diagnosed Leonberger dog can be explained by the ARHGEF10 or GJA9 variants, and we assume that additional genetic heterogeneity in this condition exists in the breed.

An ARHGEF10 deletion is highly associated with a juvenile-onset inherited polyneuropathy in Leonberger and Saint Bernard dogs.

An inherited polyneuropathy (PN) observed in Leonberger dogs has clinical similarities to a genetically heterogeneous group of peripheral neuropathies termed Charcot-Marie-Tooth (CMT) disease in humans. The Leonberger disorder is a severe, juvenile-onset, chronic, progressive, and mixed PN, characterized by exercise intolerance, gait abnormalities and muscle atrophy of the pelvic limbs, as well as inspiratory stridor and dyspnea. We mapped a PN locus in Leonbergers to a 250 kb region on canine chromosome 16 (Praw = 1.16×10-10, Pgenome, corrected = 0.006) utilizing a high-density SNP array. Within this interval is the ARHGEF10 gene, a member of the rho family of GTPases known to be involved in neuronal growth and axonal migration, and implicated in human hypomyelination. ARHGEF10 sequencing identified a 10 bp deletion in affected dogs that removes four nucleotides from the 3'-end of exon 17 and six nucleotides from the 5'-end of intron 17 (c.1955_1958+6delCACGGTGAGC). This eliminates the 3'-splice junction of exon 17, creates an alternate splice site immediately downstream in which the processed mRNA contains a frame shift, and generates a premature stop codon predicted to truncate approximately 50% of the protein. Homozygosity for the deletion was highly associated with the severe juvenile-onset PN phenotype in both Leonberger and Saint Bernard dogs. The overall clinical picture of PN in these breeds, and the effects of sex and heterozygosity of the ARHGEF10 deletion, are less clear due to the likely presence of other forms of PN with variable ages of onset and severity of clinical signs. This is the first documented severe polyneuropathy associated with a mutation in ARHGEF10 in any species.

A mutation in the FAM83G gene in dogs with hereditary footpad hyperkeratosis (HFH).

Hereditary footpad hyperkeratosis (HFH) represents a palmoplantar hyperkeratosis, which is inherited as a monogenic autosomal recessive trait in several dog breeds, such as e.g. Kromfohrländer and Irish Terriers. We performed genome-wide association studies (GWAS) in both breeds. In Kromfohrländer we obtained a single strong association signal on chromosome 5 (p(raw) = 1.0×10(-13)) using 13 HFH cases and 29 controls. The association signal replicated in an independent cohort of Irish Terriers with 10 cases and 21 controls (p(raw) = 6.9×10(-10)). The analysis of shared haplotypes among the combined Kromfohrländer and Irish Terrier cases defined a critical interval of 611 kb with 13 predicted genes. We re-sequenced the genome of one affected Kromfohrländer at 23.5× coverage. The comparison of the sequence data with 46 genomes of non-affected dogs from other breeds revealed a single private non-synonymous variant in the critical interval with respect to the reference genome assembly. The variant is a missense variant (c.155G>C) in the FAM83G gene encoding a protein with largely unknown function. It is predicted to change an evolutionary conserved arginine into a proline residue (p.R52P). We genotyped this variant in a larger cohort of dogs and found perfect association with the HFH phenotype. We further studied the clinical and histopathological alterations in the epidermis in vivo. Affected dogs show a moderate to severe orthokeratotic hyperplasia of the palmoplantar epidermis. Thus, our data provide the first evidence that FAM83G has an essential role for maintaining the integrity of the palmoplantar epidermis.

Molecular studies of phenotype variation in canine RPGR-XLPRA1.

PURPOSE: Canine X-linked progressive retinal atrophy 1 (XLPRA1) caused by a mutation in retinitis pigmentosa (RP) GTPase regulator (RPGR) exon ORF15 showed significant variability in disease onset in a colony of dogs that all inherited the same mutant X chromosome. Defective protein trafficking has been detected in XLPRA1 before any discernible degeneration of the photoreceptors. We hypothesized that the severity of the photoreceptor degeneration in affected dogs may be associated with defects in genes involved in ciliary trafficking. To this end, we examined six genes as potential disease modifiers. We also examined the expression levels of 24 genes involved in ciliary trafficking (seven), visual pathway (five), neuronal maintenance genes (six), and cellular stress response (six) to evaluate their possible involvement in early stages of the disease. METHODS: Samples from a pedigree derived from a single XLPRA1-affected male dog outcrossed to unrelated healthy mix-bred or purebred females were used for immunohistochemistry (IHC), western blot, mutational and haplotype analysis, and gene expression (GE). Cell-specific markers were used to examine retinal remodeling in the disease. Single nucleotide polymorphisms (SNPs) spanning the entire RPGR interacting and protein trafficking genes (RAB8A, RPGRIP1L, CEP290, CC2D2A, DFNB31, and RAB11B) were genotyped in the pedigree. Quantitative real-time PCR (qRT-PCR) was used to examine the expression of a total of 24 genes, including the six genes listed. RESULTS: Examination of cryosections from XLPRA1-affected animals of similar age (3-4 years) with different disease severity phenotype revealed mislocalization of opsins and upregulation of the Müller cell gliosis marker GFAP. Four to ten haplotypes per gene were identified in RAB8A, RPGRIP1L, CEP290, CC2D2A, DFNB31, and RAB11B for further assessment as potential genetic modifiers of XLPRA1. No correlation was found between the haplotypes and disease severity. During mutational analysis, several new variants, including a single intronic mutation in RAB8A and three mutations in exon 3 of DFNB31 were described (c.970G>A (V324I), c.978T>C (G326=), and c.985G>A (A329T)). Expression analysis of stress response genes in 16-week-old predisease XLPRA1 retinas revealed upregulation of GFAP but not HSPA5, DDIT3, HSPA4, HSP90B1, or HIF1A. Western blot analysis confirmed GFAP upregulation. In the same predisease group, no significant differences were found in the expression of 18 selected genes (RHO, OPN1LW, OPN1MW, RLBP1, RPGRORF15, RAB8A, RPGRIP1L, CEP290, CC2D2A, DFNB31, RAB11B, CRX, RCVRN, PVALB, CALB1, FGFR1, NTRK2, and NTRK3) involved in neuronal function. CONCLUSIONS: Lack of association between haplotypes of RAB8A, RPGRIP1L, CEP290, CC2D2A, DFNB31, and RAB11B and the disease phenotype suggests that these genes are not genetic modifiers of XLPRA1. Upregulation of GFAP, an established indicator of the Müller cell gliosis, manifests as an important early feature of the disease.

Molecular Structural Analysis of Porcine CMAH-Native Ligand Complex and High Throughput Virtual Screening to Identify Novel Inhibitors.

Porcine meat is the most consumed red meat worldwide. Pigs are also vital tools in biological and medical research. However, xenoreactivity between porcine's N-glycolylneuraminic acid (Neu5Gc) and human anti-Neu5Gc antibodies poses a significant challenge. On the one hand, dietary Neu5Gc intake has been connected to particular human disorders. On the other hand, some pathogens connected to pig diseases have a preference for Neu5Gc. The Cytidine monophospho-N-acetylneuraminic acid hydroxylase (CMAH) catalyses the conversion of N-acetylneuraminic acid (Neu5Ac) to Neu5Gc. In this study, we predicted the tertiary structure of CMAH, performed molecular docking, and analysed the protein-native ligand complex. We performed a virtual screening from a drug library of 5M compounds and selected the two top inhibitors with Vina scores of -9.9 kcal/mol for inhibitor 1 and -9.4 kcal/mol for inhibitor 2. We further analysed their pharmacokinetic and pharmacophoric properties. We conducted stability analyses of the complexes with molecular dynamic simulations of 200 ns and binding free energy calculations. The overall analyses revealed the inhibitors' stable binding, which was further validated by the MMGBSA studies. In conclusion, this result may pave the way for future studies to determine how to inhibit CMAH activities. Further in vitro studies can provide in-depth insight into these compounds' therapeutic potential.

An In Silico Functional Analysis of Non-Synonymous Single-Nucleotide Polymorphisms of Bovine CMAH Gene and Potential Implication in Pathogenesis.

The sugar molecule N-glycolylneuraminic acid (Neu5Gc) is one of the most common sialic acids discovered in mammals. Cytidine monophospho-N-acetylneuraminic acid hydroxylase (CMAH) catalyses the conversion of N-acetylneuraminic acid (Neu5Ac) to Neu5Gc, and it is encoded by the CMAH gene. On the one hand, food metabolic incorporation of Neu5Gc has been linked to specific human diseases. On the other hand, Neu5Gc has been shown to be highly preferred by some pathogens linked to certain bovine diseases. We used various computational techniques to perform an in silico functional analysis of five non-synonymous single-nucleotide polymorphisms (nsSNPs) of the bovine CMAH (bCMAH) gene identified from the 1000 Bull Genomes sequence data. The c.1271C>T (P424L) nsSNP was predicted to be pathogenic based on the consensus result from different computational tools. The nsSNP was also predicted to be critical based on sequence conservation, stability, and post-translational modification site analysis. According to the molecular dynamic simulation and stability analysis, all variations promoted stability of the bCMAH protein, but mutation A210S significantly promoted CMAH stability. In conclusion, c.1271C>T (P424L) is expected to be the most harmful nsSNP among the five detected nsSNPs based on the overall studies. This research could pave the way for more research associating pathogenic nsSNPs in the bCMAH gene with diseases.

Improving the annotation of the cattle genome by annotating transcription start sites in a diverse set of tissues and populations using Cap Analysis Gene Expression sequencing.

Understanding the genomic control of tissue-specific gene expression and regulation can help to inform the application of genomic technologies in farm animal breeding programs. The fine mapping of promoters [transcription start sites (TSS)] and enhancers (divergent amplifying segments of the genome local to TSS) in different populations of cattle across a wide diversity of tissues provides information to locate and understand the genomic drivers of breed- and tissue-specific characteristics. To this aim, we used Cap Analysis Gene Expression (CAGE) sequencing, of 24 different tissues from 3 populations of cattle, to define TSS and their coexpressed short-range enhancers (<1 kb) in the ARS-UCD1.2_Btau5.0.1Y reference genome (1000bulls run9) and analyzed tissue and population specificity of expressed promoters. We identified 51,295 TSS and 2,328 TSS-Enhancer regions shared across the 3 populations (dairy, beef-dairy cross, and Canadian Kinsella composite cattle from 2 individuals, 1 of each sex, per population). Cross-species comparative analysis of CAGE data from 7 other species, including sheep, revealed a set of TSS and TSS-Enhancers that were specific to cattle. The CAGE data set will be combined with other transcriptomic information for the same tissues to create a new high-resolution map of transcript diversity across tissues and populations in cattle for the BovReg project. Here we provide the CAGE data set and annotation tracks for TSS and TSS-Enhancers in the cattle genome. This new annotation information will improve our understanding of the drivers of gene expression and regulation in cattle and help to inform the application of genomic technologies in breeding programs.

A missense mutation in the SERPINH1 gene in Dachshunds with osteogenesis imperfecta.

Osteogenesis imperfecta (OI) is a hereditary disease occurring in humans and dogs. It is characterized by extremely fragile bones and teeth. Most human and some canine OI cases are caused by mutations in the COL1A1 and COL1A2 genes encoding the subunits of collagen I. Recently, mutations in the CRTAP and LEPRE1 genes were found to cause some rare forms of human OI. Many OI cases exist where the causative mutation has not yet been found. We investigated Dachshunds with an autosomal recessive form of OI. Genotyping only five affected dogs on the 50 k canine SNP chip allowed us to localize the causative mutation to a 5.82 Mb interval on chromosome 21 by homozygosity mapping. Haplotype analysis of five additional carriers narrowed the interval further down to 4.74 Mb. The SERPINH1 gene is located within this interval and encodes an essential chaperone involved in the correct folding of the collagen triple helix. Therefore, we considered SERPINH1 a positional and functional candidate gene and performed mutation analysis in affected and control Dachshunds. A missense mutation (c.977C>T, p.L326P) located in an evolutionary conserved domain was perfectly associated with the OI phenotype. We thus have identified a candidate causative mutation for OI in Dachshunds and identified a fifth OI gene.

An Overview of the Importance and Value of Porcine Species in Sialic Acid Research.

Humans frequently interact with pigs, whose meat is also one of the primary sources of animal protein. They are one of the main species at the center of sialic acid (Sia) research. Sias are sugars at terminals of glycoconjugates, are expressed at the cell surfaces of mammals, and are important in cellular interactions. N-glycolylneuraminic acid (Neu5Gc) and N-acetylneuraminic acid (Neu5Ac) are notable Sias in mammals. Cytidine monophospho-N-acetylneuraminic acid hydroxylase (CMAH) encodes the CMAH enzyme that biosynthesizes Neu5Gc. Although humans cannot endogenously synthesize Neu5Gc due to the inactivation of this gene by a mutation, Neu5Gc can be metabolically incorporated into human tissues from red meat consumption. Interactions between Neu5Gc and human anti-Neu5Gc antibodies have been associated with certain diseases and disorders. In this review, we summarized the sialic acid metabolic pathway, its regulation and link to viral infections, as well as the importance of the pig as a model organism in Sia research, making it a possible source of Neu5Gc antigens affecting human health. Future research in solving the structures of crucial enzymes involved in Sia metabolism, as well as their regulation and interactions with other enzymes, especially CMAH, could help to understand their function and reduce the amount of Neu5Gc.

Mitochondrial DNA Variation Contributes to the Aptitude for Dressage and Show Jumping Ability in the Holstein Horse Breed.

Maternal lineages are considered an important factor in breeding. Mitochondrial DNA (mtDNA) is maternally inherited and plays an important role in energy metabolism. It has already been associated with energy consumption and performances, e.g., stamina in humans and racehorses. For now, corresponding studies are lacking for sport performance of warmblood breeds. MtDNA sequences were available for 271 Holstein mares from 75 maternal lineages. As all mares within a lineage showed identical haplotypes regarding the non-synonymous variants, we expanded our data set by also including non-sequenced mares and assigning them to the lineage-specific haplotype. This sample consisting of 6334 to 16,447 mares was used to perform mitochondrial association analyses using breeding values (EBVs) estimated on behalf of the Fédération Équestre Nationale (FN) and on behalf of the Holstein Breeding Association (HOL). The association analyses revealed 20 mitochondrial SNPs (mtSNPs) significantly associated with FN-EBVs and partly overlapping 20 mtSNPs associated with HOL-EBVs. The results indicated that mtDNA contributes to performance differences between maternal lineages. Certain mitochondrial haplogroups were associated with special talents for dressage or show jumping. The findings encourage to set up innovative genetic evaluation models that also consider information on maternal lineages.

Allele-biased expression of the bovine APOB gene associated with the cholesterol deficiency defect suggests cis-regulatory enhancer effects of the LTR retrotransposon insertion.

The insertion of an endogenous retroviral long terminal repeat (LTR) sequence into the bovine apolipoprotein B (APOB) gene is causal to the inherited genetic defect cholesterol deficiency (CD) observed in neonatal and young calves. Affected calves suffer from developmental abnormalities, symptoms of incurable diarrhoea and often die within weeks to a few months after birth. Neither the detailed effects of the LTR insertion on APOB expression profile nor the specific mode of inheritance nor detailed phenotypic consequences of the mutation are undisputed. In our study, we analysed German Holstein dairy heifers at the peak of hepatic metabolic load and exposed to an additional pathogen challenge for clinical, metabolic and hepatic transcriptome differences between wild type (CDF) and heterozygote carriers of the mutation (CDC). Our data revealed that a divergent allele-biased expression pattern of the APOB gene in heterozygous CDC animals leads to a tenfold higher expression of exons upstream and a decreased expression of exons downstream of the LTR insertion compared to expression levels of CDF animals. This expression pattern could be a result of enhancer activity induced by the LTR insertion, in addition to a previously reported artificial polyadenylation signal. Thus, our data support a regulatory potential of mobile element insertions. With regard to the phenotype generated by the LTR insertion, heterozygote CDC carriers display significantly differential hepatic expression of genes involved in cholesterol biosynthesis and lipid metabolism. Phenotypically, CDC carriers show a significantly affected lipomobilization compared to wild type animals. These results reject a completely recessive mode of inheritance for the CD defect, which should be considered for selection decisions in the affected population. Exemplarily, our results illustrate the regulatory impact of mobile element insertions not only on specific host target gene expression but also on global transcriptome profiles with subsequent biological, functional and phenotypic consequences in a natural in-vivo model of a non-model mammalian organism.

Single-cell RNA sequencing of freshly isolated bovine milk cells and cultured primary mammary epithelial cells.

Bovine mammary function at molecular level is often studied using mammary tissue or primary bovine mammary epithelial cells (pbMECs). However, bulk tissue and primary cells are heterogeneous with respect to cell populations, adding further transcriptional variation in addition to genetic background. Thus, understanding of the variation in gene expression profiles of cell populations and their effect on function are limited. To investigate the mononuclear cell composition in bovine milk, we analyzed a single-cell suspension from a milk sample. Additionally, we harvested cultured pbMECs to characterize gene expression in a homogeneous cell population. Using the Drop-seq technology, we generated single-cell RNA datasets of somatic milk cells and pbMECs. The final datasets after quality control filtering contained 7,119 and 10,549 cells, respectively. The pbMECs formed 14 indefinite clusters displaying intrapopulation heterogeneity, whereas the milk cells formed 14 more distinct clusters. Our datasets constitute a molecular cell atlas that provides a basis for future studies of milk cell composition and gene expression, and could serve as reference datasets for milk cell analysis.

Exploring the Origin and Relatedness of Maternal Lineages Through Analysis of Mitochondrial DNA in the Holstein Horse.

Maternal lineages are important for the breeding decision in the Holstein horse breed. To investigate the genetic diversity of the maternal lineages and the relationships between founder mares, the maternal inherited mitochondrial genome (except the repetitive part of the non-coding region) of 271 mares representing 75 lineages was sequenced. The sequencing predominantly revealed complete homology in the nucleotide sequences between mares from one lineage with exceptions in 13 lineages, where differences in one to three positions are probably caused by de novo mutations or alternate fixation of heteroplasmy. We found 78 distinct haplotypes that have not yet been described in other breeds. Six of these occurred in two or three different lineages indicating a common ancestry. Haplotypes can be divided into eight clusters with all mares from one lineage belonging to the same cluster. Within a cluster, the average number of pairwise differences ranged from zero to 16.49 suggesting close maternal relationships between these mares. The results showed that the current breeding population originated from at least eight ancestral founder mares.

Coat Color Roan Shows Association with KIT Variants and No Evidence of Lethality in Icelandic Horses.

Roan (Rn) horses show a typical seasonal change of color. Their body is covered with colored and white hair. We performed a descriptive statistical analysis of breeding records of Icelandic horses to challenge the hypothesis of roan being lethal in utero under homozygous condition. The roan to non-roan ratio of foals from roan × roan matings revealed homozygous roan Icelandic horses to be viable. Even though roan is known to be inherited in a dominant mode and epistatic to other coat colors, the causative mutation is still unknown. Nevertheless, an association between roan phenotype and the KIT gene was shown for different horse breeds. In the present study, we identified KIT variants by Sanger sequencing, and show that KIT is also associated with roan in the Icelandic horse breed.

Candidate Genetic Modifiers for RPGR Retinal Degeneration.

Purpose: To define genetic variants associated with variable severity of X-linked progressive retinal atrophy 1 (XLPRA1) caused by a five-nucleotide deletion in canine RPGR exon ORF15. Methods: A genome-wide association study (GWAS) was performed in XLPRA1 phenotype informative pedigree. Whole genome sequencing (WGS) was used for mutational analysis of genes within the candidate genomic region. Retinas of normal and mutant dogs were used for gene expression, gene structure, and RNA duplex analyses. Results: GWAS followed by haplotype phasing identified an approximately 4.6 Mb candidate genomic interval on CFA31 containing seven protein-coding genes expressed in retina (ROBO1, ROBO2, RBM11, NRIP1, HSPA13, SAMSN1, and USP25). Furthermore, we identified and characterized two novel lncRNAs, ROBO1-AS and ROBO2-AS, that display overlapping gene organization with axon guidance pathway genes ROBO1 and ROBO2, respectively, producing sense-antisense gene pairs. Notably, ROBO1-AS and ROBO2-AS act in cis to form lncRNA/mRNA duplexes with ROBO1 and ROBO2, respectively, suggesting important roles for these lncRNAs in the ROBO regulatory network. A subsequent WGS identified candidate genes within the genomic region on CFA31 that might be implicated in modifying severity of XLPRA1. This approach led to discovery of genetic variants in ROBO1, ROBO1-AS, ROBO2-AS, and USP25 that are strongly associated with the XLPRA1 moderate phenotype. Conclusions: The study provides new insights into the genetic basis of phenotypic variation in severity of RPGRorf15-associated retinal degeneration. Our findings suggest an important role for ROBO pathways in disease progression further expanding on our previously reported changes of ROBO1 expression in XLPRA1 retinas.

CCDC66 frameshift variant associated with a new form of early-onset progressive retinal atrophy in Portuguese Water Dogs.

Aberrant photoreceptor function or morphogenesis leads to blinding retinal degenerative diseases, the majority of which have a genetic aetiology. A variant in PRCD previously identified in Portuguese Water Dogs (PWDs) underlies prcd (progressive rod-cone degeneration), an autosomal recessive progressive retinal atrophy (PRA) with a late onset at 3-6 years of age or older. Herein, we have identified a new form of early-onset PRA (EOPRA) in the same breed. Pedigree analysis suggested an autosomal recessive inheritance. Four PWD full-siblings affected with EOPRA diagnosed at 2-3 years of age were genotyped (173,661 SNPs) along with 2 unaffected siblings, 2 unaffected parents, and 15 unrelated control PWDs. GWAS, linkage analysis and homozygosity mapping defined a 26-Mb candidate region in canine chromosome 20. Whole-genome sequencing in one affected dog and its obligatory carrier parents identified a 1 bp insertion (CFA20:g.33,717,704_33,717,705insT (CanFam3.1); c.2262_c.2263insA) in CCDC66 predicted to cause a frameshift and truncation (p.Val747SerfsTer8). Screening of an extended PWD population confirmed perfect co-segregation of this genetic variant with the disease. Western blot analysis of COS-1 cells transfected with recombinant mutant CCDC66 expression constructs showed the mutant transcript translated into a truncated protein. Furthermore, in vitro studies suggest that the mutant CCDC66 is mislocalized to the nucleus relative to wild type CCDC66. CCDC66 variants have been associated with inherited retinal degenerations (RDs) including canine and murine ciliopathies. As genetic variants affecting the primary cilium can cause ciliopathies in which RD may be either the sole clinical manifestation or part of a syndrome, our findings further support a role for CCDC66 in retinal function and viability, potentially through its ciliary function.