Dataset Information

Genome-wide association studies for methane emission and ruminal volatile fatty acids using Holstein cattle sequence data.

ABSTRACT:

Background

Methane emission by ruminants has contributed considerably to the global warming and understanding the genomic architecture of methane production may help livestock producers to reduce the methane emission from the livestock production system. The goal of our study was to identify genomic regions affecting the predicted methane emission (PME) from volatile fatty acids (VFAs) indicators and VFA traits using imputed whole-genome sequence data in Iranian Holstein cattle.

Results

Based on the significant-association threshold (p -?8), 33 single nucleotide polymorphisms (SNPs) were detected for PME per kg milk (n?=?2), PME per kg fat (n?=?14), and valeric acid (n?=?17). Besides, 69 genes were identified for valeric acid (n?=?18), PME per kg milk (n?=?4) and PME per kg fat (n?=?47) that were located within 1?Mb of significant SNPs. Based on the gene ontology (GO) term analysis, six promising candidate genes were significantly clustered in organelle organization (GO:0004984, p?=?3.9?×?10^-?2) for valeric acid, and 17 candidate genes significantly clustered in olfactory receptors activity (GO:0004984, p?=?4?×?10^-?10) for PME traits. Annotation results revealed 31 quantitative trait loci (QTLs) for milk yield and its components, body weight, and residual feed intake within 1?Mb of significant SNPs.

Conclusions

Our results identified 33 SNPs associated with PME and valeric acid traits, as well as 17 olfactory receptors activity genes for PME traits related to feed intake and preference. Identified SNPs were close to 31 QTLs for milk yield and its components, body weight, and residual feed intake traits. In addition, these traits had high correlations with PME trait. Overall, our findings suggest that marker-assisted and genomic selection could be used to improve the difficult and expensive-to-measure phenotypes such as PME. Moreover, prediction of methane emission by VFA indicators could be useful for increasing the size of reference population required in genome-wide association studies and genomic selection.

SUBMITTER: Jalil Sarghale A

PROVIDER: S-EPMC7684878 | biostudies-literature | 2020 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Genome-wide association studies for methane emission and ruminal volatile fatty acids using Holstein cattle sequence data.

Jalil Sarghale Ali A Moradi Shahrebabak Mohammad M Moradi Shahrebabak Hossein H Nejati Javaremi Ardeshir A Saatchi Mahdi M Khansefid Majid M Miar Younes Y

BMC genetics 20201123 1

<h4>Background</h4>Methane emission by ruminants has contributed considerably to the global warming and understanding the genomic architecture of methane production may help livestock producers to reduce the methane emission from the livestock production system. The goal of our study was to identify genomic regions affecting the predicted methane emission (PME) from volatile fatty acids (VFAs) indicators and VFA traits using imputed whole-genome sequence data in Iranian Holstein cattle.<h4>Resul ...[more]

PMID: 33228565

Similar Datasets

Project description:Sole ulcers (SUs) and white line disease (WLD) are two common noninfectious claw lesions (NICL) that arise due to a compromised horn production and are frequent causes of lameness in dairy cattle, imposing welfare and profitability concerns. Low to moderate heritability estimates of SU and WLD susceptibility indicate that genetic selection could reduce their prevalence. To identify the susceptibility loci for SU, WLD, SU and/or WLD, and any type of noninfectious claw lesion, genome-wide association studies (GWAS) were performed using generalized linear mixed model (GLMM) regression, chunk-based association testing (CBAT), and a random forest (RF) approach. Cows from five commercial dairies in California were classified as controls having no lameness records and ≥6 years old (n = 102) or cases having SU (n = 152), WLD (n = 117), SU and/or WLD (SU + WLD, n = 198), or any type of noninfectious claw lesion (n = 217). The top single nucleotide polymorphisms (SNPs) were defined as those passing the Bonferroni-corrected suggestive and significance thresholds in the GLMM analysis or those that a validated RF model considered important. Effects of the top SNPs were quantified using Bayesian estimation. Linkage disequilibrium (LD) blocks defined by the top SNPs were explored for candidate genes and previously identified, functionally relevant quantitative trait loci. The GLMM and CBAT approaches revealed the same regions of association on BTA8 for SU and BTA13 common to WLD, SU + WLD, and NICL. These SNPs had effects significantly different from zero, and the LD blocks they defined explained a significant amount of phenotypic variance for each dataset (6.1-8.1%, p < 0.05), indicating the small but notable contribution of these regions to susceptibility. These regions contained candidate genes involved in wound healing, skin lesions, bone growth and mineralization, adipose tissue, and keratinization. The LD block defined by the most significant SNP on BTA8 for SU included a SNP previously associated with SU. The RF models were overfitted, indicating that the SNP effects were very small, thereby preventing meaningful interpretation of SNPs and any downstream analyses. These findings suggested that variants associated with various physiological systems may contribute to susceptibility for NICL, demonstrating the complexity of genetic predisposition.

