Genetic Parameters of Semen Traits and Their Correlations with Conformation Traits in Chinese Holstein Bulls.

The elite bull plays an extremely important role in the genetic progression of the dairy cow population. The previous results indicated the potential positive relationship of large scrotal circumference (SC) with improved semen volume, concentration, and motility. In order to improve bull's semen quantity and quality by selection, it is necessary to estimate the genetic parameters of semen traits and their correlations with other conformation traits such as SC that could be used for an indirect selection. In this study, the genetic parameters of seven semen traits (n = 66,260) and nine conformation traits (n = 3,642) of Holstein bulls (n = 453) were estimated by using the bivariate repeatability animal model with the average information-restricted maximum likelihood (AI-REML) approach. The results showed that the estimated heritabilities of semen traits ranged from 0.06 (total number of motile sperm, TNMS) to 0.37 (percentage of abnormal sperm, PAS) and conformation traits ranged from 0.23 (pin width, PW) to 0.69 (hip height, HH). The highest genetic correlations were found between semen volume per ejaculation (SVPE), semen concentration per ejaculation (SCPE), total number of sperm (TNS), and TNMS traits that were 0.97, 0.98, 1.00, and 0.99, respectively. Phenotypic correlations between SC and SVPE, SCPE, TNS, and TNMS were 0.35, 0.35, 0.48, and 0.42, respectively. In summary, the moderate or high heritability of semen traits indicates that genetic improvement of semen quality by selection is feasible, where SC could be a useful trait for indirect selection or as correlated information to improve semen quantity and production in the practical bull breeding programs.

Genetic variation in the human leukocyte antigen region confers susceptibility to Clostridioides difficile infection.

Clostridioides difficile (C. diff.) infection (CDI) is a leading cause of hospital acquired diarrhea in North America and Europe and a major cause of morbidity and mortality. Known risk factors do not fully explain CDI susceptibility, and genetic susceptibility is suggested by the fact that some patients with colons that are colonized with C. diff. do not develop any infection while others develop severe or recurrent infections. To identify common genetic variants associated with CDI, we performed a genome-wide association analysis in 19,861 participants (1349 cases; 18,512 controls) from the Electronic Medical Records and Genomics (eMERGE) Network. Using logistic regression, we found strong evidence for genetic variation in the DRB locus of the MHC (HLA) II region that predisposes individuals to CDI (P > 1.0 × 10-14; OR 1.56). Altered transcriptional regulation in the HLA region may play a role in conferring susceptibility to this opportunistic enteric pathogen.

The single-cell transcriptome and chromatin accessibility datasets of peripheral blood mononuclear cells in Chinese holstein cattle.

OBJECTIVES: This study was performed in the frame of a more extensive study dedicated to the integrated analysis of the single-cell transcriptome and chromatin accessibility datasets of peripheral blood mononuclear cells (PBMCs) with a large-scale GWAS of 45 complex traits in Chinese Holstein cattle. Lipopolysaccharide (LPS) is a crucial mediator of chronic inflammation to modulate immune responses. PBMCs include primary T and B cells, natural killer (NK) cells, monocytes (Mono), and dendritic cells (DC). How LPS stimulates PBMCs at the single-cell level in dairy cattle remains largely unknown. DATA DESCRIPTION: We sequenced 30,756 estimated single cells and mapped 26,141 of them (96.05%) with approximately 60,075 mapped reads per cell after quality control for four whole-blood treatments (no, 2 h, 4 h, and 8 h LPS) by single-cell RNA sequencing (scRNA-seq) and single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq). Finally, 7,107 (no), 9,174 (2 h), 6,741 (4 h), and 3,119 (8 h) cells were generated with ~ 15,000 total genes in the whole population. Therefore, the single-cell transcriptome and chromatin accessibility datasets in this study enable a further understanding of the cell types and functions of PBMCs and their responses to LPS stimulation in vitro.

Preliminary Study on Gene Regulation and its Pathways in Chinese Holstein Cows with Clinical Mastitis Caused by Staphylococcus Aureus.

