Deletion of one of the duplicated Hsp70 genes causes hereditary myopathy of diaphragmatic muscles in Holstein-Friesian cattle

A 44-kb deleted-type copy number variation is associated with decreasing complement component activity and calf mortality in Japanese Black cattle

BACKGROUND

Calf mortality generally occurs in calves prior to weaning, which is a serious problem in cattle breeding. Several causative variants of monogenic Mendelian disorders in calf mortality have been identified, whereas genetic factors affecting the susceptibility of calves to death are not well known. To identify variants associated with calf mortality in Japanese Black cattle, we evaluated calf mortality as a categorical trait with a threshold model and performed a genome-wide copy number variation (CNV) association study on calf mortality.

RESULTS

We identified a 44-kb deleted-type CNV ranging from 103,317,687 to 103,361,802 bp on chromosome 5, which was associated with the mortality of 1-180-day-old calves. The CNV harbored C1RL, a pseudogene, and an IncRNA localized in the C1R and C1S gene cluster, which is a component of the classical complement activation pathway for immune complexes for infectious pathogens. The average complement activity in CNVR_221 homozygotes at postnatal day 7 was significantly lower than that of wild-type animals and heterozygotes. The frequency of the risk allele in dead calves suffering from diarrhea and pneumonia and in healthy cows was 0.35 and 0.28, respectively (odds ratio = 2.2, P = 0.016), suggesting that CNVR_221 was associated with the mortality of Japanese Black calves suffering from an infectious disease.

CONCLUSIONS

This study identified a deleted-type CNV associated with the mortality of 1-180-day-old calves. The complement activity in CNVR_221 homozygotes was significantly lower than that in heterozygotes and wild type animals. The frequency of the risk allele was higher in dead calves suffering from an infectious disease than in healthy cows. These results suggest that the existence of CNVR_221 in calves could be attributed to a reduction in complement activity, which in turn leads to susceptibility to infections. Thus, the risk allele could serve as a useful marker to reduce the mortality of infected Japanese Black calves.

Stress Effects on Meat Quality: A Mechanistic Perspective

Stress inevitably occurs from the farm to abattoir in modern livestock husbandry. The effects of stress on the behavioral and physiological status and ultimate meat quality have been well documented. However, reports on the mechanism of stress effects on physiological and biochemical changes and their consequent effects on meat quality attributes have been somewhat disjointed and limited. Furthermore, the causes of variability in meat quality traits among different animal species, muscle fibers within an animal, and even positions within a piece of meat in response to stress are still not entirely clear. This review 1st summarizes the primary stress factors, including heat stress, preslaughter handling stress, oxidative stress, and other stress factors affecting animal welfare; carcass quality; and eating quality. This review further delineates potential stress-induced pathways or mediators, including AMP-activated protein kinase-mediated energy metabolism, crosstalk among calcium signaling pathways and reactive oxygen species, protein modification, apoptosis, calpain and cathepsin proteolytic systems, and heat shock proteins that exert effects that cause biochemical changes during the early postmortem period and affect the subsequent meat quality. To obtain meat of high quality, further studies are needed to unravel the intricate mechanisms involving the aforementioned signaling pathways or mediators and their crosstalk.

Characterization of a deletion in the Hsp70 cluster in the bovine reference genome

The 70 kilodalton heat shock proteins (Hsp70) are highly conserved molecular chaperones which have a crucial role in the stress response of the cell. In mammals, the Hsp70 proteins are encoded by a cluster of three genes: HSPA1A, HSPA1B and HSPA1L. In bovines, this cluster is located on chromosome 23 downstream of the major histocompatibility complex (BoLA). We detected inconsistencies in the location of markers on the Hsp70 genes reported in the literature that pointed to a potential deletion in the bovine reference genome UMD 3.1.1. An in silico analysis of the bovine genomic region of the Hsp70 cluster, using available information from public databases, confirmed the existence of a deletion of 11.1-kb spanning the HSPA1B gene and the intergenic region between HSPA1B and HSPA1A. Although we originally considered this an assembly error, it is most likely a particular condition of L1 Dominette 01449, the cow sequenced in the Bovine Genome Project. Moreover, we suggest a new classification of bovine Hsp70 sequences reported in NCBI and a reassignment of the location of SNPs from dbSNP that map to the deletion on BTA23. We also compared the location of selected transcription factor binding sites on the promoters of HSPA1A and HSPA1B. The results generated in the present work could be helpful to refine the reference genome of an important livestock species and also to understand the role and the regulation of the bovine Hsp70 genes.

