Rapid Discovery of De Novo Deleterious Mutations in Cattle Enhances the Value of Livestock as Model Species

PAX3 haploinsufficiency in Maine Coon cats with dominant blue eyes and hearing loss resembling the human Waardenburg syndrome

This study investigated the dominant blue eyes (DBE) trait linked to hearing impairment and variable white spotting in Maine Coon cats. Fifty-eight animals descending from 2 different DBE lineages, the Dutch and the Topaz lines, were sampled. They comprised 48 cats from the Dutch bloodline, including 9 green-eyed and 31 blue-eyed cats, with some individuals exhibiting signs of deafness, and 8 stillborn kittens. Samples from the Topaz lineage included 10 blue-eyed animals. A brainstem auditory evoked response test revealed a reduced to absent response to auditory stimuli and absent physiological waveforms in all of the 8 examined DBE animals. We sequenced the genome of 2 affected cats from the Dutch line and searched for variants in 19 candidate genes for the human Waardenburg syndrome and pigmentary disorders. This search yielded 9 private protein-changing candidate variants in the genes PAX3, EDN3, KIT, OCA2, SLC24A5, HERC2, and TYRP1. The genotype-phenotype cosegregation was observed for the PAX3 variant within all animals from the Dutch lineage. The mutant allele was absent from 461 control genomes and 241 additionally genotyped green-eyed Maine Coons. We considered the PAX3 variant as the most plausible candidate-a heterozygous nonsense single base pair substitution in exon 6 of PAX3 (NC_051841.1:g.205,787,310G>A, XM_019838731.3:c.937C>T, XP_019694290.1:p.Gln313*), predicted to result in a premature stop codon. PAX3 variants cause auditory-pigmentary syndrome in humans, horses, and mice. Together with the comparative data from other species, our findings strongly suggest PAX3:c.937C>T (OMIA:001688-9685) as the most likely candidate variant for the DBE, deafness, and minimal white spotting in the Maine Coon Dutch line. Finally, we propose the designation of DBERE (Rociri Elvis Dominant Blue Eyes) allele in the domestic cat.

Naturally occurring genetic diseases caused by de novo variants in domestic animals

With the advent of next-generation sequencing, an increasing number of cases of de novo variants in domestic animals have been reported in scientific literature primarily associated with clinically severe phenotypes. The emergence of new variants at each generation is a crucial aspect in understanding the pathology of early-onset diseases in animals and can provide valuable insights into similar diseases in humans. With the aim of collecting deleterious de novo variants in domestic animals, we searched the scientific literature and compiled reports on 42 de novo variants in 31 genes in domestic animals. No clear disease-associated phenotype has been established in humans for three of these genes (NUMB, ANKRD28 and KCNG1). For the remaining 28 genes, a strong similarity between animal and human phenotypes was recognized from available information in OMIM and OMIA, revealing the importance of comparative studies and supporting the use of domestic animals as natural models for human diseases, in line with the One Health approach.

A recurrent de novo missense mutation in COL1A1 causes osteogenesis imperfecta type II and preterm delivery in Normande cattle

BACKGROUND

Nine male and eight female calves born to a Normande artificial insemination bull named "Ly" were referred to the French National Observatory of Bovine Abnormalities for multiple fractures, shortened gestation, and stillbirth or perinatal mortality.

RESULTS

Using Illumina BovineSNP50 array genotypes from affected calves and 84 half-sib controls, the associated locus was mapped to a 6.5-Mb interval on chromosome 19, assuming autosomal inheritance with germline mosaicism. Subsequent comparison of the whole-genome sequences of one case and 5116 control genomes, followed by genotyping in the affected pedigree, identified a de novo missense substitution within the NC1 domain of the COL1A1 gene (Chr19 g.36,473,965G > A; p.D1412N) as unique candidate variant. Interestingly, the affected residue was completely conserved among 243 vertebrate orthologs, and the same substitution in humans has been reported to cause type II osteogenesis imperfecta (OI), a connective tissue disorder that is characterized primarily by bone deformity and fragility. Moreover, three COL1A1 mutations have been described to cause the same syndrome in cattle. Necropsy, computed tomography, radiology, and histology confirmed the diagnosis of type II OI, further supporting the causality of this variant. In addition, a detailed analysis of gestation length and perinatal mortality in 1387 offspring of Ly and more than 160,000 progeny of 63 control bulls allowed us to statistically confirm in a large pedigree the association between type II OI and preterm delivery, which is probably due to premature rupture of fetal membranes and has been reported in several isolated cases of type II OI in humans and cattle. Finally, analysis of perinatal mortality rates and segregation distortion supported a low level of germ cell mosaicism in Ly, with an estimate of 4.5% to 7.7% of mutant sperm and thus 63 to 107 affected calves born. These numbers contrast with the 17 cases reported and raise concerns about the underreporting of congenital defects to heredo-surveillance platforms, even for textbook genetic syndromes.

CONCLUSIONS

In conclusion, we describe a large animal model for a recurrent substitution in COL1A1 that is responsible for type II OI in humans. More generally, this study highlights the utility of such datasets and large half-sib families available in livestock species to characterize sporadic genetic defects.

SNPs in microRNA seed region and impact of miR-375 in concurrent regulation of multiple lipid accumulation-related genes

Bovine intramuscular fat (IMF), commonly referred to as marbling, is regulated by lipid metabolism, which includes adipogenesis, lipogenesis, glycerolipid synthesis, and lipolysis. In recent years, breeding researchers have identified single nucleotide polymorphisms (SNPs) as useful marker-assisted selection tools for improving marbling scores in national breeding programs. These included causal SNPs that induce phenotypic variation. MicroRNAs (miRNAs) are small highly conserved non-coding RNA molecules that bind to multiple non-coding regions. They are involved in post-transcriptional regulation. Multiple miRNAs may regulate a given target. Previously, three SNPs in the GPAM 3' UTR and four miRNAs were identified through in silico assays. The aim of this study is to verify the binding ability of the four miRNAs to the SNPs within the 3'UTR of GPAM, and to identify the regulatory function of miR-375 in the expression of genes related to lipid metabolism in mammalian adipocytes. It was verified that the four miRNAs bind to the GPAM 3'UTR, and identified that the miR-375 sequence is highly conserved. Furthermore, it was founded that miR-375 upregulated the GPAM gene, C/EBPα, PPARγ and lipid metabolism-related genes and promoted lipid droplet accumulation in 3T3-L1 cells. In conclusion, these results suggest that miR-375 is a multifunctional regulator of multiple lipid metabolism-related genes and may aid in obesity research as a biomarker.

A novel imprinted locus on bovine chromosome 18 homologous with human chromosome 16q24.1

Genomic imprinting is an epigenetic regulation mechanism in mammals resulting in the parentally dependent monoallelic expression of genes. Imprinting disorders in humans are associated with several congenital syndromes and cancers and remain the focus of many medical studies. Cattle is a better model organism for investigating human embryo development than mice. Imprinted genes usually cluster on chromosomes and are regulated by different methylation regions (DMRs) located in imprinting control regions that control gene expression in cis. There is an imprinted locus on human chromosome 16q24.1 associated with congenital lethal developmental lung disease in newborns. However, genomic imprinting on bovine chromosome 18, which is homologous with human chromosome 16 has not been systematically studied. The aim of this study was to analyze the allelic expressions of eight genes (CDH13, ATP2C2, TLDC1, COTL1, CRISPLD2, ZDHHC7, KIAA0513, and GSE1) on bovine chromosome 18 and to search the DMRs associated gene allelic expression. Three transcript variants of the ZDHHC7 gene (X1, X2, and X5) showed maternal imprinting in bovine placentas. In addition, the monoallelic expression of X2 and X5 was tissue-specific. Five transcripts of the KIAA0513 gene showed tissue- and isoform-specific monoallelic expression. The CDH13, ATP2C2, and TLDC1 genes exhibited tissue-specific imprinting, however, COTL1, CRISLPLD2, and GSE1 escaped imprinting. Four DMRs, established after fertilization, were found in this region. Two DMRs were located between the ZDHHC7 and KIAA0513 genes, and two were in exon 1 of the CDH13 and ATP2C2 genes, respectively. The results from this study support future studies on the molecular mechanism to regulate the imprinting of candidate genes on bovine chromosome 18.

A missense mutation in the highly conserved TNF-like domain of Ectodysplasin A is the candidate causative variant for X-linked hypohidrotic ectodermal dysplasia in Limousin cattle: Clinical, histological, and molecular analyses

