Integrating Sequence-based GWAS and RNA-Seq Provides Novel Insights into the Genetic Basis of Mastitis and Milk Production in Dairy Cattle.
Ontology highlight
ABSTRACT:
SUBMITTER: Fang L
PROVIDER: S-EPMC5372096 | biostudies-literature | 2017 Mar
REPOSITORIES: biostudies-literature
Publications
Integrating Sequence-based GWAS and RNA-Seq Provides Novel Insights into the Genetic Basis of Mastitis and Milk Production in Dairy Cattle.
Fang Lingzhao L Sahana Goutam G Su Guosheng G Yu Ying Y Zhang Shengli S Lund Mogens Sandø MS Sørensen Peter P
Scientific reports 20170330
Connecting genome-wide association study (GWAS) to biological mechanisms underlying complex traits is a major challenge. Mastitis resistance and milk production are complex traits of economic importance in the dairy sector and are associated with intra-mammary infection (IMI). Here, we integrated IMI-relevant RNA-Seq data from Holstein cattle and sequence-based GWAS data from three dairy cattle breeds (i.e., Holstein, Nordic red cattle, and Jersey) to explore the genetic basis of mastitis resist ...[more]
Similar Datasets
Project description:BackgroundKetosis is a common metabolic disease during the transition period in dairy cattle, resulting in long-term economic loss to the dairy industry worldwide. While genetic selection of resistance to ketosis has been adopted by many countries, the genetic and biological basis underlying ketosis is poorly understood.ResultsWe collected a total of 24 blood samples from 12 Holstein cows, including 4 healthy and 8 ketosis-diagnosed ones, before (2 weeks) and after (5 days) calving, respectively. We then generated RNA-Sequencing (RNA-Seq) data and seven blood biochemical indicators (bio-indicators) from leukocytes and plasma in each of these samples, respectively. By employing a weighted gene co-expression network analysis (WGCNA), we detected that 4 out of 16 gene-modules, which were significantly engaged in lipid metabolism and immune responses, were transcriptionally (FDR < 0.05) correlated with postpartum ketosis and several bio-indicators (e.g., high-density lipoprotein and low-density lipoprotein). By conducting genome-wide association signal (GWAS) enrichment analysis among six common health traits (ketosis, mastitis, displaced abomasum, metritis, hypocalcemia and livability), we found that 4 out of 16 modules were genetically (FDR < 0.05) associated with ketosis, among which three were correlated with postpartum ketosis based on WGCNA. We further identified five candidate genes for ketosis, including GRINA, MAF1, MAFA, C14H8orf82 and RECQL4. Our phenome-wide association analysis (Phe-WAS) demonstrated that human orthologues of these candidate genes were also significantly associated with many metabolic, endocrine, and immune traits in humans. For instance, MAFA, which is involved in insulin secretion, glucose response, and transcriptional regulation, showed a significantly higher association with metabolic and endocrine traits compared to other types of traits in humans.ConclusionsIn summary, our study provides novel insights into the molecular mechanism underlying ketosis in cattle, and highlights that an integrative analysis of omics data and cross-species mapping are promising for illustrating the genetic architecture underpinning complex traits.
| S-EPMC7367229 | biostudies-literature
Project description:BackgroundSubclinical intramammary infection (IMI) represents a significant problem in maintaining dairy cows' health. Disease severity and extent depend on the interaction between the causative agent, environment, and host. To investigate the molecular mechanisms behind the host immune response, we used RNA-Seq for the milk somatic cells (SC) transcriptome profiling in healthy cows (n = 9), and cows naturally affected by subclinical IMI from Prototheca spp. (n = 11) and Streptococcus agalactiae (S. agalactiae; n = 11). Data Integration Analysis for Biomarker discovery using Latent Components (DIABLO) was used to integrate transcriptomic data and host phenotypic traits related to milk composition, SC composition, and udder health to identify hub variables for subclinical IMI detection.ResultsA total of 1,682 and 2,427 differentially expressed genes (DEGs) were identified when comparing Prototheca spp. and S. agalactiae to healthy animals, respectively. Pathogen-specific pathway analyses evidenced that Prototheca's infection upregulated antigen processing and lymphocyte proliferation pathways while S. agalactiae induced a reduction of energy-related pathways like the tricarboxylic acid cycle, and carbohydrate and lipid metabolism. The integrative analysis of commonly shared DEGs between the two pathogens (n = 681) referred to the core-mastitis response genes, and phenotypic data evidenced a strong covariation between those genes and the flow cytometry immune cells (r2 = 0.72), followed by the udder health (r2 = 0.64) and milk quality parameters (r2 = 0.64). Variables with r ≥ 0.90 were used to build a network in which the top 20 hub variables were identified with the Cytoscape cytohubba plug-in. The genes in common between DIABLO and cytohubba (n = 10) were submitted to a ROC analysis which showed they had excellent predictive performances in terms of discriminating healthy and mastitis-affected animals (sensitivity > 0.89, specificity > 0.81, accuracy > 0.87, and precision > 0.69). Among these genes, CIITA could play a key role in regulating the animals' response to subclinical IMI.ConclusionsDespite some differences in the enriched pathways, the two mastitis-causing pathogens seemed to induce a shared host immune-transcriptomic response. The hub variables identified with the integrative approach might be included in screening and diagnostic tools for subclinical IMI detection.
| S-EPMC10320993 | biostudies-literature
Project description:BackgroundImproving resistance to mastitis, one of the costliest diseases in dairy production, has become an important objective in dairy cattle breeding. However, mastitis resistance is influenced by many genes involved in multiple processes, including the response to infection, inflammation, and post-infection healing. Low genetic heritability, environmental variations, and farm management differences further complicate the identification of links between genetic variants and mastitis resistance. Consequently, studies of the genetics of variation in mastitis resistance in dairy cattle lack agreement about the responsible genes.ResultsWe associated 15,552,968 imputed whole-genome sequencing markers for 5147 Nordic Holstein cattle with mastitis resistance in a genome-wide association study (GWAS). Next, we augmented P-values for markers in genes in the associated regions using Gene Ontology terms, Kyoto Encyclopedia of Genes and Genomes pathway analysis, and mammalian phenotype database. To confirm results of gene-based analyses, we used gene expression data from E. coli-challenged cow udders. We identified 22 independent quantitative trait loci (QTL) that collectively explained 14% of the variance in breeding values for resistance to clinical mastitis (CM). Using association test statistics with multiple pieces of independent information on gene function and differential expression during bacterial infection, we suggested putative causal genes with biological relevance for 12 QTL affecting resistance to CM in dairy cattle.ConclusionCombining information on the nearest positional genes, gene-based analyses, and differential gene expression data from RNA-seq, we identified putative causal genes (candidate genes with biological evidence) in QTL for mastitis resistance in Nordic Holstein cattle. The same strategy can be applied for other traits.
| S-EPMC6127918 | biostudies-literature
Project description:In recent years, studies on the biological mechanisms underlying complex traits have been facilitated by innovations in high-throughput genotyping technology. We conducted a weighted single-step genome-wide association study (WssGWAS) to evaluate backfat thickness, carcass weight, eye muscle area, marbling score, and yearling weight in a cohort of 1540 Hanwoo beef cattle using BovineSNP50 BeadChip. The WssGWAS uncovered thirty-three genomic regions that explained more than 1% of the additive genetic variance, mostly located on chromosomes 6 and 14. Among the identified window regions, seven quantitative trait loci (QTL) had pleiotropic effects and twenty-six QTL were trait-specific. Significant pathways implicated in the measured traits through Gene Ontology (GO) term enrichment analysis included the following: lipid biosynthetic process, regulation of lipid metabolic process, transport or localization of lipid, regulation of growth, developmental growth, and multicellular organism growth. Integration of GWAS results of the studied traits with pathway and network analyses facilitated the exploration of the respective candidate genes involved in several biological functions, particularly lipid and growth metabolism. This study provides novel insight into the genetic bases underlying complex traits and could be useful in developing breeding schemes aimed at improving growth and carcass traits in Hanwoo beef cattle.
| S-EPMC7601430 | biostudies-literature
Project description:BackgroundUrinary nitrogen leakage is an environmental concern in dairy cattle. Selection for reduced urinary nitrogen leakage may be done using indicator traits such as milk urea nitrogen (MUN). The result of a previous study indicated that the genetic correlation between MUN in Australia (AUS) and MUN in New Zealand (NZL) was only low to moderate (between 0.14 and 0.58). In this context, an alternative is to select sequence variants based on genome-wide association studies (GWAS) with a view to improve genomic prediction accuracies. A GWAS can also be used to detect quantitative trait loci (QTL) associated with MUN. Therefore, our objectives were to perform within-country GWAS and a meta-GWAS for MUN using records from up to 33,873 dairy cows and imputed whole-genome sequence data, to compare QTL detected in the GWAS for MUN in AUS and NZL, and to use sequence variants selected from the meta-GWAS to improve the prediction accuracy for MUN based on a joint AUS-NZL reference set.ResultsUsing the meta-GWAS, we detected 14 QTL for MUN, located on chromosomes 1, 6, 11, 14, 19, 22, 26 and the X chromosome. The three most significant QTL encompassed the casein genes on chromosome 6, PAEP on chromosome 11 and DGAT1 on chromosome 14. We selected 50,000 sequence variants that had the same direction of effect for MUN in AUS and MUN in NZL and that were most significant in the meta-analysis for the GWAS. The selected sequence variants yielded a genetic correlation between MUN in AUS and MUN in NZL of 0.95 and substantially increased prediction accuracy in both countries.ConclusionsOur results demonstrate how the sharing of data between two countries can increase the power of a GWAS and increase the accuracy of genomic prediction using a multi-country reference population and sequence variants selected based on a meta-GWAS.
| S-EPMC8858489 | biostudies-literature
Project description:Bovine mastitis continues to be a complex disease associated with significant economic loss in dairy industries worldwide. The incidence rate of subclinical mastitis (IRSCM) can show substantial variation among different farms; however, the milk microbiota, which have a direct influence on bovine mammary gland health, have never been associated with the IRSCM. Here, we aimed to use high-throughput DNA sequencing to describe the milk microbiota from two dairy farms with different IRSCMs and to identify the predominant mastitis pathogens along with commensal or potential beneficial bacteria. Our study showed that Klebsiella, Escherichia-Shigella, and Streptococcus were the mastitis-causing pathogens in farm A (with a lower IRSCM), while Streptococcus and Corynebacterium were the mastitis-causing pathogens in farm B (with a higher IRSCM). The relative abundance of all pathogens in farm B (22.12%) was higher than that in farm A (9.82%). However, the genus Bacillus was more prevalent in farm A. These results may be helpful for explaining the lower IRSCM in farm A. Additionally, the gut-associated genera Prevotella, Ruminococcus, Bacteroides, Rikenella, and Alistipes were prevalent in all milk samples, suggesting gut bacteria can be one of the predominant microbial contamination in milk. Moreover, Listeria monocytogenes (a foodborne pathogen) was found to be prevalent in farm A, even though it had a lower IRSCM. Overall, our study showed complex diversity between the milk microbiota in dairy farms with different IRSCMs. This suggests that variation in IRSCMs may not only be determined by the heterogeneity and prevalence of mastitis-causing pathogens but also be associated with potential beneficial bacteria. In the future, milk microbiota should be considered in bovine mammary gland health management. This would be helpful for both the establishment of a targeted mastitis control system and the control of the safety and quality of dairy products.
| S-EPMC6232673 | biostudies-literature
Project description:BackgroundStaphylococcus aureus is a common cause of clinical mastitis (CM) in dairy cattle. Optimizing the bovine mammary gland microbiota to resist S. aureus colonization is a growing area of research. However, the details of the interbacterial interactions between S. aureus and commensal bacteria, which would be required to manipulate the microbiome to resist infection, are still unknown. This study aims to characterize changes in the bovine milk bacterial community before, during, and after S. aureus CM and to compare bacterial communities present in milk between infected and healthy quarters.MethodsWe collected quarter-level milk samples from 698 Holstein dairy cows over an entire lactation. A total of 11 quarters from 10 cows were affected by S. aureus CM and milk samples from these 10 cows (n = 583) regardless of health status were analyzed by performing 16S rRNA gene amplicon sequencing.ResultsThe milk microbiota of healthy quarters was distinguishable from that of S. aureus CM quarters two weeks before CM diagnosis via visual inspection. Microbial network analysis showed that 11 OTUs had negative associations with OTU0001 (Staphylococcus). A low diversity or dysbiotic milk microbiome did not necessarily correlate with increased inflammation. Specifically, Staphylococcus xylosus, Staphylococcus epidermidis, and Aerococcus urinaeequi were each abundant in milk from the quarters with low levels of inflammation.ConclusionOur results show that the udder microbiome is highly dynamic, yet a change in the abundance in certain bacteria can be a potential indicator of future S. aureus CM. This study has identified potential prophylactic bacterial species that could act as a barrier against S. aureus colonization and prevent future instances of S. aureus CM.
| S-EPMC9701008 | biostudies-literature
Project description:Mastitis is one of the most impacting diseases in dairy farming, and its sensitive and specific detection is therefore of the greatest importance. The clinical evaluation of udder and mammary secretions is typically combined with the milk Somatic Cell Count (SCC) and often accompanied by its bacteriological culture to identify the causative microorganism. In a constant search for improvement, several non-enzymatic milk proteins, including milk amyloid A (M-SAA), haptoglobin (HP), cathelicidin (CATH), and lactoferrin (LF), have been investigated as alternative biomarkers of mastitis for their relationship with mammary gland inflammation, and immunoassay techniques have been developed for detection with varying degrees of success. To provide a general overview of their implementation in the different dairy species, we carried out a systematic review of the scientific literature using the PRISMA (Preferred Reporting Items for Systematic Reviews and Meta-Analysis) guidelines. Our review question falls within the type "Diagnostic test accuracy questions" and aims at answering the diagnostic question: "Which are the diagnostic performances of mastitis protein biomarkers investigated by immunoassays in ruminant milk?". Based on 13 keywords combined into 42 searches, 523 manuscripts were extracted from three scientific databases. Of these, 33 passed the duplicate removal, title, abstract, and full-text screening for conformity to the review question and document type: 78.8% investigated cows, 12.1% sheep, 9.1% goats, and 6.1% buffaloes (some included more than one dairy species). The most frequently mentioned protein was M-SAA (48.5%), followed by HP (27.3%), CATH (24.2%) and LF (21.2%). However, the large amount of heterogeneity among studies in terms of animal selection criteria (45.5%), index test (87.9%), and standard reference test (27.3%) resulted in a collection of data not amenable to meta-analysis, a common finding illustrating how important it is for case definitions and other criteria to be standardized between studies. Therefore, results are presented according to the SWiM (Synthesis Without Meta-analysis) guidelines. We summarize the main findings reported in the 33 selected articles for the different markers and report their results in form of comparative tables including sample selection criteria, marker values, and diagnostic performances, where available. Finally, we report the study limitations and bias assessment findings.
| S-EPMC9165246 | biostudies-literature
Project description:BackgroundPeople are paying more attention to the healthy and balanced diet with the improvement of their living standards. Milk fatty acids (FAs) have been reported that they were related to some atherosclerosis and coronary heart diseases in human. In our previous genome-wide association study (GWAS) on milk FAs in dairy cattle, 83 genome-wide significant single nucleotide polymorphisms (SNPs) were detected. Among them, two SNPs, ARS-BFGL-NGS-109493 and BTA-56389-no-rs associated with C18index (P = 0.0459), were located in the upstream of 1-acylglycerol-3-phosphate O-acyltransferase 3 (AGPAT3) gene. AGPAT3 is involved in glycerol-lipid, glycerol-phospholipid metabolism and phospholipase D signaling pathways. Hence, it was inferred as a candidate gene for milk FAs. The aim of this study was to further confirm the genetic effects of the AGPAT3 gene on milk FA traits in dairy cattle.ResultsThrough re-sequencing the complete coding region, and 3000 bp of 5' and 3' regulatory regions of the AGPAT3 gene, a total of 17 SNPs were identified, including four in 5' regulatory region, one in 5' untranslated region (UTR), three in introns, one in 3' UTR, and eight in 3' regulatory region. By the linkage disequilibrium (LD) analysis with Haploview4.1 software, two haplotype blocks were observed that were formed by four and 12 identified SNPs, respectively. Using SAS9.2, we performed single locus-based and haplotype-based association analysis on 24 milk FAs in 1065 Chinese Holstein cows, and discovered that all the SNPs and the haplotype blocks were significantly associated with C6:0, C8:0 and C10:0 (P < 0.0001-0.0384). Further, with Genomatix, we predicted that four SNPs in 5' regulatory region (g.146702957G > A, g.146704373A > G, g.146704618A > G and g.146704699G > A) changed the transcription factor binding sites (TFBSs) for transcription factors SMARCA3, REX1, VMYB, BRACH, NKX26, ZBED4, SP1, USF1, ARNT and FOXA1. Out of them, two SNPs were validated to impact transcriptional activity by performing luciferase assay that the alleles A of both SNPs, g.146704373A > G and g.146704618A > G, increased the transcriptional activities of AGPAT3 promoter compared with alleles G (P = 0.0004).ConclusionsIn conclusion, our findings first demonstrated the significant genetic associations of the AGPAT3 gene with milk FAs in dairy cattle, and two potential causal mutations were detected.
| S-EPMC7849138 | biostudies-literature
Project description:BACKGROUND: One major problem in dairy cattle husbandry is the prevalence of udder infections. In today's breeding programmes, top priority is being given to making animal evaluation more cost-effective and reliable and less time-consuming. We proposed tumor necrosis factor ? (TNF-?), lactoferrin (LTF) and macrophage-expressed lysozyme (mLYZ) genes as potential DNA markers in the improvement of immunity to mastitis.This study included 588 Polish Holstein-Friesian cows kept on one farm located in the north-western region of Poland. All clinical cases of mastitis in the herd under study were recorded by a qualified veterinarian employed by the farm. The following indicators were applied to determine udder immunity to mastitis in the cows under study: morbidity rate (MR), duration of mastitis (DM) and extent of mastitis (EM). TNF-?, mLYZ and LTF genotypes were identified by real-time PCR method, using SimpleProbe technology. Due to the very low frequency of mLYZ allele T, the gene was excluded from further analysis.A statistical analysis of associations between TNF-? and LTF genes and immunity to mastitis were performed using three models: 1) a parity-averaged model including only additive effects of the genes; 2) a parity-averaged model including both additive and epistatic effects of the genes; and 3) a parity-specific model including only additive effects of the genes. RESULTS: With the first and second models it was revealed that the genes effects on the applied indicators of immunity to mastitis were non-significant whereas with the third one the effects were found to be statistically significant. Particularly noteworthy was the finding that the effects of TNF-? and LTF varied depending on age (parity). The alleles which were linked to high immunity to mastitis in lower parities appeared to be less favourable in higher parities. CONCLUSIONS: These interactions might be related to inflamm-ageing, that is an increased susceptibility to infection due to immune system deregulation that progresses with age. Such pattern of interactions makes it impossible to use the genes in question in marker-assisted selection aimed at reducing heritable susceptibility to mastitis. This is because the immune mechanisms behind resistance to infections proved to be too complex.
| S-EPMC3682883 | biostudies-other