- Research article
- Open Access
A 44-kb deleted-type copy number variation is associated with decreasing complement component activity and calf mortality in Japanese Black cattle
BMC Genomics volume 22, Article number: 107 (2021)
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.
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.
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.
Calf health is one of the most important components in cattle production, both in terms of animal welfare and profitability. In both dairy and beef cattle, however, calf mortality, which generally ranges from 2 to 15%, mainly occurs in calves prior to weaning [1,2,3,4]. In Japanese Black cattle, the calf mortality rate during the first 180 days, at which point the calves are weaned, is 3.3% (16,292 dead animals/year) ; this has mainly been attributed to gastrointestinal and respiratory infectious diseases, such as diarrhea and pneumonia . Recently, several studies discovered a causal variant of single gene disorders (monogenic Mendelian disorders) in calf mortality using single-nucleotide polymorphism arrays and next-generation sequencing technologies [7,8,9,10,11,12], whereas genetic factors affecting calves’ susceptibility to disease and infections and their association with calf mortality are unknown. Information on the molecular mechanisms underlying the role of infectious diseases in calf mortality provides clues both for the treatment of the disease and breeding, thereby effectively leading to improved calf health.
Copy number variation (CNV) is defined as deletion or duplication regions of the genome that range from 1 kb to several Mb, as reviewed in . CNVs cause diseases through various mechanisms, such as gene dosage alternation and gene structure disturbance, as reviewed in . We previously reported genome-wide CNV maps for Japanese Black cattle  and found a CNV related to embryonic mortality in Japanese Black cattle . Thus, the CNV maps also enabled us to search the CNVs associated with the mortality of Japanese Black calves.
In current study, we evaluated calf mortality as a categorical trait with a threshold model [17, 18] and performed a genome-wide CNV association study using these traits. We identified a deleted-type CNV that was associated with decreased complement component activity and calf mortality at the first 180 days (6 months) after birth.
Results and discussion
CNVR_221 on chromosome 5 was associated with calf mortality in Japanese Black cattle
The mortality of 1–180-day-old Japanese Black calves was previously estimated from 40,412 calf records in a threshold model using Bayesian analysis with Gibbs sampling [18,19,20]. The mortality rate of 1–180-day-old calves was 10.33%, and the estimates of direct heritability of calf mortality were 0.165 .
We detected 861 CNV regions (CNVRs), defined as the union area of overlapping CNVs identified in at least two animals , in the autosomes of 1481 Japanese Black cattle . Of the 1481 animals, 791 cows were evaluated for the breeding value of calf mortality (calves that died aged 1–180 d) (Additional file 1). Previous reports have indicated that deleted-type CNVs are connected to diseases in cattle [22,23,24,25,26]; thus, in this study, we evaluated the association between deleted-type CNVRs and calf mortality. Of the 861 CNVRs, 116 deleted-type CNVRs (Additional file 2: Table S1), with a minor allele frequency greater than 0.01 in 1481 animals  (Additional file 2: Table S1), were used in the association study. To identify CNVR associated with calf mortality, we performed a genome-wide CNV association study with the breeding value of calf mortality of 791 cows, in GEMMA software . We found that CNVR_221 was associated with calf mortality, reaching the Bonferroni-corrected threshold (P < 4.31 × 10 − 4, Fig. 1a). CNVR_221 was located within a 44,115 bp window from 103,317,687 to 103,361,802 bp on bovine chromosome 5 (P = 2.94 × 10− 4) (Table 1). Of note, the surrounding SNPs at a distance of approximately ±200 kb from CNVR_221 were not associated with the traits studied (Fig. 1b) and were not in linkage disequilibrium (LD) with CNVR_221 (Additional file 2: Table S2).
The mean log R ratio of six SNPs, which were continuously located within a 44,115 bp window between 103,317,687 bp and 103,361,802 bp on chromosome 5, was low (Fig. 2a). Whole genome resequencing and quantitative polymerase chain reaction (PCR) (qPCR) were conducted to verify CNVR calls using the SNP array as two independent experimental validations. Alignment of sequence reads resulted in the disappearance of read counts in the CNVR_221 region from five CNVR_221 homozygotes (Fig. 2c). The copy number calculated by qPCR was approximately one and zero in the CNVR_221 heterozygotes and the homozygotes, respectively (Fig. 2d), which concurred with the expected copy number estimated by the SNP array.
We found that CNVR_221 harbored three transcripts: C1RL (complement C1r subcomponent: LOC112446756), pseudogene (LOC783686), and IncRNA (LOC101902570) (Fig. 2b, Additional file 2: Table S3). CNVR was localized on the C1R and C1S gene cluster (Fig. 2b, Additional file 2: Table S3), which is a component of the classical complement activation pathway for immune complexes for infectious pathogens. C1R and C1S have two CUB, a single EGF, and two CCP modules followed by a serine protease domain  (Additional file 2: Table S4). The structure of C1RL (complement C1r subcomponent like: LOC112446756, XM_024992686.1 was the Refseq in model) was similar to that of C1R, whereas C1RL had no CCP modules and serine protease domain (Additional file 2: Table S4). Although the expression of the C1RL transcript was confirmed by RNAseq data (Additional file 3) , it is uncertain whether C1RL is a functional protein involved in complement activity.
CNVR_221 created a C1RL (Complement C1r subcomponent like: LOC112446756) null allele
CNVR_221, which contained the full length of C1RL (LOC112446756) (Fig. 2b), resulted in the C1RL null allele (Fig. 2c, d). To examine the effect of the 44 kb deletion on the expressions, we conducted RT-qPCR using total RNA from the livers of cattle, as complement genes are mainly expressed in the liver [30,31,32]. The level of expression in the heterozygous animals was significantly lower than that in the wild-type animals (P = 0.002) (Fig. 3a). The expression in the homozygous animals was undetectable (Fig. 3a), suggesting that the C1RL (LOC112446756) mRNA was expressed in the liver of wild-type animals and C1RL (LOC112446756) mRNA derived from the CNVR_221 allele could be not expressed; thus, CNVR_221 generated null alleles of C1RL (LOC112446756).
C1R was closely located on the telomeric side at a distance of approximately 15 kbp from CNVR_221 (Fig. 2b), which could include regulatory elements of C1R expression. The level of expression in the CNVR_221 homozygotes was not significantly different from that in heterozygotes and wild-type animals (Fig. 3b), indicating that CNVR_221 did not affect C1R expression.
The loss of C1RL (LOC112446756) in calves decreased the complement activity
C1r and C1s, the serine protease subunit, are components of the classical complement activation pathway, reviewed in . The activation of the classical pathway is primed by antigen-antibody complexes. The complexes trigger the assembly and activation of the C1 complex, C1q–C1r–C1s, causing the activation of the C1r and C1s. Following the activation of the components, later components of the complement cascade from the membrane attack complex (MAC), which creates membrane pores for cell lysis and facilitates phagocytosis through opsonization. Previous studies in humans reported that C1r and C1s deficient patients had recurrent infections, and result in death at juvenile stages due to severe infection [33,34,35,36]. These results suggest that the deficiency of C1RL (LOC112446756) in CNVR_221 calves could be explained by a reduction in complement activity, which in turn leads to recurrent infections.
To determine whether the complement activity was affected in CNVR_221 calves, the complement activity was examined at postnatal day 7 to exclude as much as possible the effects of environmental factors after birth. The average complement activity in CNVR_221 homozygotes at postnatal day 7 was 0.7- and 0.75-fold lower than that of wild animals (P < 0.01) and heterozygotes (P < 0.05), respectively (Fig. 4). The assembly of the C1s-C1r-C1r-C1s involves Ca2+-dependent C1r-C1s interactions between the CUB-EGF-CUB modules [28, 37,38,39]. Although C1RL (LOC112446756) only has CUB-EGF-CUB modules (Additional file 2: Table S4), the protein could affect the assembly of the complement complex via an unknown mechanism. These findings suggest that the complement activity in CNVR_221 homozygous calves slightly decreases and might lead to susceptibility to disease and infection. As the effect of the loss of C1RL (LOC112446756) on the complement activity is currently unknown, further research into the functions of C1RL (LOC112446756) on complement activity is required.
The frequency of CNVR_221 in dead calves that suffered from infectious disease
To determine the CNVR_221 associated with infectious disease and calf mortality, we genotyped 246 dead calves that had suffered from diarrhea or pneumonia, respectively (Table S5), and 287 healthy cows that had given birth at least once, as control animals. To evaluate the population stratification, principal component analysis (PCA) was applied to the numerator relationship matrix of dead calves and healthy cows. The PCA plot showed no obvious distinction between the dead calves and the healthy cows (Additional file 4), suggesting that there was no population structure and none of the tests led to a bias for the genotyping analysis. The frequency of CNVR_221 in dead calves that suffered from an infectious disease and in healthy cows was 0.35 and 0.28, respectively (odds ratio = 2.2, Table 2). Additionally, we detected an excess of CNVR_221 in dead calves compared to healthy cows (P = 0.016; two-sided Fisher’s exact test, Table 2), suggesting that CNVR_221 was associated with mortality in Japanese Black calves suffering from an infectious disease. Notably, the CNVR_221 homozygous healthy cows were found in the population at a frequency of 0.08 (Table 2), indicating that the homozygote had survival ability under certain conditions; thus, CNVR_221 could be a quantitative trait locus (QTL) affecting the susceptibility of infected calves to death. In addition to the classical pathway, the complement pathway is also activated throughout lectin and alternative pathways, as reviewed in ; thus, some individuals with classical complement pathway deficiency could be rescued from infectious diseases by other complement pathways. As it is not known whether a slight reduction in the complement pathway in CNVR_221 homozygotes is directly connected to the onset of diarrhea or pneumonia, further investigation into the relationship between the complement activity and the CNVR_221 genotype, and the pathogenesis of infectious disease in calves will be required in many cases.
In this study, we evaluated calf mortality as a categorical trait with a threshold model and searched the CNVs associated with calf mortality in Japanese Black cattle. We identified a deleted-type CNVR, with decreasing complement activity, which was associated with calf mortality and susceptibility to infection. The mortality rate of Japanese Black calves during the first 180 days was 3.3%, which was mainly attributed to gastrointestinal and respiratory infectious diseases. Thus, CNVR could serve as a useful marker of the mortality of Japanese Black calves suffering from infections.
All animal experiments were conducted according to the guidelines for the care and use of laboratory animals of the Shirakawa Institute of Animal Genetics, Northern Center of Agricultural Technology, General Technological Center of Hyogo Prefecture for Agriculture, and the University of the Ryukyus (A2018128). We obtained written agreement from the cattle farms to use the samples and data. The husbandry conditions of the animals were according to animal welfare guidelines in Japan for beef cattle . A total of 44,001 Japanese Black cattle were used in this study, which were obtained from the Shirakawa Institute of Animal Genetics, Okayama University, and Hokubu Agricultural Technology Institute. After collecting the blood samples, the animals were released and retained in the farms.
Phenotypic data and breeding value of calf mortality
The data consisted of 40,412 records of animals . The threshold statistical model for calf mortality included the fixed effects of sex, parity, breeding farm, and combination of mating year and month as covariances:
where Uijkl is an unobservable underlying continuous variable (liability), μ is the overall mean, SEXi is the fixed effect of sex, PARITYj is the fixed effect of parity, FARMk is the fixed effect of farm, YMl is the fixed effect of month and year of birth, aijkl is the direct additive genetic effect as a random effect, and eijkl is the residual effect. Yijkl is the observation ijkl for calf mortality; the trait was treated as a categorical trait: if the calf died, it was assigned a value of 2; if the cow lived it was assigned a value of 1. The observed binary response was assigned values of 1 and 0 if the liability was above a fixed threshold (t) and below t, respectively. Genetic parameters were estimated using the Gibbs sampling procedure. The THRGIBBS1F90 program was used to fit a threshold animal model for embryonic mortality [19, 20].
SNP genotyping and CNVR detection
Blood was collected from 791 Japanese Black cows. Genomic DNA was extracted using the Easy-DNA kit (Invitrogen, Cat. #K1800–01). All samples were genotyped using the BovineHD BeadChip Array (Illumina, Cat. #WG-450-1002) . SNPs genotype calling were performed using GenomeStudio V2011 (Illumina, version 1.9.4). The ARS-UCD1.2 assembly was used to map the SNP position . CNVR was detected as previously reported , using PennCNV software (version June 2011) .
Genome-wide CNV association study for calf mortality
Among the 861 CNVRs, 116 deleted-type CNVRs (Additional file 2: Table S1), which constituted a minor allele frequency greater than 0.01 in 1481 animals  (Additional file 2: Table S1), were used in the association study. Association study was performed on the breeding value of calf mortality of 791 cows using GEMMA software . SNP association analyses were performed for 791 samples using SNPs in the BovineHD BeadChip Array. The estimated power in the genome-wide association study (GWAS) was 0.88 (Sample size = 791, Heritability = 0.165 , a= 0.000431, variance of the off-diagonals of the grm = 0.0026) using GCTA-GREML Power Calculator .
Whole genome resequencing
Genomic DNA was sheared using a DNA shearing system (Covaris, Cat. #S220) prior to library preparation. Genomic DNA (200 ng) from six Japanese Black cattle (one wild/wild homozygote and five CNVR_221/CNVR_221 homozygotes) was used to generate sequence libraries, which were prepared using TruSeq Nano DNA library prep kit (Illumina, Cat. #20015964) and IDT for Illumina TruSeq DNA UD Indexes (Illumina, Cat. # 20020590) according to the manufacturer’s instructions. The sequence data were generated as 2 × 101 bp reads using the NovaSeq 6000 SP Reagent Kit (Illumina, Cat. # 20027465) and Novaseq6000, and processing and base calling were performed using Illumina Real-Time Analysis 3.
The sequence reads (fastq file) for each animal were aligned to ARS-UCD1.2 genome assembly  using bwa-0.7.17 . The bam files were created using Samtools [46, 47], and duplicate reads were removed using Picard . The realignment around indels and recalibration were performed using GATK-220.127.116.11 . SNPs and indels were called using Haplotype caller by GATK. The data from the WGS were deposited in the Wagyu genome database (WGDB) of the Japan Livestock Technology Association (Yushima, Bunyouku, Tokyo 113–0034, Japan) and were managed by the WGDB consortium.
Real-time quantitative PCR validation of CNVR
For real-time quantitative PCR (qPCR), we extracted genomic DNA from the blood of 28 Japanese Black cattle (two wild/wild homozygotes, 13 wild/CNVR_221 heterozygotes, and 13 CNVR_221/CNVR_221 homozygotes). Real-time qPCR was performed for CNVR validation using an ABI StepOnePlus Real Time PCR system (Applied Biosystems, Cat. #4376600). Primers and probes were designed for CNVR_221 (Additional file 2: Table S6). The basic transcription factor 3 gene (BTF3), which served as an internal qPCR standard for both copies at a locus (2n) [15, 50].
Expression analysis using qPCR
We collected the livers of Japanese Black cattle (five wild/wild homozygotes, seven wild/CNVR_221 heterozygotes, and seven CNVR_221/CNVR_221 homozygotes) from a slaughterhouse (Hyogo prefecture). For RT qPCR, we isolated the total RNA from the livers using RNeasy mini kits (QIAGEN, Cat. #74104). cDNA was synthesized from 50 ng RNA using ReverTra Ace-α with random primers (TOYOBO Cat. #FSK-101). Bovine C1RL (LOC112446756) and C1R were amplified using the primers described in Additional file 1: Table S7. PCR products using the primer pairs for C1RL (LOC112446756) and C1R were cloned into the pCRII-TOPO vector (Invitrogen, Cat. #K460040), and the sequences were confirmed using M13_R primers. Real-time PCR was performed on an ABI StepOnePlus Real Time PCR system using the comparative Ct method with glyceraldehyde-3-phosphate dehydrogenase (GAPD) as the internal control.
Measurement of complement activity
To determine whether the complement activities were affected in CNVR_221 calves, blood sera were collected postnatally (day 7) from 25 Japanese Black cattle (five wild/wild homozygotes, six wild/CNVR_221 heterozygotes, and 14 CNVR_221/CNVR_221 homozygotes). The sera were stored at 80 °C until the assay was performed. The 50% hemolytic complement activity  was measured in the serum using a complement activity assay kit (Denka-Seiken, Cat. #CH50) according to the manufacturer’s instructions. Hemolytic sheep erythrocytes, which were pre-coated with anti-sheep erythrocyte antibody, were measured at 541 nm with a plate spectrometer (BioRad, Cat. #iMark).
CNVR_221 frequency in the Japanese Black cattle population
For the surveillance of the allele frequency, we randomly selected 246 dead calves that had suffered from diarrhea or pneumonia, as diagnosed by veterinarians. We used 287 healthy cows that had given birth at least once as control animals. To evaluate the population stratification, PCA was applied to the numerator relationship matrix. Real-time qPCR was performed for CNVR_221 genotyping using an ABI StepOnePlus Real Time PCR system. Primers and probes were designed for CNVR_221 (Additional file 2: Table S6).
Availability of data and materials
The data sets supporting the results of this article are included within the article and its additional file. The data from WGS that support the findings of this study are available from the Wagyu genome database (WGDB) of the Japan Livestock Technology Association (Yushima, Bunyouku, Tokyo 113–0034, Japan) and were managed by the WGDB consortium. Data are available from the authors upon reasonable request and with permission from WGDB. Refseq ID and protein name are based on Reference Sequence Database (RefSeq) of the National Center of Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov/refseq/).
Copy number variation
Genetic relationship matrix
Genome-wide association study
Membrane attack complex
Principal component analysis
Polymerase chain reaction
Single nucleotide polymorphism
Cundiff LV, MacNeil MD, Gregory KE, Koch RM. Between- and within-breed genetic analysis of calving traits and survival to weaning in beef cattle. J Anim Sci. 1986;63(1):27–33.
Erf DF, Hansen LB, Neitzel RR. Inheritance of calf mortality for Brown Swiss cattle. J Dairy Sci. 1990;73(4):1130–4.
Azzam SM, Kinder JE, Nielsen MK, Werth LA, Gregory KE, Cundiff LV, Koch RM. Environmental effects on neonatal mortality of beef calves. J Anim Sci. 1993;71(2):282–90.
Hansen M, Madsen P, Jensen J, Pedersen J, Christensen LG. Genetic parameters of postnatal mortality in Danish Holstein calves. J Dairy Sci. 2003;86(5):1807–17.
National Livestock Breeding Center: Cattle Individual Identification Number Information: Notification information. http://www.id.nlbc.go.jp/data/toukei.html (2017) Accessed 13 Jul 2017.
Ministry of Agriculture Forestry and Fisheries: Statistical Tables of Livestock Mutual Relief (in Japanese). http://www.maff.go.jp/j/tokei/kouhyou/katiku_kyosai/ (2017) Accessed 13 Jul 2017.
Jung S, Pausch H, Langenmayer MC, Schwarzenbacher H, Majzoub-Altweck M, Gollnick NS, Fries R. A nonsense mutation in PLD4 is associated with a zinc deficiency-like syndrome in Fleckvieh cattle. BMC Genomics. 2014;15:623.
Sartelet A, Druet T, Michaux C, Fasquelle C, Geron S, Tamma N, Zhang Z, Coppieters W, Georges M, Charlier C. A splice site variant in the bovine RNF11 gene compromises growth and regulation of the inflammatory response. PLoS Genet. 2012;8(3):e1002581.
Hirano T, Kobayashi N, Matsuhashi T, Watanabe D, Watanabe T, Takasuga A, Sugimoto M, Sugimoto Y. Mapping and exome sequencing identifies a mutation in the IARS gene as the cause of hereditary perinatal weak calf syndrome. PLoS One. 2013;8(5):e64036.
Venhoranta H, Pausch H, Flisikowski K, Wurmser C, Taponen J, Rautala H, Kind A, Schnieke A, Fries R, Lohi H, et al. In frame exon skipping in UBE3B is associated with developmental disorders and increased mortality in cattle. BMC Genomics. 2014;15:890.
Menzi F, Besuchet-Schmutz N, Fragniere M, Hofstetter S, Jagannathan V, Mock T, Raemy A, Studer E, Mehinagic K, Regenscheit N, et al. A transposable element insertion in APOB causes cholesterol deficiency in Holstein cattle. Anim Genet. 2016;47(2):253–7.
Kipp S, Segelke D, Schierenbeck S, Reinhardt F, Reents R, Wurmser C, Pausch H, Fries R, Thaller G, Tetens J, et al. Identification of a haplotype associated with cholesterol deficiency and increased juvenile mortality in Holstein cattle. J Dairy Sci. 2016;99(11):8915–31.
Feuk L, Carson AR, Scherer SW. Structural variation in the human genome. Nat Rev Genet. 2006;7(2):85–97.
Bickhart DM, Liu GE. The challenges and importance of structural variation detection in livestock. Front Genet. 2014;5:37.
Sasaki S, Watanabe T, Nishimura S, Sugimoto Y. Genome-wide identification of copy number variation using high-density single-nucleotide polymorphism array in Japanese black cattle. BMC Genet. 2016;17(1):26.
Sasaki S, Ibi T, Akiyama T, Fukushima M, Sugimoto Y. Loss of maternal ANNEXIN A10 via a 34-kb deleted-type copy number variation is associated with embryonic mortality in Japanese black cattle. BMC Genomics. 2016;17(1):968.
Wright S. An analysis of variability in number of digits in an inbred strain of Guinea pigs. Genetics. 1934;19(6):506–36.
Ibi T, Kahi AK, Hirooka H. Genetic parameters of postnatal mortality and birth weight in Japanese black calves. Anim Sci J. 2014;86(1):25–30.
Lee D, Misztal I, Bertrand JK, Rekaya R. National evaluation for calving ease, gestation length and birth weight by linear and threshold model methodologies. J Appl Genet. 2002;43(2):209–16.
Tsuruta S, Misztal I. THRGIBBS1F90 for estimation of variance components with threshold-linear model. In: Proceedings of 8th World Congress on Genetics Applied to Livestock Production: 2006. Brazil: Belo Horizonte; 2006. p. 1156–65.
Redon R, Ishikawa S, Fitch KR, Feuk L, Perry GH, Andrews TD, Fiegler H, Shapero MH, Carson AR, Chen W, et al. Global variation in copy number in the human genome. Nature. 2006;444(7118):444–54.
Kadri NK, Sahana G, Charlier C, Iso-Touru T, Guldbrandtsen B, Karim L, Nielsen US, Panitz F, Aamand GP, Schulman N, et al. A 660-kb deletion with antagonistic effects on fertility and milk production segregates at high frequency in nordic red cattle: additional evidence for the common occurrence of balancing selection in livestock. PLoS Genet. 2014;10(1):e1004049.
Hirano T, Kobayashi N, Itoh T, Takasuga A, Nakamaru T, Hirotsune S, Sugimoto Y. Null mutation of PCLN-1/Claudin-16 results in bovine chronic interstitial nephritis. Genome Res. 2000;10(5):659–63.
Drogemuller C, Distl O, Leeb T. Partial deletion of the bovine ED1 gene causes anhidrotic ectodermal dysplasia in cattle. Genome Res. 2001;11(10):1699–705.
Sugimoto M, Furuoka H, Sugimoto Y. Deletion of one of the duplicated Hsp70 genes causes hereditary myopathy of diaphragmatic muscles in Holstein-Friesian cattle. Anim Genet. 2003;34(3):191–7.
Meyers SN, McDaneld TG, Swist SL, Marron BM, Steffen DJ, O'Toole D, O'Connell JR, Beever JE, Sonstegard TS, Smith TP. A deletion mutation in bovine SLC4A2 is associated with osteopetrosis in red Angus cattle. BMC Genomics. 2010;11:337.
Zhou X, Stephens M. Genome-wide efficient mixed-model analysis for association studies. Nat Genet. 2012;44(7):821–4.
Thielens NM, Bersch B, Hernandez JF, Arlaud GJ. Structure and functions of the interaction domains of C1r and C1s: keystones of the architecture of the C1 complex. Immunopharmacology. 1999;42(1–3):3–13.
National Center for Biotechnology Information: RNA-seq exon coverage, aggregate (filtered), NCBI Bos taurus Annotation Release 106. https://www.ncbi.nlm.nih.gov/gene/?term=XM_024992686.1 (2020) Accessed 14 Nov 2020.
Fagerberg L, Hallstrom BM, Oksvold P, Kampf C, Djureinovic D, Odeberg J, Habuka M, Tahmasebpoor S, Danielsson A, Edlund K, et al. Analysis of the human tissue-specific expression by genome-wide integration of transcriptomics and antibody-based proteomics. Mol Cell Proteomics. 2014;13(2):397–406.
Lizio M, Abugessaisa I, Noguchi S, Kondo A, Hasegawa A, Hon CC, de Hoon M, Severin J, Oki S, Hayashizaki Y, et al. Update of the FANTOM web resource: expansion to provide additional transcriptome atlases. Nucleic Acids Res. 2019;47(D1):D752–8.
Lizio M, Harshbarger J, Shimoji H, Severin J, Kasukawa T, Sahin S, Abugessaisa I, Fukuda S, Hori F, Ishikawa-Kato S, et al. Gateways to the FANTOM5 promoter level mammalian expression atlas. Genome Biol. 2015;16:22.
Lintner KE, Wu YL, Yang Y, Spencer CH, Hauptmann G, Hebert LA, Atkinson JP, Yu CY. Early components of the complement classical activation pathway in human systemic autoimmune diseases. Front Immunol. 2016;7:36.
Loos M, Heinz HP. Component deficiencies. 1. The first component: C1q, C1r, C1s. Prog Allergy. 1986;39:212–31.
Lee SL, Wallace SL, Barone R, Blum L, Chase PH. Familial deficiency of two subunits of the first component of complement. C1r and C1s associated with a lupus erythematosus-like disease. Arthritis Rheum. 1978;21(8):958–67.
Day NK, Geiger H, Stroud R, DeBracco M, Mancaido B, Windhorst D, Good RA. C1r deficiency: an inborn error associated with cutaneous and renal disease. J Clin Invest. 1972;51(5):1102–8.
Thielens NM, Enrie K, Lacroix M, Jaquinod M, Hernandez JF, Esser AF, Arlaud GJ. The N-terminal CUB-epidermal growth factor module pair of human complement protease C1r binds Ca2+ with high affinity and mediates Ca2+−dependent interaction with C1s. J Biol Chem. 1999;274(14):9149–59.
Tsai SW, Poon PH, Schumaker VN. Expression and characterization of a 159 amino acid, N-terminal fragment of human complement component C1s. Mol Immunol. 1997;34(18):1273–80.
Busby TF, Ingham KC. NH2-terminal calcium-binding domain of human complement C1s- mediates the interaction of C1r- with C1q. Biochemistry. 1990;29(19):4613–8.
Ministry of Agriculture Forestry and Fisheries: Animal Welfare in Japan. https://www.maff.go.jp/j/chikusan/sinko/animal_welfare.html (2020) Accessed 25 July 2020.
Rincon G, Weber KL, Eenennaam AL, Golden BL, Medrano JF. Hot topic: performance of bovine high-density genotyping platforms in Holsteins and jerseys. J Dairy Sci. 2011;94(12):6116–21.
National Center for Biotechnology Information Support Center: Bos taurus assembly ARS-UCD1.2. https://www.ncbi.nlm.nih.gov/assembly/GCF_002263795.1/ (2020) Accessed 14 Apr 2020.
Wang K, Li M, Hadley D, Liu R, Glessner J, Grant SF, Hakonarson H, Bucan M. PennCNV: an integrated hidden Markov model designed for high-resolution copy number variation detection in whole-genome SNP genotyping data. Genome Res. 2007;17(11):1665–74.
Visscher PM, Hemani G, Vinkhuyzen AA, Chen GB, Lee SH, Wray NR, Goddard ME, Yang J. Statistical power to detect genetic (co) variance of complex traits using SNP data in unrelated samples. PLoS Genet. 2014;10(4):e1004269.
Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25(14):1754–60.
Li H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics. 2011;27(21):2987–93.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. Genome project data processing S: the sequence alignment/map format and SAMtools. Bioinformatics. 2009;25(16):2078–9.
PICARD. http://sourceforge.net/projects/picard/ Accessed 7 May 2019.
Genome Analysis Toolkit (GATK). https://www.broadinstitute.org/gatk/ (2019) Accessed 7 May 2019.
Bae JS, Cheong HS, Kim LH, NamGung S, Park TJ, Chun JY, Kim JY, Pasaje CF, Lee JS, Shin HD. Identification of copy number variations and common deletion polymorphisms in cattle. BMC Genomics. 2010;11:232.
Costabile M. Measuring the 50% haemolytic complement (CH50) activity of serum. J Vis Exp. 2010;37.
We would like to thank the members of the WGDB consortium. We would like to thank Kiyoumi Imamura, Takako Arauchi, Kanako Yoshinari, Tomomi Uechi for their technical assistance and other lab members for their generous support and valuable suggestions.
This work was supported by grants from the Project of the NARO Bio-oriented Technology Research Advancement Institution (Research program on development of innovative technology) (27008 BC) (SS, YM, YY, YS1, and YS2) and the Japan Racing and Livestock Promotion (SS, YS2). 1Yoshikazu Sugimoto, 2Yutaka Suzuki.
Ethics approval and consent to participate
All animal experiments were conducted according to the guidelines for the care and use of laboratory animals of the Shirakawa Institute of Animal Genetics, Northern Center of Agricultural Technology, General Technological Center of Hyogo Prefecture for Agriculture, and the University of the Ryukyus. This research was approved by Shirakawa Institute of Animal Genetics Committee on Animal Research (H21–2) and the Committees on Ethics in Animal Experiments in the University of the Ryukyus (A2018128). We obtained written agreement from the cattle owners to use the samples and data.
Consent for publication
The authors declare no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The estimated breeding value (EBV) distribution of 791 cows’ 1–180-day-old calves that died
Features of deleted-type CNVRs on autosomes in this study. Positions are based on the ARS-UCD1.2 assembly of the bovine genome. NA: 12 SNPs in BovineHD beadchip (UMD3.1 assembly) were not assigned by the ARS-UCD1.2 assembly. Table S2 Features of CNVR_221 and the surrounding SNPs. The SNP positions are based on the ARS-UCD1.2 assembly of the bovine genome. Table S3 Detailed features of genes in CNVR_221 and in close proximity to CNVR_221. The SNP positions are based on the ARS-UCD1.2 assembly of the bovine genome. Table S4 Detailed features of the domain of complement C1-related protein in CNVR_221 and in close proximity to CNVR_221. Domains were determined by the SMART protein sequence analysis tool. Table S5 Descriptive statistics of dead calves that suffered from diarrhea or pneumonia. Table S6 Primers and probes for CNVR_221 and BTF3. The SNP positions are based on the ARS-UCD1.2 assembly of the bovine genome. Table S7 Primers and probes for C1R subcomponent-like and C1R.
RNAseq data of the C1RL transcript (LOC112446756, XM_024992686.1) from National Center for Biotechnology Information Bos taurus Annotation release 106. Available from URL:https://www.ncbi.nlm.nih.gov/gene/?term=XM_024992686.1
Principal component analysis of dead calves and healthy cows. The first twenty eigenvectors and their values were calculated using an estimated numerator relationship matrix from 246 dead calves (magenta) and 287 healthy cows (blue). Scatter plots of PC1 vs. PC2 (A) and PC1 vs. PC3 (B).
About this article
Cite this article
Sasaki, S., Miki, Y., Ibi, T. et al. A 44-kb deleted-type copy number variation is associated with decreasing complement component activity and calf mortality in Japanese Black cattle. BMC Genomics 22, 107 (2021). https://doi.org/10.1186/s12864-021-07415-6
- Postnatal mortality
- Copy number variation
- Beef cattle