Introduction: Clinical mastitis (CM) is one of the most common diseases of dairy cows globally, has a complex aetiology and recurs easily. Staphylococcus aureus is a frequently isolated pathogen responsible for bovine mastitis and remains difficult to eradicate. Material and Methods: To characterise the transcriptional profiles of dairy cows infected by S. aureus, we performed an RNA-seq analysis of peripheral blood leukocytes in lactating Chinese Holstein dairy cows with CM and did the same with healthy cows' samples as controls. Results: A total of 4,286 genes were detected in the CM cases infected with S. aureus which were differentially expressed compared to the controls, 3,085 of which were upregulated, the remainder being downregulated. Notably, we observed that some differentially expressed genes (DEGs) had strong protein-protein interaction. Of these, six downregulated DEGs (AKR1C4, PTGS2, HNMT, EPHX2, CMBL, and IDH1) were involved in the metabolic pathway, while eight upregulated DEGs (VWF, GP9, MYLK, GP6, F2RL3, ITGB3, GP5, and PRKG1) were associated with the platelet activation pathway. Conclusion: The transcriptome dataset of CM cases would be a valuable resource for clinical guidance on anti-inflammatory medication and for deeper understanding of the biological processes of CM response to S. aureus infection, and it would enable us to identify specific genes for diagnostic markers and possibly for targeted therapy.

Single-cell transcriptomic and chromatin accessibility analyses of dairy cattle peripheral blood mononuclear cells and their responses to lipopolysaccharide.

BACKGROUND: Gram-negative bacteria are important pathogens in cattle, causing severe infectious diseases, including mastitis. Lipopolysaccharides (LPS) are components of the outer membrane of Gram-negative bacteria and crucial mediators of chronic inflammation in cattle. LPS modulations of bovine immune responses have been studied before. However, the single-cell transcriptomic and chromatin accessibility analyses of bovine peripheral blood mononuclear cells (PBMCs) and their responses to LPS stimulation were never reported. RESULTS: We performed single-cell RNA sequencing (scRNA-seq) and single-cell sequencing assay for transposase-accessible chromatin (scATAC-seq) in bovine PBMCs before and after LPS treatment and demonstrated that seven major cell types, which included CD4 T cells, CD8 T cells, and B cells, monocytes, natural killer cells, innate lymphoid cells, and dendritic cells. Bioinformatic analyses indicated that LPS could increase PBMC cell cycle progression, cellular differentiation, and chromatin accessibility. Gene analyses further showed significant changes in differential expression, transcription factor binding site, gene ontology, and regulatory interactions during the PBMC responses to LPS. Consistent with the findings of previous studies, LPS induced activation of monocytes and dendritic cells, likely through their upregulated TLR4 receptor. NF-κB was observed to be activated by LPS and an increased transcription of an array of pro-inflammatory cytokines, in agreement that NF-κB is an LPS-responsive regulator of innate immune responses. In addition, by integrating LPS-induced differentially expressed genes (DEGs) with large-scale GWAS of 45 complex traits in Holstein, we detected trait-relevant cell types. We found that selected DEGs were significantly associated with immune-relevant health, milk production, and body conformation traits. CONCLUSION: This study provided the first scRNAseq and scATAC-seq data for cattle PBMCs and their responses to the LPS stimulation to the best of our knowledge. These results should also serve as valuable resources for the future study of the bovine immune system and open the door for discoveries about immune cell roles in complex traits like mastitis at single-cell resolution.

The clinical utility of polygenic risk scores in genomic medicine practices: a systematic review.

Genomic medicine aims to improve health using the individual genomic data of people to inform care. While clinical utility of genomic medicine in many monogenic, Mendelian disorders is amply demonstrated, clinical utility is less evident in polygenic traits, e.g., coronary artery disease or breast cancer. Polygenic risk scores (PRS) are subsets of individual genotypes designed to capture heritability of common traits, and hence to allow the stratification of risk of the trait in a population. We systematically reviewed the PubMed database for unequivocal evidence of clinical utility of polygenic risk scores, using stringent inclusion and exclusion criteria. While we identified studies demonstrating clinical validity in conditions where medical intervention based on a PRS is likely to benefit patient outcome, we did not identify a single study demonstrating unequivocally such a benefit, i.e. clinical utility. We conclude that while the routine use of PRSs hold great promise, translational research is still needed before they should enter mainstream clinical practice.

Penetrance of Breast Cancer Susceptibility Genes From the eMERGE III Network.

Background: Unbiased estimates of penetrance are challenging but critically important to make informed choices about strategies for risk management through increased surveillance and risk-reducing interventions. Methods: We studied the penetrance and clinical outcomes of 7 breast cancer susceptibility genes (BRCA1, BRCA2, TP53, CHEK2, ATM, PALB2, and PTEN) in almost 13 458 participants unselected for personal or family history of breast cancer. We identified 242 female participants with pathogenic or likely pathogenic variants in 1 of the 7 genes for penetrance analyses, and 147 women did not previously know their genetic results. Results: Out of the 147 women, 32 women were diagnosed with breast cancer at an average age of 52.8 years. Estimated penetrance by age 60 years ranged from 17.8% to 43.8%, depending on the gene. In clinical-impact analysis, 42.3% (95% confidence interval = 31.3% to 53.3%) of women had taken actions related to their genetic results, and 2 new breast cancer cases were identified within the first 12 months after genetic results disclosure. Conclusions: Our study provides population-based penetrance estimates for the understudied genes CHEK2, ATM, and PALB2 and highlights the importance of using unselected populations for penetrance studies. It also demonstrates the potential clinical impact of genetic testing to improve health care through early diagnosis and preventative screening.

Loci identified by a genome-wide association study of carotid artery stenosis in the eMERGE network.

Carotid artery atherosclerotic disease (CAAD) is a risk factor for stroke. We used a genome-wide association (GWAS) approach to discover genetic variants associated with CAAD in participants in the electronic Medical Records and Genomics (eMERGE) Network. We identified adult CAAD cases with unilateral or bilateral carotid artery stenosis and controls without evidence of stenosis from electronic health records at eight eMERGE sites. We performed GWAS with a model adjusting for age, sex, study site, and genetic principal components of ancestry. In eMERGE we found 1793 CAAD cases and 17,958 controls. Two loci reached genome-wide significance, on chr6 in LPA (rs10455872, odds ratio [OR] (95% confidence interval [CI]) = 1.50 (1.30-1.73), p = 2.1 × 10-8 ) and on chr7, an intergenic single nucleotide variant (SNV; rs6952610, OR (95% CI) = 1.25 (1.16-1.36), p = 4.3 × 10-8 ). The chr7 association remained significant in the presence of the LPA SNV as a covariate. The LPA SNV was also associated with coronary heart disease (CHD; 4199 cases and 11,679 controls) in this study (OR (95% CI) = 1.27 (1.13-1.43), p = 5 × 10-5 ) but the chr7 SNV was not (OR (95% CI) = 1.03 (0.97-1.09), p = .37). Both variants replicated in UK Biobank. Elevated lipoprotein(a) concentrations ([Lp(a)]) and LPA variants associated with elevated [Lp(a)] have previously been associated with CAAD and CHD, including rs10455872. With electronic health record phenotypes in eMERGE and UKB, we replicated a previously known association and identified a novel locus associated with CAAD.

Donkey genomes provide new insights into domestication and selection for coat color.

Current knowledge about the evolutionary history of donkeys is still incomplete due to the lack of archeological and whole-genome diversity data. To fill this gap, we have de novo assembled a chromosome-level reference genome of one male Dezhou donkey and analyzed the genomes of 126 domestic donkeys and seven wild asses. Population genomics analyses indicate that donkeys were domesticated in Africa and conclusively show reduced levels of Y chromosome variability and discordant paternal and maternal histories, possibly reflecting the consequences of reproductive management. We also investigate the genetic basis of coat color. While wild asses show diluted gray pigmentation (Dun phenotype), domestic donkeys display non-diluted black or chestnut coat colors (non-Dun) that were probably established during domestication. Here, we show that the non-Dun phenotype is caused by a 1 bp deletion downstream of the TBX3 gene, which decreases the expression of this gene and its inhibitory effect on pigment deposition.

Impact of the Order of Legendre Polynomials in Random Regression Model on Genetic Evaluation for Milk Yield in Dairy Cattle Population.