Project description:Urban sludge produced on wastewater treatment plants are often contaminated by organic pollutants such as polycyclic aromatic hydrocarbons (PAH). Their removal under methanogenic conditions was already reported, but the factors influencing this removal remain unclear. Here, we determined the influence of microbial communities on PAH removal under controlled physico-chemical conditions. Twelve mesophilic anaerobic digesters were inoculated with three microbial communities extracted from ecosystems with contrasting pollution histories: a PAH contaminated soil, a PCB contaminated sediment and a low contaminated anaerobic sludge. These anaerobic digesters were operated during 100 days in continuous mode. A sterilised activated sludge, spiked with 13 PAH at concentrations usually encountered in full-scale wastewater treatment plants, was used as substrate. The dry matter and volatile solid degradation, the biogas production rate and composition, the volatile fatty acids (VFA) production and the PAH removals were monitored. Bacterial and archaeal communities were compared in abundance (qPCR), in community structure (SSCP fingerprinting) and in dominant microbial species (454-pyrosequencing). The bioreactors inoculated with the community extracted from low contaminated anaerobic sludge showed the greater methane production. The PAH removals ranged from 10% to 30%, respectively, for high and low molecular weight PAH, whatever the inoculums tested, and were highly correlated with the dry matter and volatile solid removals. The microbial community structure and diversity differed with the inoculum source; this difference was maintained after the 100 days of digestion. However, the PAH removal was not correlated to these diverse structures and diversities. We hence obtained three functional stable consortia with two contrasted metabolic activities, and three different pictures of microbial diversity, but similar PAH and matter removals. These results confirm that PAH removal depends on the molecule type and on the solid matter removal. But, as PAH elimination is similar whether the solid substrate is degraded into VFA or into methane, it seems that the fermentative communities are responsible for their elimination.

Project description:Mitigation of greenhouse gas emissions is relevant for reducing the environmental impact of ruminant production. In this study, the rumen microbiome from Holstein cows was characterized through a combination of 16S rRNA gene and shotgun metagenomic sequencing. Methane production (CH4 ) and dry matter intake (DMI) were individually measured over 4-6 weeks to calculate the CH4 yield (CH4 y = CH4 /DMI) per cow. We implemented a combination of clustering, multivariate and mixed model analyses to identify a set of operational taxonomic unit (OTU) jointly associated with CH4 y and the structure of ruminal microbial communities. Three ruminotype clusters (R1, R2 and R3) were identified, and R2 was associated with higher CH4 y. The taxonomic composition on R2 had lower abundance of Succinivibrionaceae and Methanosphaera, and higher abundance of Ruminococcaceae, Christensenellaceae and Lachnospiraceae. Metagenomic data confirmed the lower abundance of Succinivibrionaceae and Methanosphaera in R2 and identified genera (Fibrobacter and unclassified Bacteroidales) not highlighted by metataxonomic analysis. In addition, the functional metagenomic analysis revealed that samples classified in cluster R2 were overrepresented by genes coding for KEGG modules associated with methanogenesis, including a significant relative abundance of the methyl-coenzyme M reductase enzyme. Based on the cluster assignment, we applied a sparse partial least-squares discriminant analysis at the taxonomic and functional levels. In addition, we implemented a sPLS regression model using the phenotypic variation of CH4 y. By combining these two approaches, we identified 86 discriminant bacterial OTUs, notably including families linked to CH4 emission such as Succinivibrionaceae, Ruminococcaceae, Christensenellaceae, Lachnospiraceae and Rikenellaceae. These selected OTUs explained 24% of the CH4 y phenotypic variance, whereas the host genome contribution was ~14%. In summary, we identified rumen microbial biomarkers associated with the methane production of dairy cows; these biomarkers could be used for targeted methane-reduction selection programmes in the dairy cattle industry provided they are heritable.