Mild Aerobic Exercise Training Hardly Affects the Diaphragm of mdx Mice

In the mdx mice model of Duchenne Muscular Dystrophy (DMD), mild endurance exercise training positively affected limb skeletal muscles, whereas few and controversial data exist on the effects of training on the diaphragm. The diaphragm was examined in mdx (C57BL/10ScSn-Dmdmdx) and wild-type (WT, C57BL/10ScSc) mice under sedentary conditions (mdx-SD, WT-SD) and during mild exercise training (mdx-EX, WT-EX). At baseline, and after 30 and 45 days (training: 5 d/wk for 6 weeks), diaphragm muscle morphology and Cx39 protein were assessed. In addition, tissue levels of the chaperonins Hsp60 and Hsp70 and the p65 subunit of nuclear factor-kB (NF-kB) were measured in diaphragm, gastrocnemius, and quadriceps in each experimental group at all time points. Although morphological analysis showed unchanged total area of necrosis/regeneration in the diaphragm after training, there was a trend for larger areas of regeneration than necrosis in the diaphragm of mdx-EX compared to mdx-SD mice. However, the levels of Cx39, a protein associated with active regeneration in damaged muscle, were similar in the diaphragm of mdx-EX and mdx-SD mice. Hsp60 significantly decreased at 45 days in the diaphragm, but not in limb muscles, in both trained and sedentary mdx compared to WT mice. In limb muscles, but not in the diaphragm, Hsp70 and NF-kB p65 levels were increased in mdx mice irrespective of training at 30 and 45 days. Therefore, the diaphragm of mdx mice showed little inflammatory and stress responses over time, and appeared hardly affected by mild endurance training. J. Cell. Physiol. 232: 2044-2052, 2017. © 2016 Wiley Periodicals, Inc.

Loss of maternal ANNEXIN A10 via a 34-kb deleted-type copy number variation is associated with embryonic mortality in Japanese Black cattle

Background

Conception is a fundamental trait for successful cattle reproduction. However, conception rates in Japanese Black cattle have been gradually declining over the last two decades. Although conception failures are mainly caused by embryonic mortality, the role of maternal genetic factors in the process remains unknown. Copy number variation (CNV), defined as large-scale genomic structural variants, contributes to several genetic disorders. To identify CNV associated with embryonic mortality in Japanese Black cattle, we evaluated embryonic mortality as a categorical trait with a threshold model and conducted a genome-wide CNV association study for embryonic mortality using 791 animals.

Results

We identified a deleted-type CNV ranging from 378,127 to 412,061 bp on bovine chromosome 8, which was associated with embryonic mortality at 30–60 days after artificial insemination (AI). The CNV harbors exon 2 to 6 of ANNEXIN A10 (ANXA10). Analysis of sequence traces from the CNV identified that 63 bp reads bridging the breakpoint were present on both sides of the CNV, indicating that the CNV was generated by non-allelic homologous recombination using the 63 bp homologous sequences. Western blot analysis showed that the CNV results in a null allele of ANXA10. This association was replicated using a sample population size of 2552 animals. To elucidate the function of ANXA10 in vivo, we generated Anxa10 null mice using the CRISPR/Cas9 system. Crossbreeding experiments showed that litter size from crosses of both Anxa10^-/- and Anxa10^+/- females had fewer pups than did Anxa10^+/+ females, and embryos of Anxa10^-/- females died between implantation stages E4.5 and E12.5. These results indicate that loss of maternal Anxa10 causes embryonic mortality.

Conclusions

This study identified a deleted-type CNV encompassing ANXA10 in cows that was associated with embryonic mortality at 30–60 days after AI. Using a mouse model, we confirmed that litter sizes were smaller in crosses of both Anxa10^-/- and Anxa10^+/- females relative to those of wild females. These results indicate that ANXA10 is a maternal factor that is critical for embryo development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3312-z) contains supplementary material, which is available to authorized users.

Genome-wide patterns of copy number variation in the Chinese yak genome

Background

Copy number variation (CNV) represents an important source of genetic divergence that can produce drastic phenotypic differences and may therefore be subject to selection during domestication and environmental adaptation. To investigate the evolutionary dynamics of CNV in the yak genome, we used a read depth approach to detect CNV based on genome resequencing data from 14 wild and 65 domestic yaks and determined CNV regions related to domestication and adaptations to high-altitude.