Ectodysplasin A related hypohidrotic ectodermal dysplasia (XLHED) is a well-studied fetal developmental disorder in mammals that mainly affects ectodermal structures. It has been identified in a variety of species, including mice, rats, dogs, cattle, and humans. Here, we report the clinical, histological, and molecular biological analyses of a case of XLHED in Limousin cattle. An affected Limousin calf showed pathognomonic signs of ectodermal dysplasia, i.e. sparse hair and characteristic dental aplasia. Histopathologic comparison of hairy and glabrous skin and computed tomography of the mandible confirmed the phenotypic diagnosis. In addition, a keratoconjunctivitis sicca was noted in one eye, which was also confirmed histopathologically. To identify the causative variant, we resequenced the bovine X-chromosomal ectodysplasin A gene (EDA) of the affected calf and compared the sequences to the bovine reference genome. A single missense variant (rs439722471) at position X:g.80411716T>C (ARS-UCD1.3) was identified. The variant resulted in an amino acid substitution from glutamic acid to glycine within the highly conserved TNF-like domain. To rule out the possibility that the variant was relatively common in the cattle population we genotyped 2,016 individuals including 40% Limousin cattle by fluorescence resonance energy transfer analysis. We also tested 5,116 multibreed samples from Run9 of the 1000 Bull Genomes Project for the said variant. The variant was not detected in any of the cattle tested, confirming the assumption that it was the causative variant. This is the first report of Ectodysplasin A related hypohidrotic ectodermal dysplasia in Limousin cattle and the description of a novel causal variant in cattle.

A Missense Mutation in the Collagen Triple Helix of EDA Is Associated with X-Linked Recessive Hypohidrotic Ectodermal Dysplasia in Fleckvieh Cattle

Mutations within the ectodysplasin A (EDA) gene have been associated with congenital hypotrichosis and anodontia (HAD/XHED) in humans, mice, dogs and cattle. We identified a three-generation family of Fleckvieh cattle with male calves exhibiting clinical and histopathological signs consistent with an X-linked recessive HAD (XHED). Whole genome and Sanger sequencing of cDNA showed a perfect association of the missense mutation g.85716041G>A (ss2019497443, rs1114816375) within the EDA gene with all three cases following an X-linked recessive inheritance, but normal EDAR and EDARADD. This mutation causes an exchange of glycine (G) with arginine (R) at amino acid position 227 (p.227G>R) in the second collagen triple helix repeat domain of EDA. The EDA variant was associated with a significant reduction and underdevelopment of hair follicles along with a reduced outgrowth of hairs, a complete loss of seromucous nasolabial and mucous tracheal and bronchial glands and a malformation of and reduction in number of teeth. Thermostability of EDA G227R was reduced, consistent with a relatively mild hair and tooth phenotype. However, incisors and canines were more severely affected in one of the calves, which correlated with the presence of a homozygous missense mutation of RNF111 (g.51306765T>G), a putative candidate gene possibly associated with tooth number in EDA-deficient Fleckvieh calves.

Large-scale detection and characterization of interchromosomal rearrangements in normozoospermic bulls using massive genotype and phenotype data sets

In this paper, we developed a highly sensitive approach to detect interchromosomal rearrangements in cattle by searching for abnormal linkage disequilibrium patterns between markers located on different chromosomes in large paternal half-sib families genotyped as part of routine genomic evaluations. We screened 5571 families of artificial insemination sires from 15 breeds and revealed 13 putative interchromosomal rearrangements, 12 of which were validated by cytogenetic analysis and long-read sequencing. These consisted of one Robertsonian fusion, 10 reciprocal translocations, and the first case of insertional translocation reported in cattle. Taking advantage of the wealth of data available in cattle, we performed a series of complementary analyses to define the exact nature of these rearrangements, investigate their origins, and search for factors that may have favored their occurrence. We also evaluated the risks to the livestock industry and showed significant negative effects on several traits in the sires and in their balanced or aneuploid progeny compared with wild-type controls. Thus, we present the most comprehensive and thorough screen for interchromosomal rearrangements compatible with normal spermatogenesis in livestock species. This approach is readily applicable to any population that benefits from large genotype data sets, and will have direct applications in animal breeding. Finally, it also offers interesting prospects for basic research by allowing the detection of smaller and rarer types of chromosomal rearrangements than GTG banding, which are interesting models for studying gene regulation and the organization of genome structure.

Review: Genetic mutations affecting bull fertility

Cattle are a well-suited "model organism" to study the genetic underpinnings of variation in male reproductive performance. The adoption of artificial insemination and genomic prediction in many cattle breeds provide access to microarray-derived genotypes and repeated measurements for semen quality and insemination success in several thousand bulls. Similar-sized mapping cohorts with phenotypes for male fertility are not available for most other species precluding powerful association testing. The repeated measurements of the artificial insemination bulls' semen quality enable the differentiation between transient and biologically relevant trait fluctuations, and thus, are an ideal source of phenotypes for variance components estimation and genome-wide association testing. Genome-wide case-control association testing involving bulls with either aberrant sperm quality or low insemination success revealed several causal recessive loss-of-function alleles underpinning monogenic reproductive disorders. These variants are routinely monitored with customised genotyping arrays in the male selection candidates to avoid the use of subfertile or infertile bulls for artificial insemination and natural service. Genome-wide association studies with quantitative measurements of semen quality and insemination success revealed quantitative trait loci for male fertility, but the underlying causal variants remain largely unknown. Moreover, these loci explain only a small part of the heritability of male fertility. Integrating genome-wide association studies with gene expression and other omics data from male reproductive tissues is required for the fine-mapping of candidate causal variants underlying variation in male reproductive performance in cattle.

Deleterious mutation load in the admixed mice population

Deleterious mutation loads are known to correlate negatively with effective population size (Ne). Due to this reason, previous studies observed a higher proportion of harmful mutations in small populations than that in large populations. However, the mutational load in an admixed population that derived from introgression between individuals from two populations with vastly different Ne is not known. We investigated this using the whole genome data from two subspecies of the mouse (Mus musculus castaneus and Mus musculus musculus) with significantly different Ne. We used the ratio of diversities at nonsynonymous and synonymous sites (dN/dS) to measure the harmful mutation load. Our results showed that this ratio observed for the admixed population was intermediate between those of the parental populations. The dN/dS ratio of the hybrid population was significantly higher than that of M. m. castaneus but lower than that of M. m. musculus. Our analysis revealed a significant positive correlation between the proportion of M. m. musculus ancestry in admixed individuals and their dN/dS ratio. This suggests that the admixed individuals with high proportions of M. m. musculus ancestry have large dN/dS ratios. We also used the proportion of deleterious nonsynonymous SNVs as a proxy for deleterious mutation load, which also produced similar results. The observed results were in concordance with those expected by theory. We also show a shift in the distribution of fitness effects of nonsynonymous SNVs in the admixed genomes compared to the parental populations. These findings suggest that the deleterious mutation load of the admixed population is determined by the proportion of the ancestries of the subspecies. Therefore, it is important to consider the status and the level of genetic admixture of the populations whilst estimating the mutation loads.

A large deletion encompassing exon 2 of the ectodysplasin A (EDA) gene in a British blue crossbred calf with hypohidrotic ectodermal dysplasia

Background

Hypohidrotic ectodermal dysplasia (HED) is a congenital syndrome of mammals affecting organs and tissues of ectodermal origin characterized by absence or hypoplasia of hair, teeth, and eccrine glands. The disorder has been reported in several species, including humans, mice, dogs and cattle, associated with variants in genes affecting the ectodysplasin pathway, including the X-linked ectodysplasin A (EDA) gene. Until now, nine pathogenic variants have been found in the bovine EDA gene. Here we report a novel variant in EDA in a crossbreed male Belgian Blue calf with HED, and provide an overview of the phenotypic and allelic heterogeneity of EDA-related forms of HED in cattle.

Case presentation

A 45-day-old male crossbreed British Blue calf was referred with congenital hypotrichosis, oligodontia and omphalitis. On histopathological examination of the nasal planum, nasolabial glands and ducts were not observed. The density of hair follicles was low, and they were small, with a predominance of telogen-phase hairs, and some serocellular crusts. The phenotype of the calf resembled that of HED. Whole-genome sequencing (WGS) was performed and revealed a 21,899 base-pair deletion encompassing the coding exon 2 of EDA, predicted to result in an altered transcript and aberrant protein.

Conclusions

The clinicopathological and genetic findings were consistent with a case of X-linked HED. A very similar EDA deletion has been previously reported in a family of Holstein cattle with HED. The newly identified hemizygous EDA loss-of-function variant is certainly pathogenic and therefore is the genetic cause for the observed phenotype. This case report provides an additional example of the potential of WGS-based precise diagnostics in livestock species such as cattle to increase the diagnostic yield in rare diseases.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13028-022-00641-2.