The random regression test-day model has become the most commonly adopted model for routine genetic evaluations in dairy populations, which allows accurately accounting for genetic and environmental effects over lactation. The objective of this study was to explore appropriate random regression test-day models for genetic evaluation of milk yield in a Holstein population with a relatively small size, which is the common situation in regional or independent breeding companies to preform genetic evaluation. Data included 419,567 test-day records from 54,417 cows from the first lactation. Variance components and breeding values were estimated using a random regression test-day model with different orders (from first to fifth) of Legendre polynomials (LP) and accounted for homogeneous or heterogeneous residual variance across the lactation. Models were compared based on Akaike information criterion (AIC), Bayesian information criterion (BIC), and predictive ability. In general, models with a higher order of LP showed better goodness of fit based on AIC and BIC values. However, models with third, fourth, and fifth order of LP led to similar estimates of genetic parameters and predictive ability. Models with assumption of heterogeneous residual variances achieved better goodness of fit but yielded similar predictive ability compared with those with assumption of homogeneous residual variances. Therefore, a random regression model with third order of LP is suggested to be an appropriate model for genetic evaluation of milk yield in local Chinese Holstein populations.

Identification of a dynamic system model for a building and heating system including heat pump and thermal energy storage.

Controllers employing optimal control strategies will path the way to enable flexible operations in future power grids. As buildings will increasingly act as prosumers in future power grids, optimal control of buildings' energy consumption will play a major role in providing flexible operations. Optimal controllers such as model predictive controller are able to manage buildings' operations and to optimise their energy consumption. For online optimisation, model predictive controller requires a model of the energy system. The more accurate the system model represents the system dynamics, the more accurate the model predictive controller predicts the future states of the energy system while optimising its energy consumption. In this article, we present a system model that can be used in online MPC, including dynamic programming as optimisation strategy. The system model is validated using a building and heating system, including heat pump and thermal energy storage. The following bullet points summarise the main requirements for the configuration of the system model:•The system model performs fast with low computational effort in less than 1 s;•The system model can be implemented in online MPC;•The system model accurately represents the dynamic behaviour.

Machine learning algorithms applied to a prediction of personal overall thermal comfort using skin temperatures and occupants' heating behavior.

Thermal comfort modeling has been of interest in built environment research for decades. Mostly the modeling approaches focused on an average response of a large group of building occupants. Recently, the focus has been shifted towards personal comfort models that predict individuals' thermal comfort responses. Currently, thermal comfort responses are collected from the occupants via survey. This study explored if the thermal comfort of individuals could be predicted using machine learning algorithms while relaying on the set of collected inputs from an experiment. The model was developed using experimental data including collected from a previously performed experiment in the climate chamber. Two different approaches based on the output data (thermal sensation and thermal comfort votes) and five different sets of input variables were explored. The algorithms tested were Support Vector Machine with four different Kernel functions (Linear, Quadratic, Cubic and Gaussian) and Ensemble Algorithms (Boosted trees, Bagged trees and RUSBoosted trees). The combination of occupants' heating behavior with a personal comfort system (PCS), skin temperatures, time and environmental data were used for the development of personal comfort models to predict individuals' thermal preference. The study investigated the novel combination of inputs such as the use of skin temperature and settings of the personalized heating system as parameters in predicting personal thermal comfort. The results showed that personal comfort models among all tested approaches and subjects showed the best median accuracy of 0.84 using RUSBoosted trees. Individually looking, the approach using thermal sensation output produced better prediction accuracy. On the other hand, the models based on inputs that consisted of PCS control behavior and mean and hand skin temperatures produced the best prediction accuracy when assessing all tested algorithms. The main limitation of the study is the number of test subjects, and further recommendation is to perform more experiments.

GWAS and enrichment analyses of non-alcoholic fatty liver disease identify new trait-associated genes and pathways across eMERGE Network.