Results

We identified 2,634 CNV regions (CNVRs) comprising a total of 153 megabases (5.7 % of the yak genome) and 3,879 overlapping annotated genes. Comparison between domestic and wild yak populations identified 121 potentially selected CNVRs, harboring genes related to neuronal development, reproduction, nutrition and energy metabolism. In addition, we found 85 CNVRs that are significantly different between domestic yak living in high- and low-altitude areas, including three genes related to hypoxia response and six related to immune defense. This analysis shows that genic CNVs may play an important role in phenotypic changes during yak domestication and adaptation to life at high-altitude.

Conclusions

We present the first refined CNV map for yak along with comprehensive genomic analysis of yak CNV. Our results provide new insights into the genetic basis of yak domestication and adaptation to living in a high-altitude environment, as well as a valuable genetic resource that will facilitate future CNV association studies of important traits in yak and other bovid species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2702-6) contains supplementary material, which is available to authorized users.

Genome-wide identification of copy number variation using high-density single-nucleotide polymorphism array in Japanese Black cattle

Background

Copy number variation (CNV) is an important source of genetic variability associated with phenotypic variation and disease susceptibility. Comprehensive genome-wide CNV maps provide valuable information for genetic and functional studies. To identify CNV in Japanese Black cattle, we performed a genome-wide autosomal screen using genomic data from 1,481 animals analyzed with the Illumina Bovine High-Density (HD) BeadChip Array (735,293 single-nucleotide polymorphisms (SNPs) with an average marker interval of 3.4 kb on the autosomes).

Results

We identified a total of 861 CNV regions (CNVRs) across all autosomes, which covered 43.65 Mb of the UMD3.1 genome assembly and corresponded to 1.74 % of the 29 bovine autosomes. Overall, 35 % of the CNVRs were present at a frequency of > 1 % in 1,481 animals. The estimated lengths of CNVRs ranged from 1.1 kb to 1.4 Mb, with an average of 50.7 kb. The average number of CNVR events per animal was 35. Comparisons with previously reported cattle CNV showed that 72 % of the CNVR calls detected in this study were within or overlapped with known CNVRs. Experimentally, three CNVRs were validated using quantitative PCR, and one CNVR was validated using PCR with flanking primers for the deleted region. Out of the 861 CNVRs, 390 contained 717 Ensembl-annotated genes significantly enriched for stimulus response, cellular defense response, and immune response in the Gene Ontology (GO) database. To associate genes contained in CNVRs with phenotypes, we converted 560 bovine Ensembl gene IDs to their 438 orthologous associated mouse gene IDs, and 195 of these mouse orthologous genes were categorized into 1,627 phenotypes in the Mouse Genome Informatics (MGI) database.

Conclusions

We identified 861 CNVRs in 1,481 Japanese Black cattle using the Illumina BovineHD BeadChip Array. The genes contained in CNVRs were characterized using GO analysis and the mouse orthologous genes were characterized using the MGI database. The comprehensive genome-wide CNVRs map will facilitate identification of genetic variation and disease-susceptibility alleles in Japanese Black cattle.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0335-z) contains supplementary material, which is available to authorized users.

Genome-Wide Study of Structural Variants in Bovine Holstein, Montbéliarde and Normande Dairy Breeds

High-throughput sequencing technologies have offered in recent years new opportunities to study genome variations. These studies have mostly focused on single nucleotide polymorphisms, small insertions or deletions and on copy number variants. Other structural variants, such as large insertions or deletions, tandem duplications, translocations, and inversions are less well-studied, despite that some have an important impact on phenotypes. In the present study, we performed a large-scale survey of structural variants in cattle. We report the identification of 6,426 putative structural variants in cattle extracted from whole-genome sequence data of 62 bulls representing the three major French dairy breeds. These genomic variants affect DNA segments greater than 50 base pairs and correspond to deletions, inversions and tandem duplications. Out of these, we identified a total of 547 deletions and 410 tandem duplications which could potentially code for CNVs. Experimental validation was carried out on 331 structural variants using a novel high-throughput genotyping method. Out of these, 255 structural variants (77%) generated good quality genotypes and 191 (75%) of them were validated. Gene content analyses in structural variant regions revealed 941 large deletions removing completely one or several genes, including 10 single-copy genes. In addition, some of the structural variants are located within quantitative trait loci for dairy traits. This study is a pan-genome assessment of genomic variations in cattle and may provide a new glimpse into the bovine genome architecture. Our results may also help to study the effects of structural variants on gene expression and consequently their effect on certain phenotypes of interest.

Inverse relationships between biomarkers and beef tenderness according to contractile and metabolic properties of the muscle

Previous proteomic analyses established a list of proteins biomarkers of beef tenderness. The present study quantified the relative abundance of 21 of these proteins by dot-blot technique in the Longissimus thoracis and Semitendinosus muscles of 71 young bulls from three breeds: Aberdeen Angus (AA), Limousin (LI), and Blond d'Aquitaine (BA). For both muscles overall tenderness was estimated by sensory analysis; shear force was measured with a Warner-Bratzler instrument, and an index combining sensory and mechanical measurements was calculated. Multiple regressions based on relative abundances of these proteins were used to propose equations of prediction of the three evaluations of tenderness. Hsp70-1B appeared to be a good biomarker of low tenderness in the three breeds and in the two muscles. Proteins such as lactate dehydrogenase-B, myosin heavy chain IIx, and small heat shock proteins (Hsp27, Hsp20, and αB-crystallin) were related to tenderness but inversely according to the muscle and breed. The results demonstrate that prediction of tenderness must take into account muscle characteristics and animal type.

Identification and differential expression of human carcinoma-associated antigens in hepatocellular carcinoma tissues

This study was designed to identify and verify hepatocellular carcinoma (HCC)-associated human carcinoma antigens (HCAs) that may be useful as tumor markers for HCC. We found that BCE075 and BCD021 anti-HCA antibodies were immunostained in the liver tissue samples and showed specific staining. Their expression was increased in HCC compared with normal liver tissues (P = 0.008). Immunoprecipitation and mass spectrometry analyses of the proteins precipitated by these two antibodies were identified to be cytoskeleton-associated protein 4 (CLIMP63) and brain-type glycogen phosphorylase (PYGB). This study demonstrated that HCC tissues expressed specific HCA glycoproteins, suggesting that our mouse monoclonal anti-HCA antibodies could be useful for immunohistochemical analysis of HCA expression as potential biomarkers for HCC diagnosis.

[Bovine disease-related DNA mutations and their genetic control strategies in breeding for disease resistance]

Bovine genomic DNA mutations and their genetic effects on gene expression and protein function influence disease susceptibility and resistance of cattle. The genetic loci related to cattle diseases are mainly divided into two types: single-locus-disease genes and multigenic-disease loci. The single-locus-disease genes are called causal mutations; their genetic basis is simply and normally detected in the coding and non-coding regions inducing substitution of amino acid, premature termination of translation, and complete deletion of entire exon(s). In contrast, the genetic basis of disease related to multiple genes is more complex since susceptibility or resistance of these diseases is affected by the interactions among host, pathogen, and environment. This article reviewed current research and application of the major diseases of cattle con-trolled by single gene or polygenic genes. The genetic control strategies of effective identification and control of these dis-eases in bovine breeding and production were also analyzed.

Variant in the 5' untranslated region of insulin-like growth factor 1 receptor is associated with susceptibility to mastitis in cattle

Mastitis is a common infectious disease of the mammary gland and generates large losses in the dairy industry. By means of positional cloning and functional analysis techniques, we here show that insulin-like growth factor 1 receptor (IGF1R) can possibly mediate susceptibility to mastitis through autophagy. Scanning the whole genome of cows (Bos taurus) that were susceptible or resistant to mastitis in the half-sib families revealed that susceptible cows had a relatively long stretch of cytosine residues (C stretch) in the 5' untranslated region of IGF1R. The forebrain embryonic zinc finger-like (FEZL) transcription factor, which was previously identified as a factor controlling mastitis resistance in the same half-sib families, bound the C stretch of IGF1R. The susceptible type of FEZL with a glycine stretch containing 13 glycines (13G) and the longer C stretch of IGF1R together enhanced expression of IGF1R. Enhancing IGF1R inhibited autophagy in response to Streptococcus agalactiae invasion of mammary epithelial cells, whereas treatment with rapamycin, a known inducer of autophagy, rescued it. Cows carrying the variant combination of 13GFEZL might be more susceptible to mastitis as the result of impaired autophagy. Our results suggest that IGF1R could control innate immunity in mammals and serve as a potential tool for preventing mastitis.

Identification and functional characterization of an aggregation domain in long myosin light chain kinase

The functions of long smooth muscle myosin light chain kinase (L-MLCK), a molecule with multiple domains, are poorly understood. To examine the existence of further potentially functional domains in this molecule, we analyzed its amino acid sequence with a tango program and found a putative aggregation domain located at the 4Ig domain of the N-terminal extension. To verify its aggregation capability in vitro, expressible truncated L-MLCK variants driven by a cytomegalovirus promoter were transfected into cells. As anticipated, only the overexpression of the 4Ig fragment led to particle formation in Colon26 cells. These particles contained 4Ig polymers and actin. Analysis with detergents demonstrated that the particles shared features in common with aggregates. Thus, we conclude that the 4Ig domain has a potent aggregation ability. To further examine this aggregation domain in vivo, eight transgenic mouse lines expressing the 4Ig domain (4Ig lines) were generated. The results showed that the transgenic mice had typical aggregation in the thigh and diaphragm muscles. Histological examination showed that 7.70 +/- 1.86% of extensor digitorum longus myofibrils displayed aggregates with a 36.44% reduction in myofibril diameter, whereas 65.13 +/- 3.42% of diaphragm myofibrils displayed aggregates and the myofibril diameter was reduced by 43.08%. Electron microscopy examination suggested that the aggregates were deposited at the mitochondria, resulting in structural impairment. As a consequence, the oxygen consumption of mitochondria in the affected muscles was also reduced. Macrophenotypic analysis showed the presence of muscular degeneration characterized by a reduction in force development, faster fatigue, decreased myofibril diameters, and structural alterations. In summary, our study revealed the existence of a novel aggregation domain in L-MLCK and provided a direct link between L-MLCK and aggregation. The possible significance and mechanism underlying the aggregation-based pathological processes mediated by L-MLCK are also discussed.

A critical analysis of disease-associated DNA polymorphisms in the genes of cattle, goat, sheep, and pig

Genetic variations through their effects on gene expression and protein function underlie disease susceptibility in farm animal species. The variations are in the form of single nucleotide polymorphisms, deletions/insertions of nucleotides or whole genes, gene or whole chromosomal rearrangements, gene duplications, and copy number polymorphisms or variants. They exert varying degrees of effects on gene action, such as substitution of an amino acid for another, shift in reading frame and premature termination of translation, and complete deletion of entire exon(s) or gene(s) in diseased individuals. These factors influence gene function by affecting mRNA splicing pattern or by altering/eliminating protein function. Elucidating the genetic bases of diseases under the control of many genes is very challenging, and it is compounded by several factors, including host x pathogen x environment interactions. In this review, the genetic variations that underlie several diseases of livestock (under monogenic and polygenic control) are analyzed. Also, factors hampering research efforts toward identification of genetic influences on animal disease identification and control are highlighted. A better understanding of the factors analyzed could be better harnessed to effectively identify and control, genetically, livestock diseases. Finally, genetic control of animal diseases can reduce the costs associated with diseases, improve animal welfare, and provide healthy animal products to consumers, and should be given more attention.

YAC/BAC contig spanning the MHC class III region of cattle

A contig of the class III region of the bovine major histocompatibility complex (MHC) was established from bacterial and yeast artificial chromosomes using PCR and BAC-end sequencing. The marker content of individual clones was determined by gene and BAC-end specific PCR, and the location of genes and BAC-ends was confirmed analyzing somatic hybrid cells. A comparative analysis indicated that the content and order of MHC class III genes is strongly conserved between cattle and other mammalian species. Fluorescence in situ hybridization localized the bovine class III region to BTA23q21-->q22. The results show that the collection of sequenced BAC-ends is a powerful resource for generating high-resolution comparative chromosome maps.

Heat shock proteins and neuromuscular disease

The heat shock proteins are families of proteins with known activities that include chaperoning nascent peptides within the cell and cytoprotection. Most work on the nervous system has related to the role of heat shock proteins in neuroprotection from either hypoxic-ischemic or traumatic injury. The role of these proteins during normal physiological activity and injury is still under investigation. Heat shock proteins in neuromuscular disease have been investigated to some extent but were largely neglected until recently. The goal of this review is to summarize the evidence linking heat shock proteins with neuromuscular disease and to provide some insight into the roles or functions of these proteins in disease states.