Locating a novel autosomal recessive genetic variant in the cattle glucokinase gene using only WGS data from three cases and six carriers

New Mendelian genetic conditions, which adversely affect livestock, arise all the time. To manage them effectively, some methods need to be devised that are quick and accurate. Until recently, finding the causal genomic site of a new autosomal recessive genetic disease has required a two-stage approach using single-nucleotide polymorphism (SNP) chip genotyping to locate the region containing the new variant. This region is then explored using fine-mapping methods to locate the actual site of the new variant. This study explores bioinformatic methods that can be used to identify the causative variants of recessive genetic disorders with full penetrance with just nine whole genome-sequenced animals to simplify and expedite the process to a one-step procedure. Using whole genome sequencing of only three cases and six carriers, the site of a novel variant causing perinatal mortality in Irish moiled calves was located. Four methods were used to interrogate the variant call format (VCF) data file of these nine animals, they are genotype criteria (GCR), autozygosity-by-difference (ABD), variant prediction scoring, and registered SNP information. From more than nine million variants in the VCF file, only one site was identified by all four methods (Chr4: g.77173487A>T (ARS-UCD1.2 (GCF_002263795.1)). This site was a splice acceptor variant located in the glucokinase gene (GCK). It was verified on an independent sample of animals from the breed using genotyping by polymerase chain reaction at the candidate site and autozygosity-by-difference using SNP-chips. Both methods confirmed the candidate site. Investigation of the GCR method found that sites meeting the GCR were not evenly spread across the genome but concentrated in regions of long runs of homozygosity. Locating GCR sites was best performed using two carriers to every case, and the carriers should be distantly related to the cases, within the breed concerned. Fewer than 20 animals need to be sequenced when using the GCR and ABD methods together. The genomic site of novel autosomal recessive Mendelian genetic diseases can be located using fewer than 20 animals combined with two bioinformatic methods, autozygosity-by-difference, and genotype criteria. In many instances it may also be confirmed with variant prediction scoring. This should speed-up and simplify the management of new genetic diseases to a single-step process.

Molecular genetic variation of animals and plants under domestication

Domesticated plants and animals played crucial roles as models for evolutionary change by means of natural selection and for establishing the rules of inheritance, originally proposed by Charles Darwin and Gregor Mendel, respectively. Here, we review progress that has been made during the last 35 y in unraveling the molecular genetic variation underlying the stunning phenotypic diversity in crops and domesticated animals that inspired Mendel and Darwin. We notice that numerous domestication genes, crucial for the domestication process, have been identified in plants, whereas animal domestication appears to have a polygenic background with no obvious "domestication genes" involved. Although model organisms, such as Drosophila and Arabidopsis, have replaced domesticated species as models for basic research, the latter are still outstanding models for evolutionary research because phenotypic change in these species represents an evolutionary process over thousands of years. A consequence of this is that some alleles contributing to phenotypic diversity have evolved by accumulating multiple changes in the same gene. The continued molecular characterization of crops and farm animals with ever sharper tools is essential for future food security.

Identification of mutant gene for Black crystal coat and non-allelic gene interactions in Neogale vison

Sable (Martes zibellina) and American mink (Neogale vison) are valuable species characterized by a variety of coat colour produced on fur farms. Black crystal fur phenotype is Mendelian codominant trait: heterozygous animals (C^r/ +) have white guard hairs scattered predominantly on the spine and the head, while homozygous (C^r/C^r) minks have coats resembling the Himalayan (c^h/c^h) or white Hedlund (h/h) types. It is one of the most recent of more than 35 currently known phenotypic traits of fur colour in American mink. Black crystal fur phenotype was first described in 1984 in the Russian population of mink, which had undergone selection for domestic defensive response to humans. Here, we performed whole-genome sequencing of American mink with C^r/C^r phenotype. We identified a missense mutation in the gene encoding the α-COP subunit of the COPI complex (COPA). The COPI complex mediates retrograde trafficking from the Golgi system to the endoplasmic reticulum and sorting of transmembrane proteins. We observed an interaction between a newly identified mutation in the COPA gene and a mutation in the microphthalmia-associated transcription factor (MITF), the latter mutation led to the formation of the white Hedlund (h/h) phenotype. Double heterozygotes for these mutations have an entirely white coat and a black-eyed phenotype similar to the phenotype of C^r/C^r or h/h minks. Our data could be useful for tracking economically valuable fur traits in mink breeding programs to contribute to global fur production.

Severe osteogenesis imperfecta caused by CREB3L1 mutation in a cat

We examined the clinical features and pathology, and identified the causative mutation, of osteogenesis imperfecta in a 2-mo-old kitten with growth retardation and abnormal gait. Blood and radiographic examinations were performed on presentation. Radiographs revealed decreased opacity of numerous bones. Fractures were observed in some long bones, including femur and tibia. Histologic examination of the tibia showed decreased osteoid and osteoblasts at the primary spongiosa extending from the growth plate. The periosteum was thickened, and cortical bone and osteoblasts were decreased. Consequently, osteogenesis imperfecta was diagnosed. Genomic DNA and total RNA were extracted from the skin and used for PCR. Whole-genome sequencing identified a 2-bp deletion (c.370_371delTG; p.C124fs), which resulted in a homozygous frameshift mutation on exon 3 of CREB3L1. This mutation introduced a premature stop codon, suggesting production of the truncated protein without a functional domain as a transcription factor for expression of COL1A1 mRNA. This error may have affected collagen fibril formation, leading to the development of osteogenesis imperfecta.

Identification of two new recessive MC1R alleles in red-coloured Evolèner cattle and other breeds

Sequence variations in the melanocortin-1 receptor (MC1R) gene are associated with melanism in different animal species. Six functionally relevant alleles have been described in cattle to date. In a hypothesis-free approach we performed a genome-wide allelic association study with black, red and wild-coloured cattle of three Alpine cattle breeds (Eringer, Evolèner and Valdostana), revealing a single significant association signal close to the MC1R gene. We searched for candidate causative variants by sequencing the entire coding sequence and identified two novel protein-changing variants. We propose designating the mutant alleles at MC1R:c.424C>T as e^v1 and at MC1R:c.263G>A as e^v2 . Both affect conserved amino acid residues in functionally important transmembrane domains (p.Arg142Cys and p.Ser88Asn). Both alleles segregate predominantly in the Swiss Evolèner breed. They occur in other European cattle breeds such as Abondance and Rotes Höhenvieh as well. We observed almost perfect association between the MC1R genotypes and the coat colour phenotype in a cohort of 513 black, red and wild-coloured cattle. Animals carrying two copies of MC1R loss-of-function alleles or that were compound heterozygous for e, e^v1 , or e^v2 have a red to dark red (chestnut-like red) coat colour. These findings expand the spectrum of causal MC1R variants causing recessive red in cattle.

A 1-bp deletion in bovine QRICH2 causes low sperm count and immotile sperm with multiple morphological abnormalities

BACKGROUND

Semen quality and insemination success are monitored in artificial insemination bulls to ensure high male fertility rates. Only ejaculates that fulfill minimum quality requirements are processed and eventually used for artificial inseminations. We examined 70,990 ejaculates from 1343 Brown Swiss bulls to identify bulls from which all ejaculates were rejected due to low semen quality. This procedure identified a bull that produced 12 ejaculates with an aberrantly small number of sperm (0.2 ± 0.2 × 109 sperm per mL) which were mostly immotile due to multiple morphological abnormalities.

RESULTS

The genome of this bull was sequenced at a 12× coverage to investigate a possible genetic cause. Comparing the sequence variant genotypes of this bull with those from 397 fertile bulls revealed a 1-bp deletion in the coding sequence of the QRICH2 gene which encodes the glutamine rich 2 protein, as a compelling candidate causal variant. This 1-bp deletion causes a frameshift in translation and a premature termination codon (ENSBTAP00000018337.1:p.Cys1644AlafsTer52). The analysis of testis transcriptomes from 76 bulls showed that the transcript with the premature termination codon is subject to nonsense-mediated mRNA decay. The 1-bp deletion resides in a 675-kb haplotype that includes 181 single nucleotide polymorphisms (SNPs) from the Illumina BovineHD Bead chip. This haplotype segregates at a frequency of 5% in the Brown Swiss cattle population. Our analysis also identified another bull that carried the 1-bp deletion in the homozygous state. Semen analyses from the second bull confirmed low sperm concentration and immotile sperm with multiple morphological abnormalities that primarily affect the sperm flagellum and, to a lesser extent, the sperm head.

CONCLUSIONS

A recessive loss-of-function allele of the bovine QRICH2 gene likely causes low sperm concentration and immotile sperm with multiple morphological abnormalities. Routine sperm analyses unambiguously identify homozygous bulls for this allele. A direct gene test can be implemented to monitor the frequency of the undesired allele in cattle populations.

Genomic landscape, polymorphism and possible LINE-associated delivery of G-quadruplex motifs in the bovine genes

G-Quadruplex structures are non-B DNA structures that occur in regions carrying short runs of guanines. They are implicated in several biological processes including transcription, translation, replication and telomere maintenance as well as in several pathological conditions like cancer and thus they have gained the attention of the scientific community. The rise of the -omics era significantly affected the G-quadruplex research and the genome-wide characterization of G-Quadruplexes has been rendered a necessary first step towards applying genomics approaches for their study. While in human and several model organisms there is a considerable number of works studying genome-wide the DNA motifs with potential to form G-quadruplexes (G4-motifs), there is a total absence of any similar studies regarding livestock animals. The objectives of the present study were to provide a detailed characterization of the bovine genic G4-motifs' distribution and properties and to suggest a possible mechanism for the delivery of G4 motifs in the genes. Our data indicate that the distribution of G4-motifs within bovine genes and the annotation of said genes to Gene Ontology terms are similar to what is already shown for other organisms. By investigating their structural characteristics and polymorphism, it is obvious that the overall stability of the putative quadruplex structures is in line with the current notion in the G4 field. Similarly to human, the bovine G4-motifs are overrepresented in specific LINE repeat elements, the L1_BTs in the case of cattle. We highlight the potential role of these elements as vehicles for delivery of G4 motifs in the introns of the bovine genes. Lastly, it seems that a basis exists for connecting traits of agricultural importance to the genetic variation of G4 motifs, thus, the value of cattle as an interesting new model organism for G4-related genetic studies might be worth to be investigated.

Non-additive QTL mapping of lactation traits in 124,000 cattle reveals novel recessive loci

BACKGROUND

Deleterious recessive conditions have been primarily studied in the context of Mendelian diseases. Recently, several deleterious recessive mutations with large effects were discovered via non-additive genome-wide association studies (GWAS) of quantitative growth and developmental traits in cattle, which showed that quantitative traits can be used as proxies of genetic disorders when such traits are indicative of whole-animal health status. We reasoned that lactation traits in cattle might also reflect genetic disorders, given the increased energy demands of lactation and the substantial stresses imposed on the animal. In this study, we screened more than 124,000 cows for recessive effects based on lactation traits.

RESULTS

We discovered five novel quantitative trait loci (QTL) that are associated with large recessive impacts on three milk yield traits, with these loci presenting missense variants in the DOCK8, IL4R, KIAA0556, and SLC25A4 genes or premature stop variants in the ITGAL, LRCH4, and RBM34 genes, as candidate causal mutations. For two milk composition traits, we identified several previously reported additive QTL that display small dominance effects. By contrasting results from milk yield and milk composition phenotypes, we note differing genetic architectures. Compared to milk composition phenotypes, milk yield phenotypes had lower heritabilities and were associated with fewer additive QTL but had a higher non-additive genetic variance and were associated with a higher proportion of loci exhibiting dominance.

CONCLUSIONS

We identified large-effect recessive QTL which are segregating at surprisingly high frequencies in cattle. We speculate that the differences in genetic architecture between milk yield and milk composition phenotypes derive from underlying dissimilarities in the cellular and molecular representation of these traits, with yield phenotypes acting as a better proxy of underlying biological disorders through presentation of a larger number of major recessive impacts.

DSP missense variant in a Scottish Highland calf with congenital ichthyosis, alopecia, acantholysis of the tongue and corneal defects

BACKGROUND

Ichthyosis describes a localized or generalized hereditary cornification disorder caused by an impaired terminal keratinocyte differentiation resulting in excessive stratum corneum with the formation of more or less adherent scales. Ichthyosis affects humans and animals. Two rare bovine forms are reported, the severe harlequin ichthyosis and the less severe congenital ichthyosis, both characterized by a severe orthokeratotic lamellar hyperkeratosis.

RESULTS

A 2-weeks-old purebred Scottish Highland calf was referred because of a syndrome resembling congenital ichthyosis. The clinical phenotype included diffuse alopecia and a markedly lichenified skin covered with large and excessive scales. Additionally, conjunctivitis and ulceration of the cornea were noted. Post-mortem examination revealed deep fissures in the diffusely thickened tongue and histopathological findings in the skin confirmed the clinical diagnosis. Whole-genome sequencing of the affected calf and comparison of the data with control genomes was performed. A search for private variants in known candidate genes for skin phenotypes including genes related with erosive and hyperkeratotic lesions revealed a single homozygous protein-changing variant, DSP: c.6893 C>A, or p.Ala2298Asp. The variant is predicted to change a highly conserved residue in the C-terminal plakin domain of the desmoplakin protein, which represents a main intracellular component of desmosomes, important intercellular adhesion molecules in various tissues including epidermis. Sanger sequencing confirmed the variant was homozygous in the affected calf and heterozygous in both parents. Further genotyping of 257 Scottish Highland animals from Switzerland revealed an estimated allele frequency of 1.2%. The mutant allele was absent in more than 4800 controls from various other cattle breeds.

CONCLUSIONS

This study represents the first report of combined lesions compatible with congenital ichthyosis, alopecia, acantholysis of the tongue and corneal defects associated with a DSP missense variant as the most likely underlying cause. To the best of our knowledge, this study is also the first report of a DSP-related syndromic form of congenital ichthyosis in domestic animals. The results of our study enable genetic testing to avoid the unintentional occurrence of further affected cattle. The findings were added to the Online Mendelian Inheritance in Animals (OMIA) database (OMIA 002243-9913).

The intronic branch point sequence is under strong evolutionary constraint in the bovine and human genome

The branch point sequence is a cis-acting intronic motif required for mRNA splicing. Despite their functional importance, branch point sequences are not routinely annotated. Here we predict branch point sequences in 179,476 bovine introns and investigate their variability using a catalogue of 29.4 million variants detected in 266 cattle genomes. We localize the bovine branch point within a degenerate heptamer "nnyTrAy". An adenine residue at position 6, that acts as branch point, and a thymine residue at position 4 of the heptamer are more strongly depleted for mutations than coding sequences suggesting extreme purifying selection. We provide evidence that mutations affecting these evolutionarily constrained residues lead to alternative splicing. We confirm evolutionary constraints on branch point sequences using a catalogue of 115 million SNPs established from 3,942 human genomes of the gnomAD database.

Identification of candidate genes on the basis of SNP by time-lagged heat stress interactions for milk production traits in German Holstein cattle

The aim of this study was to estimate genotype by time-lagged heat stress (HS) variance components as well as main and interaction SNP-marker effects for maternal HS during the last eight weeks of cow pregnancy, considering milk production traits recorded in the offspring generation. The HS indicator was the temperature humidity index (THI) for each week. A dummy variable with the code = 1 for the respective week for THI ≥ 60 indicated HS, otherwise, for no HS, the code = 0 was assigned. The dataset included test-day and lactation production traits from 14,188 genotyped first parity Holstein cows. After genotype quality control, 41,139 SNP markers remained for the genomic analyses. Genomic animal models without (model VC_nHS) and with in-utero HS effects (model VC_wHS) were applied to estimate variance components. Accordingly, for genome-wide associations, models GWA_nHS and GWA_wHS, respectively, were applied to estimate main and interaction SNP effects. Common genomic and residual variances for the same traits were very similar from models VC_nHS and VC_wHS. Genotype by HS interaction variances varied, depending on the week with in-utero HS. Among all traits, lactation milk yield with HS from week 5 displayed the largest proportion for interaction variances (0.07). For main effects from model GWA_wHS, 380 SNPs were suggestively associated with all production traits. For the SNP interaction effects from model GWA_wHS, we identified 31 suggestive SNPs, which were located in close distance to 62 potential candidate genes. The inferred candidate genes have various biological functions, including mechanisms of immune response, growth processes and disease resistance. Two biological processes excessively represented in the overrepresentation tests addressed lymphocyte and monocyte chemotaxis, ultimately affecting immune response. The modelling approach considering time-lagged genotype by HS interactions for production traits inferred physiological mechanisms being associated with health and immunity, enabling improvements in selection of robust animals.

Whole-genome resequencing reveals genetic structure and introgression in Pudong White pigs

Pudong White (PDW) pigs, historically originating from Shanghai, are the only Chinese indigenous pigs characterised by their completely white coats, with the exception of Rongchang pigs. However, there is limited information concerning their overall genetic structure or relationship with other breeds, especially the East Chinese (ECN) and European pigs. To uncover the genetic structure, selection signatures, and potential exotic introgression in PDW pigs, we sampled 15 PDW pigs using whole-genome sequencing (~20×). We then conducted in-depth population genetic analyses in 320 pigs from 27 global pig groups, namely, European wild boars, Chinese wild boars, and outgroup. Neighbour-joining tree and principal component analysis confirmed that PDW pigs belonged to the ecotype of ECN pigs. Both f3, D-statistics, and structure analysis showed that PDW pigs shared apparent alleles with Large White (LW) pigs. Three statistics, rIBD, a haplotype heat map and copy number variation, further indicated that PDW pigs shared apparent alleles with LW pigs at the KIT Proto-Oncogene, Receptor Tyrosine Kinase (KIT) and PARG-MARCHF8 loci, suggesting that the lineage of European pigs in PDW originated from LW pigs. After further detecting the KIT mutations in different pig breeds, PDW was confirmed to have the same duplication region 1, duplication region 2, and the splicing mutation on intron 17 of KIT as LW pigs that determine the white coat colour phenotype in European white pigs. We hypothesised that LW pigs were imported to China ∼110-160 years ago according to the admixture time estimate and then crossed with ECN pigs, resulting in the introgression of the KIT alleles that produce the white coat colour phenotype in the PDW pig breed. To our knowledge, this study presents the first thorough description of the genetic structure of PDW pigs via whole-genome resequencing data; moreover, the results provide a basis for the national project for the conservation of this unique Chinese local population.

Infertility due to defective sperm flagella caused by an intronic deletion in DNAH17 that perturbs splicing

Artificial insemination in pig (Sus scrofa domesticus) breeding involves the evaluation of the semen quality of breeding boars. Ejaculates that fulfill predefined quality requirements are processed, diluted and used for inseminations. Within short time, eight Swiss Large White boars producing immotile sperm that had multiple morphological abnormalities of the sperm flagella were noticed at a semen collection center. The eight boars were inbred on a common ancestor suggesting that the novel sperm flagella defect is a recessive trait. Transmission electron microscopy cross-sections revealed that the immotile sperm had disorganized flagellar axonemes. Haplotype-based association testing involving microarray-derived genotypes at 41,094 SNPs of six affected and 100 fertile boars yielded strong association (P = 4.22 × 10-15) at chromosome 12. Autozygosity mapping enabled us to pinpoint the causal mutation on a 1.11 Mb haplotype located between 3,473,632 and 4,587,759 bp. The haplotype carries an intronic 13-bp deletion (Chr12:3,556,401-3,556,414 bp) that is compatible with recessive inheritance. The 13-bp deletion excises the polypyrimidine tract upstream exon 56 of DNAH17 (XM_021066525.1: c.8510-17_8510-5del) encoding dynein axonemal heavy chain 17. Transcriptome analysis of the testis of two affected boars revealed that the loss of the polypyrimidine tract causes exon skipping which results in the in-frame loss of 89 amino acids from DNAH17. Disruption of DNAH17 impairs the assembly of the flagellar axoneme and manifests in multiple morphological abnormalities of the sperm flagella. Direct gene testing may now be implemented to monitor the defective allele in the Swiss Large White population and prevent the frequent manifestation of a sterilizing sperm tail disorder in breeding boars.

Conservation of Imprinting and Methylation of MKRN3, MAGEL2 and NDN Genes in Cattle

Genomic imprinting is the epigenetic mechanism of transcriptional regulation that involves differential DNA methylation modification. Comparative analysis of imprinted genes between species can help us to investigate the biological significance and regulatory mechanisms of genomic imprinting. MKRN3, MAGEL2 and NDN are three maternally imprinted genes identified in the human PWS/AS imprinted locus. This study aimed to assess the allelic expression of MKRN3, MAGEL2 and NDN and to examine the differentially methylated regions (DMRs) of bovine PWS/AS imprinted domains. An expressed single-nucleotide polymorphism (SNP)-based approach was used to investigate the allelic expression of MKRN3, MAGEL2 and NDN genes in bovine adult tissues and placenta. Consistent with the expression in humans and mice, we found that the MKRN3, MAGEL2 and NDN genes exhibit monoallelic expression in bovine somatic tissues and the paternal allele expressed in the bovine placenta. Three DMRs, PWS-IC, MKRN3 and NDN DMR, were identified in the bovine PWS/AS imprinted region by analysis of the DNA methylation status in bovine tissues using the bisulfite sequencing method and were located in the promoter and exon 1 of the SNRPN gene, NDN promoter and 5' untranslated region (5'UTR) of MKRN3 gene, respectively. The PWS-IC DMR is a primary DMR inherited from the male or female gamete, but NDN and MKRN3 DMR are secondary DMRs that occurred after fertilization by examining the methylation status in gametes.

Single nucleotide polymorphisms affect miRNA target prediction in bovine

Single nucleotide polymorphisms (SNPs) can have significant effects on phenotypic characteristics in cattle. MicroRNAs (miRNAs) are small, non-coding RNAs that act as post-transcriptional regulators by binding them to target mRNAs. In the present study, we scanned ~56 million SNPs against 1,064 bovine miRNA sequences and analyzed, in silico, their possible effects on target binding prediction, primary miRNA formation, association with QTL regions and the evolutionary conservation for each SNP locus. Following target prediction, we show that 71.6% of miRNA predicted targets were altered as a consequence of SNPs located within the seed region of the mature miRNAs. Next, we identified variations in the Minimum Free Energy (MFE), which represents the capacity to alter molecule stability and, consequently, miRNA maturation. A total of 48.6% of the sequences analyzed showed values within those previously reported as sufficient to alter miRNA maturation. We have also found 131 SNPs in 46 miRNAs, with altered target prediction, occurring in QTL regions. Lastly, analysis of evolutionary conservation scores for each SNP locus suggested that they have a conserved biological function through the evolutionary process. Our results suggest that SNPs in microRNAs have the potential to affect bovine phenotypes and could be of great value for genetic improvement studies, as well as production.

X-Linked Hypohidrotic Ectodermal Dysplasia in Crossbred Beef Cattle Due to a Large Deletion in EDA

Simple Summary

Ectodermal dysplasias such as hypohidrotic ectodermal dysplasia (HED), are genetic conditions affecting the development and/or homeostasis of two or more ectodermal derivatives, including hair, teeth, nails, and eccrine glands. In particular, X-linked hypohidrotic ectodermal dysplasia-1 (ECTD1) in humans is characterized by a triad of signs comprising sparse hair, abnormal teeth, and anhidrosis or hypohidrosis. It has been reported in cattle, dogs, mice and rats. Until now, eight pathogenic variants in the bovine ectodysplasin A (EDA) gene causing ECTD1-like disorders have been found. Herein, five affected Red Angus-Simmental bull calves born over a 6-year period (2013–2019) in a single herd in the Western United States are reported showing an ECTD1-like syndrome. Calves were born with severe hypotrichosis and oligodontia. At age 1-week-old two calves died of severe pneumonia. Microscopic findings of the skin revealed small-caliber hair follicles with a mean density in flank skin slightly greater in affected animals than in control animals. Nasolabial, intranasal and tracheobronchial mucosal glands were absent, whereas olfactory glands were unaffected. Whole-genome sequencing (WGS) identified a 53 kb deletion of the X chromosome including parts of the EDA gene as well as the entire AWAT2 gene. The partial deletion of the EDA gene that is known to be associated with forms of ECTD1 is the most likely cause for the reported genodermatosis. Similar rare disorders in livestock are often not diagnosed at the molecular level due to lack of resources, short lifespan of the animals, and concerns for the producers’ reputation.

Abstract

X-linked hypohidrotic ectodermal dysplasia-1 (ECTD1) in people results in a spectrum of abnormalities, most importantly hypotrichosis, anodontia/oligodontia, and absent or defective ectodermally derived glands. Five Red Angus-Simmental calves born over a 6-year period demonstrated severe hypotrichosis and were diagnosed as affected with ECTD1-like syndrome. Two died of severe pneumonia within a week of birth. The skin of three affected calves revealed a predominance of histologically unremarkable small-caliber hair follicles. Larger follicles (>50 µm) containing medullated hairs (including guard and tactile hairs) were largely restricted to the muzzle, chin, tail, eyelids, tragus and distal portions of the limbs and tail. The mean histological density of hair follicles in flank skin of two affected calves was slightly greater than that in two unaffected calves. One affected calf was examined postmortem at 10 days of age to better characterize systemic lesions. Nasolabial, intranasal and tracheobronchial mucosal glands were absent, whereas olfactory glands were unaffected. Mandibular incisor teeth were absent. Premolar teeth were unerupted and widely spaced. Other than oligodontia, histological changes in teeth were modest, featuring multifocal disorganization of ameloblasts, new bone formation in dental alveoli, and small aggregates of osteodentin and cementum at the margins of the enamel organ. A 52,780 base pair deletion spanning six out of eight coding exons of EDA and all of AWAT2 was identified. Partial deletion of the EDA gene is the presumed basis for the reported X-chromosomal recessive inherited genodermatosis.

A De Novo Mutation in COL1A1 in a Holstein Calf with Osteogenesis Imperfecta Type II

Osteogenesis imperfecta (OI) type II is a genetic connective tissue disorder characterized by bone fragility, severe skeletal deformities and shortened limbs. OI usually causes perinatal death of affected individuals. OI type II diagnosis in humans is established by the identification of heterozygous mutations in genes coding for collagens. The purpose of this study was to characterize the pathological phenotype of an OI type II-affected neonatal Holstein calf and to identify the causative genetic variant by whole-genome sequencing (WGS). The calf had acute as well as intrauterine fractures, abnormally shaped long bones and localized arthrogryposis. Genetic analysis revealed a private heterozygous missense variant in COL1A1 (c.3917T>A) located in the fibrillar collagen NC1 domain (p.Val1306Glu) that most likely occurred de novo. This confirmed the diagnosis of OI type II and represents the first report of a pathogenic variant in the fibrillar collagen NC domain of COL1A1 associated to OI type II in domestic animals. Furthermore, this study highlights the utility of WGS-based precise diagnostics for understanding congenital disorders in cattle and the need for continued surveillance for rare lethal genetic disorders in cattle.

Brain-specific monoallelic expression of bovine UBE3A is associated with genomic position

Genomic imprinting is a rare epigenetic process in mammalian cells that leads to monoallelic expression of a gene with a parent-specific pattern. The UBE3A (ubiquitin protein ligase E3A) gene is imprinted with maternal allelic expression in the brain but biallelically expressed in all other tissues in humans. The silencing of the paternal UBE3A allele is thought to be caused by the paternally expressed antisense RNA transcript of UBE3A-ATS. The aberrant imprinted expression of the UBE3A is associated with several neurodevelopmental syndromes and psychological disorders. Cattle are a valuable model species in determining the genetic etiology of sporadic human disorder, and maternal expression of UEB3A has been revealed by next-generation sequencing study in the bovine conceptus. In this study, we investigated the allelic expression of UBE3A and UBE3A-ATS in adult bovine somatic tissues. To confirm the splicing pattern of bovine UBE3A, five 5' alternative transcripts (MT210534-MT210538) were first obtained from bovine brain tissue by RT-PCR. Based on 10 SNP genotypes, we found that the brain-specific monoallelic expression of bovine UBE3A did not occur along the entire locus, and there was a shift from biallelic expression to monoallelic expression in exon 14 of the UBE3A gene. However, the brain-specific monoallelic expression of bovine UBE3A-ATS occurred in the entire gene. These observations demonstrated that the monoallelic expression did not occur along the bovine UBE3A entire locus and was associated with the genomic position.

A de novo mutation in KRT5 in a crossbred calf with epidermolysis bullosa simplex

A 6‐day‐old Belgian Blue‐Holstein calf was referred because of a syndrome resembling epidermolysis bullosa simplex (EBS). The clinical phenotype included irregular and differently sized erosions and ulcerations spread over the body, in particular on the limbs and over bone prominences, as well as in the nasal planum and oral mucosa. Blisters were easily induced by rubbing the skin. The skin lesions displayed a clear dermal‐epidermal separation at the level of the basal cell layer. Post mortem examination revealed erosions in the pharynx, proximal esophagus, and rumen. Whole‐genome sequencing revealed a heterozygous disruptive in‐frame deletion variant in KRT5 (c.534_536delCAA). Genotyping of both parents confirmed the variant as de novo mutation. Clinicopathological and genetic findings were consistent with the diagnosis of KRT5‐related EBS providing the second example of a spontaneous mutation causing epidermolysis bullosa in cattle.

A de novo frameshift mutation in ZEB2 causes polledness, abnormal skull shape, small body stature and subfertility in Fleckvieh cattle

Polledness in cattle is an autosomal dominant trait. Previous studies have revealed allelic heterogeneity at the polled locus and four different variants were identified, all in intergenic regions. In this study, we report a case of polled bull (FV-Polled1) born to horned parents, indicating a de novo origin of this polled condition. Using 50K genotyping and whole genome sequencing data, we identified on chromosome 2 an 11-bp deletion (AC_000159.1:g.52364063_52364073del; Del11) in the second exon of ZEB2 gene as the causal mutation for this de novo polled condition. We predicted that the deletion would shorten the protein product of ZEB2 by almost 91%. Moreover, we showed that all animals carrying Del11 mutation displayed symptoms similar to Mowat-Wilson syndrome (MWS) in humans, which is also associated with genetic variations in ZEB2. The symptoms in cattle include delayed maturity, small body stature and abnormal shape of skull. This is the first report of a de novo dominant mutation affecting only ZEB2 and associated with a genetic absence of horns. Therefore our results demonstrate undoubtedly that ZEB2 plays an important role in the process of horn ontogenesis as well as in the regulation of overall development and growth of animals.

Deletion of porcine BOLL is associated with defective acrosomes and subfertility in Yorkshire boars

A recessive sperm defect of Yorkshire boars was detected more than a decade ago. Affected boars produce ejaculates that contain spermatozoa with defective acrosomes, resulting in low fertility. The acrosome defect was mapped to porcine chromosome 15 but the causal mutation has not been identified. We re-analyzed microarray-derived genotypes of affected boars and confirmed that the acrosome defect maps to a 12.24 Mb segment of porcine chromosome 15. To detect the mutation causing defective acrosomes, we sequenced the genomes of two affected and three unaffected boars to an average coverage of 11-fold. Read depth analysis revealed a 55 kb deletion that is associated with the acrosome defect. The deletion encompasses the BOLL gene encoding the boule homolog, an RNA binding protein which is an evolutionarily conserved member of the DAZ (Deleted in AZoospermia) gene family. Lack of BOLL expression causes spermatogenic arrest and sperm maturation failure in many species. Boars that carry the deletion in the homozygous state produce sperm but their acrosomes are defective, suggesting that lack of porcine BOLL compromises acrosome formation. Our findings warrant further research to investigate the role of BOLL during spermatogenesis and sperm maturation in pigs.

A large deletion in the COL2A1 gene expands the spectrum of pathogenic variants causing bulldog calf syndrome in cattle

Background

Congenital bovine chondrodysplasia, also known as bulldog calf syndrome, is characterized by disproportionate growth of bones resulting in a shortened and compressed body, mainly due to reduced length of the spine and the long bones of the limbs. In addition, severe facial dysmorphisms including palatoschisis and shortening of the viscerocranium are present. Abnormalities in the gene collagen type II alpha 1 chain (COL2A1) have been associated with some cases of the bulldog calf syndrome. Until now, six pathogenic single-nucleotide variants have been found in COL2A1. Here we present a novel variant in COL2A1 of a Holstein calf and provide an overview of the phenotypic and allelic heterogeneity of the COL2A1-related bulldog calf syndrome in cattle.

Case presentation

The calf was aborted at gestation day 264 and showed generalized disproportionate dwarfism, with a shortened compressed body and limbs, and dysplasia of the viscerocranium; a phenotype resembling bulldog calf syndrome due to an abnormality in COL2A1. Whole-genome sequence (WGS) data was obtained and revealed a heterozygous 3513 base pair deletion encompassing 10 of the 54 coding exons of COL2A1. Polymerase chain reaction analysis and Sanger sequencing confirmed the breakpoints of the deletion and its absence in the genomes of both parents.

Conclusions

The pathological and genetic findings were consistent with a case of “bulldog calf syndrome”. The identified variant causing the syndrome was the result of a de novo mutation event that either occurred post-zygotically in the developing embryo or was inherited because of low-level mosaicism in one of the parents. The identified loss-of-function variant is pathogenic due to COL2A1 haploinsufficiency and represents the first structural variant causing bulldog calf syndrome in cattle. Furthermore, this case report highlights the utility of WGS-based precise diagnostics for understanding congenital disorders in cattle and the need for continued surveillance for genetic disorders in cattle.

A genome-wide scan for candidate lethal variants in Thoroughbred horses

Domestic animal populations are often characterised by high rates of inbreeding and low effective population sizes due to selective breeding practices. These practices can result in otherwise rare recessive deleterious alleles drifting to high frequencies, resulting in reduced fertility rates. This study aimed to identify potential recessive lethal haplotypes in the Thoroughbred horse breed, a closed population that has been selectively bred for racing performance. In this study, we identified a haplotype in the LY49B gene that shows strong evidence of being homozygous lethal, despite having high frequencies of heterozygotes in Thoroughbreds and other domestic horse breeds. Variant analysis of whole-genome sequence data identified two SNPs in the 3'UTR of the LY49B gene that may result in loss of function. Analysis of transcriptomic data from equine embryonic tissue revealed that LY49B is expressed in the trophoblast during placentation stage of development. These findings suggest that LY49B may have an essential, but as yet unknown function in the implantation stage of equine development. Further investigation of this region may allow for the development of a genetic test to improve fertility rates in horse populations. Identification of other lethal variants could assist in improving natural levels of fertility in horse populations.

Osteogenesis imperfecta in a male holstein calf associated with a possible oligogenic origin

Background: Neuromusculoskeletal anomalies generally in combination with severe clinical symptoms, comprise a heterogeneous group of fairly common and mostly fatal disorders in man and animals. Osteogenesis imperfecta (OI), also known as brittle bone disease, causes bone fragility and deformity. Prominent extra-skeletal accessory manifestations of OI comprise blue/gray sclerae, hearing impairment, lung abnormalities and hypercalciuria. Cases of OI in cattle have been reported. However, no causative mutations have been identified in cattle so far.Aim: To report a possible oligogenic origin identified in a calf from clinically healthy parents suffering from OI.Materials and Methods: A neonatal embryo transfer male Holstein calf developing multiple fractures with bone tissue showing marked osteopenia was used for whole genome re-sequencing as well as its parents. In addition, 2,612 randomly chosen healthy Holstein cattle were genotyped as well as controls.Results: Sixteen candidate genes with potential protein-altering variants were selected revealing non-synonymous variants only within IFITM5 and CRTAP genes. However, in-depth gene analysis did not result in the identification of a single causative mutation in the OI calf.Conclusion: The analysis of the OI case revealed a possible oligogenic origin of the disease attributable to additive effects of three candidate genes, i.e., ABCA13, QRFPR, and IFTIM5.Clinical relevance: Most OI cases in humans and domestic animals reported so far are caused by distinct dominant or recessive monogenic mutations, therefore a potential oligogenic additive genetic effect is a novel finding. Furthermore, the case presented here demonstrates that cross-species genetic analyses might not always be straightforward.

ZNF597 is a maternally expressed imprinted gene in the Holstein breed

Genomic imprinting is an epigenetic phenomenon that leads to the preferential expression of genes from either the paternal or maternal allele. Imprinted genes play important roles in mammalian growth and development and a central role in placental function. ZNF597 and NAA60 are two paternally imprinted genes in the human ZNF597-NAA60 imprinted locus, both of which show biallelic expression in the mouse, but their imprinting status in cattle is still unknown. In this study, we examined the allelic expression of ZNF597 and NAA60 in adult bovine placental and somatic tissues. By comparing the mRNA-based genotypes with the genomic DNA-based genotypes, we identified monoallelic expression of ZNF597 in the placenta and in seven other tissues, including the cerebrum, heart, liver, spleen, lung, kidney, and muscle. Nevertheless, analysis revealed biallelic expression of the NAA60 gene in these tissues. Moreover, we tested the imprinting status of ZNF597 and confirmed that the maternal allele is expressed in the bovine placenta. To determine the role of DNA methylation in regulating monoallelic/imprinted expression of bovine ZNF597, the methylation status of two CpG-enriched regions in the bovine ZNF597-NAA60 locus was analyzed using the bisulfite sequencing method. Differentially methylated regions were detected on ten CpG loci in the bovine ZNF597 promoter region. In summary, the bovine ZNF597 gene is a maternally expressed gene, and its expression is regulated by DNA methylation, whereas the NAA60 gene is not imprinted in cattle.

Comprehensive Geno- and Phenotyping in a Complex Pedigree Including Four Different Inherited Retinal Dystrophies

Inherited retinal dystrophies (IRDs) are characterized by high clinical and genetic heterogeneity. A precise characterization is desirable for diagnosis and has impact on prognosis, patient counseling, and potential therapeutic options. Here, we demonstrate the effectiveness of the combination of in-depth retinal phenotyping and molecular genetic testing in complex pedigrees with different IRDs. Four affected Caucasians and two unaffected relatives were characterized including multimodal retinal imaging, functional testing, and targeted next-generation sequencing. A considerable intrafamilial phenotypic and genotypic heterogeneity was identified. While the parents of the index family presented with rod-cone dystrophy and ABCA4-related retinopathy, their two sons revealed characteristics in the spectrum of incomplete congenital stationary night blindness and ocular albinism, respectively. Molecular testing revealed previously described variants in RHO, ABCA4, and MITF as well as a novel variant in CACNA1F. Identified variants were verified by intrafamilial co-segregation, bioinformatic annotations, and in silico analysis. The coexistence of four independent IRDs caused by distinct mutations and inheritance modes in one pedigree is demonstrated. These findings highlight the complexity of IRDs and underscore the need for the combination of extensive molecular genetic testing and clinical characterization. In addition, a novel variant in the CACNA1F gene is reported associated with incomplete congenital stationary night blindness.

1000 Bull Genomes Project to Map Simple and Complex Genetic Traits in Cattle: Applications and Outcomes

The 1000 Bull Genomes Project is a collection of whole-genome sequences from 2,703 individuals capturing a significant proportion of the world's cattle diversity. So far, 84 million single-nucleotide polymorphisms (SNPs) and 2.5 million small insertion deletions have been identified in the collection, a very high level of genetic diversity. The project has greatly accelerated the identification of deleterious mutations for a range of genetic diseases, as well as for embryonic lethals. The rate of identification of causal mutations for complex traits has been slower, reflecting the typically small effect size of these mutations and the fact that many are likely in as-yet-unannotated regulatory regions. Both the deleterious mutations that have been identified and the mutations associated with complex trait variation have been included in low-cost SNP array designs, and these arrays are being genotyped in tens of thousands of dairy and beef cattle, enabling management of deleterious mutations in these populations as well as genomic selection.

A comprehensive biomedical variant catalogue based on whole genome sequences of 582 dogs and eight wolves

The domestic dog serves as an excellent model to investigate the genetic basis of disease. More than 400 heritable traits analogous to human diseases have been described in dogs. To further canine medical genetics research, we established the Dog Biomedical Variant Database Consortium (DBVDC) and present a comprehensive list of functionally annotated genome variants that were identified with whole genome sequencing of 582 dogs from 126 breeds and eight wolves. The genomes used in the study have a minimum coverage of 10× and an average coverage of ~24×. In total, we identified 23 133 692 single-nucleotide variants (SNVs) and 10 048 038 short indels, including 93% undescribed variants. On average, each individual dog genome carried ∼4.1 million single-nucleotide and ~1.4 million short-indel variants with respect to the reference genome assembly. About 2% of the variants were located in coding regions of annotated genes and loci. Variant effect classification showed 247 141 SNVs and 99 562 short indels having moderate or high impact on 11 267 protein-coding genes. On average, each genome contained heterozygous loss-of-function variants in 30 potentially embryonic lethal genes and 97 genes associated with developmental disorders. More than 50 inherited disorders and traits have been unravelled using the DBVDC variant catalogue, enabling genetic testing for breeding and diagnostics. This resource of annotated variants and their corresponding genotype frequencies constitutes a highly useful tool for the identification of potential variants causative for rare inherited disorders in dogs.

Genome-wide association analysis reveals QTL and candidate mutations involved in white spotting in cattle

BACKGROUND

White spotting of the coat is a characteristic trait of various domestic species including cattle and other mammals. It is a hallmark of Holstein-Friesian cattle, and several previous studies have detected genetic loci with major effects for white spotting in animals with Holstein-Friesian ancestry. Here, our aim was to better understand the underlying genetic and molecular mechanisms of white spotting, by conducting the largest mapping study for this trait in cattle, to date.

RESULTS

Using imputed whole-genome sequence data, we conducted a genome-wide association analysis in 2973 mixed-breed cows and bulls. Highly significant quantitative trait loci (QTL) were found on chromosomes 6 and 22, highlighting the well-established coat color genes KIT and MITF as likely responsible for these effects. These results are in broad agreement with previous studies, although we also report a third significant QTL on chromosome 2 that appears to be novel. This signal maps immediately adjacent to the PAX3 gene, which encodes a known transcription factor that controls MITF expression and is the causal locus for white spotting in horses. More detailed examination of these loci revealed a candidate causal mutation in PAX3 (p.Thr424Met), and another candidate mutation (rs209784468) within a conserved element in intron 2 of MITF transcripts expressed in the skin. These analyses also revealed a mechanistic ambiguity at the chromosome 6 locus, where highly dispersed association signals suggested multiple or multiallelic QTL involving KIT and/or other genes in this region.

CONCLUSIONS

Our findings extend those of previous studies that reported KIT as a likely causal gene for white spotting, and report novel associations between candidate causal mutations in both the MITF and PAX3 genes. The sizes of the effects of these QTL are substantial, and could be used to select animals with darker, or conversely whiter, coats depending on the desired characteristics.

Α de novo 3.8-Mb inversion affecting the EDA and XIST genes in a heterozygous female calf with generalized hypohidrotic ectodermal dysplasia

BACKGROUND

In mammals, hypohidrotic ectodermal dysplasia (HED) is a genetic disorder that is characterized by sparse hair, tooth abnormalities, and defects in cutaneous glands. Only four genes, EDA, EDAR, EDARADD and WNT10A account for more than 90% of HED cases, and EDA, on chromosome X, is involved in 50% of the cases. In this study, we explored an isolated case of a female Holstein calf with symptoms similar to HED.

RESULTS

Clinical examination confirmed the diagnosis. The affected female showed homogeneous hypotrichosis and oligodontia as previously observed in bovine EDAR homozygous and EDA hemizygous mutants. Under light microscopy, the hair follicles were thinner and located higher in the dermis of the frontal skin in the affected animal than in the control. Moreover, the affected animal showed a five-fold increase in the number of hair follicles and a four-fold decrease in the diameter of the pilary canals. Pedigree analysis revealed that the coefficient of inbreeding of the affected calf (4.58%) was not higher than the average population inbreeding coefficient (4.59%). This animal had ten ancestors in its paternal and maternal lineages. By estimating the number of affected cases that would be expected if any of these common ancestors carried a recessive mutation, we concluded that, if they existed, other cases of HED should have been reported in France, which is not the case. Therefore, we assumed that the causal mutation was dominant and de novo. By analyzing whole-genome sequencing data, we identified a large chromosomal inversion with breakpoints located in the first introns of the EDA and XIST genes. Genotyping by PCR-electrophoresis the case and its parents allowed us to demonstrate the de novo origin of this inversion. Finally, using various sources of information we present a body of evidence that supports the hypothesis that this mutation is responsible for a skewed inactivation of X, and that only the normal X can be inactivated.

CONCLUSIONS

In this article, we report a unique case of X-linked HED affected Holstein female calf with an assumed full inactivation of the normal X-chromosome, thus leading to a severe phenotype similar to that of hemizygous males.

Evidence for a de novo, dominant germ-line mutation causative of osteogenesis imperfecta in two Red Angus calves

A genetic disorder, osteogenesis imperfecta (OI) is broadly characterized by connective tissue abnormalities and bone fragility most commonly attributed to alterations in Type I collagen. Two Red Angus calves by the same sire presented with severe bone and dental fragility, blue sclera, and evidence of in utero fractures consistent with OI congenita. Comparative analyses with human cases suggested the OI in these calves most closely resembled that classified as OI Type II. Due to the phenotypic classification and shared paternity, a dominant, germ-line variant was hypothesized as causative although recessive genotypes were also considered due to a close relationship between the sire and dam of one calf. Whole-genome sequencing revealed the presence of a missense mutation in the alpha 1 chain of collagen Type I (COL1A1), for which both calves were heterozygous. The variant resulted in the substitution of a glycine residue with serine in the triple helical domain of the protein; in this region, glycine normally occupies every third position as is critical for correct formation of the Type I collagen molecule. Allele-specific amplification by droplet digital PCR further quantified the variant at a frequency of nearly 4.4% in the semen of the sire while it was absent in his blood, supporting the hypothesis of a de novo causative variant for which the germ line of the sire was mosaic. The identification of novel variants associated with unwanted phenotypes in livestock is critical as the high prolificacy of breeding stock has the potential to rapidly disseminate undesirable variation.

A recessive lethal chondrodysplasia in a miniature zebu family results from an insertion affecting the chondroitin sulfat domain of aggrecan

Background

Congenital skeletal malformations represent a heterogeneous group of disorders affecting bone and cartilage development. In cattle, particular chondrodysplastic forms have been identified in several miniature breeds. In this study, a phenotypic characterization was performed of an affected Miniature Zebu calf using computed tomography, necropsy and histopathological examinations, whole genome sequencing of the case and its parents on an Illumina NextSeq 500 in 2 × 150 bp paired-end mode and validation using Sanger sequencing and a Kompetitive Allele Specific PCR assay. Samples from the family of an affected Miniature Zebu with bulldog syndrome including parents and siblings, 42 healthy Miniature Zebu not related with members of the herd and 88 individuals from eight different taurine cattle breeds were available for validation.

Results

A bulldog-like Miniature Zebu calf showing a large bulging head, a short and compressed body and extremely short and stocky limbs was delivered after a fetotomy. Computed tomography and necropsy revealed severe craniofacial abnormalities including a shortening of the ventral nasal conchae, a cleft hard palate, rotated limbs as well as malformed and fused vertebrae and ribs. Histopathologic examination showed a disorganization of the physeal cartilage with disorderly arranged chondrocytes in columns and a multifocal closed epiphyseal plate. Whole-genome sequencing of this malformed Miniature Zebu calf, its dam and sire and subsequent comparative sequence analysis revealed a one base pair insertion (ACAN:c.5686insC) located within the cartilage development gene aggrecan (ACAN) exclusively homozygous in the affected calf and heterozygous in its parents. This variant was predicted to cause a frameshift (p.Val1898fsTer9) and thus a truncation of the chondroitin sulfate domain as well as a loss of the C-terminal globular domain of ACAN. It perfectly co-segregated with the lethal bulldog syndrome in Miniature Zebus.

Conclusions

We found a novel mutation in ACAN causing a recessive lethal chondrodysplasia in Miniature Zebu cattle. A diagnostic test for this mutation is now available for Miniature Zebu breeders preventing further cases of bulldog syndrome by targeted matings. To the authors’ best knowledge, this is the first case of a Miniature Zebu associated with an ACAN mutation.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0678-8) contains supplementary material, which is available to authorized users.

Genome variants associated with RNA splicing variations in bovine are extensively shared between tissues

Background

Mammalian phenotypes are shaped by numerous genome variants, many of which may regulate gene transcription or RNA splicing. To identify variants with regulatory functions in cattle, an important economic and model species, we used sequence variants to map a type of expression quantitative trait loci (expression QTLs) that are associated with variations in the RNA splicing, i.e., sQTLs. To further the understanding of regulatory variants, sQTLs were compare with other two types of expression QTLs, 1) variants associated with variations in gene expression, i.e., geQTLs and 2) variants associated with variations in exon expression, i.e., eeQTLs, in different tissues.

Results

Using whole genome and RNA sequence data from four tissues of over 200 cattle, sQTLs identified using exon inclusion ratios were verified by matching their effects on adjacent intron excision ratios. sQTLs contained the highest percentage of variants that are within the intronic region of genes and contained the lowest percentage of variants that are within intergenic regions, compared to eeQTLs and geQTLs. Many geQTLs and sQTLs are also detected as eeQTLs. Many expression QTLs, including sQTLs, were significant in all four tissues and had a similar effect in each tissue. To verify such expression QTL sharing between tissues, variants surrounding (±1 Mb) the exon or gene were used to build local genomic relationship matrices (LGRM) and estimated genetic correlations between tissues. For many exons, the splicing and expression level was determined by the same cis additive genetic variance in different tissues. Thus, an effective but simple-to-implement meta-analysis combining information from three tissues is introduced to increase power to detect and validate sQTLs. sQTLs and eeQTLs together were more enriched for variants associated with cattle complex traits, compared to geQTLs. Several putative causal mutations were identified, including an sQTL at Chr6:87392580 within the 5th exon of kappa casein (CSN3) associated with milk production traits.

Conclusions

Using novel analytical approaches, we report the first identification of numerous bovine sQTLs which are extensively shared between multiple tissue types. The significant overlaps between bovine sQTLs and complex traits QTL highlight the contribution of regulatory mutations to phenotypic variations.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4902-8) contains supplementary material, which is available to authorized users.