BACKGROUND: Non-alcoholic fatty liver disease (NAFLD) is a common chronic liver illness with a genetically heterogeneous background that can be accompanied by considerable morbidity and attendant health care costs. The pathogenesis and progression of NAFLD is complex with many unanswered questions. We conducted genome-wide association studies (GWASs) using both adult and pediatric participants from the Electronic Medical Records and Genomics (eMERGE) Network to identify novel genetic contributors to this condition. METHODS: First, a natural language processing (NLP) algorithm was developed, tested, and deployed at each site to identify 1106 NAFLD cases and 8571 controls and histological data from liver tissue in 235 available participants. These include 1242 pediatric participants (396 cases, 846 controls). The algorithm included billing codes, text queries, laboratory values, and medication records. Next, GWASs were performed on NAFLD cases and controls and case-only analyses using histologic scores and liver function tests adjusting for age, sex, site, ancestry, PC, and body mass index (BMI). RESULTS: Consistent with previous results, a robust association was detected for the PNPLA3 gene cluster in participants with European ancestry. At the PNPLA3-SAMM50 region, three SNPs, rs738409, rs738408, and rs3747207, showed strongest association (best SNP rs738409 p = 1.70 × 10- 20). This effect was consistent in both pediatric (p = 9.92 × 10- 6) and adult (p = 9.73 × 10- 15) cohorts. Additionally, this variant was also associated with disease severity and NAFLD Activity Score (NAS) (p = 3.94 × 10- 8, beta = 0.85). PheWAS analysis link this locus to a spectrum of liver diseases beyond NAFLD with a novel negative correlation with gout (p = 1.09 × 10- 4). We also identified novel loci for NAFLD disease severity, including one novel locus for NAS score near IL17RA (rs5748926, p = 3.80 × 10- 8), and another near ZFP90-CDH1 for fibrosis (rs698718, p = 2.74 × 10- 11). Post-GWAS and gene-based analyses identified more than 300 genes that were used for functional and pathway enrichment analyses. CONCLUSIONS: In summary, this study demonstrates clear confirmation of a previously described NAFLD risk locus and several novel associations. Further collaborative studies including an ethnically diverse population with well-characterized liver histologic features of NAFLD are needed to further validate the novel findings.

In silico genome-wide miRNA-QTL-SNPs analyses identify a functional SNP associated with mastitis in Holsteins.

BACKGROUND: Single-nucleotide polymorphisms (SNPs) in microRNAs (miRNAs) and their target binding sites affect miRNA function and are involved in biological processes and diseases, including bovine mastitis, a frequent inflammatory disease. Our previous study has shown that bta-miR-2899 is significantly upregulated in the mammary gland tissue of mastitis-infected cow than that of healthy cows. RESULTS: In the present study, we used a customized miRNAQTLsnp software and identified 5252 SNPs in 691 bovine pre-miRNAs, which are also located within the quantitative trait loci (QTLs) that are associated with mastitis and udder conformation-related traits. Using luciferase assay in the bovine mammary epithelial cells, we confirmed a candidate SNP (rs109462250, g. 42,198,087 G > A) in the seed region of bta-miR-2899 located in the somatic cell score (SCS)-related QTL (Chr.18: 33.9-43.9 Mbp), which affected the interaction of bta-miR-2899 and its putative target Spi-1 proto-oncogene (SPI1), a pivotal regulator in the innate and adaptive immune systems. Quantitative real-time polymerase chain reaction results showed that the relative expression of SPI1 in the mammary gland of AA genotype cows was significantly higher than that of GG genotype cows. The SNP genotypes were associated with SCS in Holstein cows. CONCLUSIONS: Altogether, miRNA-related SNPs, which influence the susceptibility to mastitis, are one of the plausible mechanisms underlying mastitis via modulating the interaction of miRNAs and immune-related genes. These miRNA-QTL-SNPs, such as the SNP (rs109462250) of bta-miR-2899 may have implication for the mastitis resistance breeding program in Holstein cattle.

The eMERGE genotype set of 83,717 subjects imputed to ~40 million variants genome wide and association with the herpes zoster medical record phenotype.

The Electronic Medical Records and Genomics (eMERGE) network is a network of medical centers with electronic medical records linked to existing biorepository samples for genomic discovery and genomic medicine research. The network sought to unify the genetic results from 78 Illumina and Affymetrix genotype array batches from 12 contributing medical centers for joint association analysis of 83,717 human participants. In this report, we describe the imputation of eMERGE results and methods to create the unified imputed merged set of genome-wide variant genotype data. We imputed the data using the Michigan Imputation Server, which provides a missing single-nucleotide variant genotype imputation service using the minimac3 imputation algorithm with the Haplotype Reference Consortium genotype reference set. We describe the quality control and filtering steps used in the generation of this data set and suggest generalizable quality thresholds for imputation and phenotype association studies. To test the merged imputed genotype set, we replicated a previously reported chromosome 6 HLA-B herpes zoster (shingles) association and discovered a novel zoster-associated loci in an epigenetic binding site near the terminus of chromosome 3 (3p29).

Enrichment sampling for a multi-site patient survey using electronic health records and census data.

Objective: We describe a stratified sampling design that combines electronic health records (EHRs) and United States Census (USC) data to construct the sampling frame and an algorithm to enrich the sample with individuals belonging to rarer strata. Materials and Methods: This design was developed for a multi-site survey that sought to examine patient concerns about and barriers to participating in research studies, especially among under-studied populations (eg, minorities, low educational attainment). We defined sampling strata by cross-tabulating several socio-demographic variables obtained from EHR and augmented with census-block-level USC data. We oversampled rarer and historically underrepresented subpopulations. Results: The sampling strategy, which included USC-supplemented EHR data, led to a far more diverse sample than would have been expected under random sampling (eg, 3-, 8-, 7-, and 12-fold increase in African Americans, Asians, Hispanics and those with less than a high school degree, respectively). We observed that our EHR data tended to misclassify minority races more often than majority races, and that non-majority races, Latino ethnicity, younger adult age, lower education, and urban/suburban living were each associated with lower response rates to the mailed surveys. Discussion: We observed substantial enrichment from rarer subpopulations. The magnitude of the enrichment depends on the accuracy of the variables that define the sampling strata and the overall response rate. Conclusion: EHR and USC data may be used to define sampling strata that in turn may be used to enrich the final study sample. This design may be of particular interest for studies of rarer and understudied populations.

Parents' attitudes toward consent and data sharing in biobanks: A multisite experimental survey.

BACKGROUND: The factors influencing parents' willingness to enroll their children in biobanks are poorly understood. This study sought to assess parents' willingness to enroll their children, and their perceived benefits, concerns, and information needs under different consent and data-sharing scenarios, and to identify factors associated with willingness. METHODS: This large, experimental survey of patients at the 11 eMERGE Network sites used a disproportionate stratified sampling scheme to enrich the sample with historically underrepresented groups. Participants were randomized to receive one of three consent and data-sharing scenarios. RESULTS: In total, 90,000 surveys were mailed and 13,000 individuals responded (15.8% response rate). 5737 respondents were parents of minor children. Overall, 55% (95% confidence interval 50-59%) of parents were willing to enroll their youngest minor child in a hypothetical biobank; willingness did not differ between consent and data-sharing scenarios. Lower educational attainment, higher religiosity, lower trust, worries about privacy, and attitudes about benefits, concerns, and information needs were independently associated with less willingness to allow their child to participate. Of parents who were willing to participate themselves, 25% were not willing to allow their child to participate. Being willing to participate but not willing to allow one's child to participate was independently associated with multiple factors, including race, lower educational attainment, lower annual household income, public health care insurance, and higher religiosity. CONCLUSIONS: Fifty-five percent of parents were willing to allow their youngest minor child to participate in a hypothetical biobank. Building trust, protecting privacy, and addressing attitudes may increase enrollment and diversity in pediatric biobanks.

The effect of the SNP g.18475 A>G in the 3'UTR of NCF4 on mastitis susceptibility in dairy cattle.

Neutrophil cytosolic factor 4 (NCF4) is a member of the nicotinamide adenine dinucleotide phosphate oxidase subunit. This protein functions as an essential factor in the host defense against the progression of bacterial infection. To explore the variability of the NCF4 gene and the susceptibility of cows to mastitis, NCF4 functional single nucleotide polymorphism (SNP) of the 3' untranslated region (3'UTR) and its targeted microRNA (miRNA) were identified. One SNP g.18475 A>G in the 3'UTR of NCF4 was found within the binding seed region of bta-miR-2426. We constructed two recombinant pMIR-REPORT™ vectors with the A or G allele in the g.18475 locus and transiently co-transfected the vectors in human embryo kidney 293T (HEK 293T) cells, along with bta-miR-2426 mimics. A luciferase assay indicated that this SNP affects the binding of NCF4 and bta-miR-2426. In addition, the association analysis results showed that cows with the GG genotype in SNP g.18475 A>G had a relatively lower SCS value than cows with the AA genotype. Finally, quantitative real-time PCR (RT-qPCR) results showed that the cows with genotype GG had a relatively higher expression of NCF4 mRNA compared to the cows with genotype AA. NCF4 expression was regulated by the miRNA-mRNA interaction mechanism, and an important role for NCF4 in mastitis susceptibility in dairy cow was suggested.