Project description:Sole ulcers (SU) and white line disease (WLD) are noninfectious claw lesions that arise due to compromised horn production and are common causes of lameness in dairy cattle, imposing welfare and profitability concerns. The low to moderate heritability estimates of SU and WLD susceptibility indicate that genetic selection could reduce their prevalence. To identify loci associated with SU and WLD susceptibility, generalized linear mixed model (GLMM) regression and random forest (RF) genome-wide association studies (GWAS) were performed. Cows from five commercial dairies in California were classified as sound controls having no instances of lameness and above six years of age (n = 102) or cases having SU (n = 152), WLD (n = 117), SU and/or WLD (SU+WLD, n = 198), or any type of noninfectious claw lesion (n = 217). Top SNPs were defined as those passing Bonferroni-corrected suggestive and significance thresholds in the GLMM analysis or those that a validated RF model considered important. Effects of top SNPs were quantified using Bayesian estimation. Linkage disequilibrium (LD) blocks defined by top SNPs were explored for candidate genes and previously identified, functionally relevant quantitative trait loci. The GLMM GWAS revealed regions of association on BTA8 for SU and BTA13 common to WLD, SU+WLD, and noninfectious claw lesions. Bayesian estimation of effect sizes indicated that these SNPs had effects significantly different from zero, indicating their small but notable contribution to susceptibility. Promising candidate genes identified in these regions were involved in wound healing, skin lesions, bone growth and mineralization, adipose tissue, and keratinization. The LD block defined by the most significant SNP on BTA8 for SU included a SNP previously reported to be associated with SU. The RF models were overfitted, indicating that SNP effects were very small and thereby preventing meaningful interpretation of SNPs with the highest importance values and downstream analyses. These findings suggested that variants associated with a variety of physiological systems may contribute to susceptibility for noninfectious claw lesions, demonstrating the complexity of genetic predisposition.

Project description:BACKGROUND: Breeding for enhanced immune response (IR) has been suggested as a tool to improve inherent animal health. Dairy cows with superior antibody-mediated (AMIR) and cell-mediated immune responses (CMIR) have been demonstrated to have a lower occurrence of many diseases including mastitis. Adaptive immune response traits are heritable, and it is, therefore, possible to breed for improved IR, decreasing the occurrence of disease. The objective of this study was to perform genome-wide association studies to determine differences in genetic profiles among Holstein cows classified as High or Low for AMIR and CMIR. From a total of 680 cows with immune response phenotypes, 163 cows for AMIR (81 High and 82 Low) and 140 for CMIR (75 High and 65 Low) were selectively genotyped using the Illumina Bovine SNP50 BeadChip. Results were validated using an unrelated population of 164 Holstein bulls IR phenotyped for AMIR and 146 for CMIR. RESULTS: A generalized quasi likelihood score method was used to determine single nucleotide polymorphisms (SNP) and chromosomal regions associated with immune response. After applying a 5% chromosomal false discovery rate, 186 SNPs were significantly associated with AMIR. The majority (93%) of significant markers were on chromosome 23, with a similar peak found in the bull population. For CMIR, 21 SNP markers remained significant. Candidate genes within 250,000 base pairs of significant SNPs were identified to determine biological pathways associated with AMIR and CMIR. Various pathways were identified, including the antigen processing and presentation pathway, important in host defense. Candidate genes included those within the bovine Major Histocompatability Complex such as BoLA-DQ, BoLA-DR and the non-classical BoLA-NC1 for AMIR and BoLA-DQ for CMIR, the complement system including C2 and C4 for AMIR and C1q for CMIR, and cytokines including IL-17A, IL17F for AMIR and IL-17RA for CMIR and tumor necrosis factor for both AMIR and CMIR. Additional genes associated with CMIR included galectins 1, 2 and 3, BCL2 and ?-defensin. CONCLUSIONS: The significant genetic variation associated with AMIR and CMIR in this study may imply feasibility to include immune response in genomic breeding indices as an approach to improve inherent animal health.

Dataset Information

Genome-wide association studies for methane emission and ruminal volatile fatty acids using Holstein cattle sequence data.

Background

Results

Conclusions

Publications

Genome-wide association studies for methane emission and ruminal volatile fatty acids using Holstein cattle sequence data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure