Effects of eight InDel variants in FHIT on milk traits in Xinjiang brown cattle

Identification and Validation of Genus/Species-Specific Short InDels in Dairy Ruminants

BACKGROUND

Over the past thirty years, the identification of species-specific molecular markers has significantly advanced our understanding of genetic diversity in both plants and animals. Among these, short InDels have emerged as vital genomic features, contributing more to sequence divergence than single nucleotide polymorphisms do in closely related species. This study aimed to identify specific InDels for Bos taurus, Bubalus bubalis, Capra hircus, and Ovis aries via an in silico approach and validated them in 400 individuals (100 for each species).

RESULTS

We identified and characterized short, specific InDels in the sequences of the CSN1S1, CSN1S2, MSTN, and PRLR genes, which can be used for species identification of Capra hircus, Ovis aries, Bos taurus, and Bubalus bubalis, respectively. We developed a Tetraplex Specific PCR assay to enable efficient discrimination among these species.

CONCLUSIONS

This study highlights the utility of InDels as biallelic, codominant markers that are cost-effective and easy to analyse, providing valuable tools for genetic diversity analysis and species identification.

Genome-Wide Association Study of Insertions and Deletions Identified Novel Loci Associated with Milk Production Traits in Dairy Cattle

Genome-wide association study (GWAS) have identified a large number of SNPs associated with milk production traits in dairy cattle. Behind SNPs, INDELs are the second most abundant genetic polymorphisms in the genome, which may exhibit an independent association with complex traits in humans and other species. However, there are no reports on GWASs of INDELs for milk production traits in dairy cattle. In this study, using imputed sequence data, we performed INDEL-based and SNP-based GWASs for milk production traits in a Holstein cattle population. We identified 58 unique significant INDELs for one or multiple traits. The majority of these INDELs are in considerable LD with nearby significant SNPs. However, through conditional association analysis, we identified nine INDELs which showed independent associations. Genomic annotations of these INDELs indicated some novel associated genes, i.e., TRNAG-CCC, EPPK1, PPM1K, PTDSS1, and mir-10163, which were not reported in previous SNP-based GWASs. Our findings suggest that INDEL-based GWASs could be valuable complement to SNP-based GWASs for milk production traits.

DNA methylation patterns in the peripheral blood of Xinjiang brown cattle with variable somatic cell counts

The use of wide-ranging dairy herd improvement (DHI) measurements has resulted in the investigation of somatic cell count (SCC) and the identification of many genes associated with mastitis resistance. In this study, blood samples of Xinjiang brown cattle with different SCCs were collected, and genome-wide DNA methylation was analyzed by MeDIP-seq. The results showed that peaks were mostly in intergenic regions, followed by introns, exons, and promoters. A total of 1,934 differentially expressed genes (DEGs) associated with mastitis resistance in Xinjiang brown cattle were identified. The enrichment of differentially methylated CpG islands of the TRAPPC9 and CD4 genes was analyzed by bisulfate genome sequencing. The methylation rate of differentially methylated CpGs was higher in the TRAPPC9 gene of cattle with clinical mastitis (mastitis group) compared with healthy cattle (control group), while methylation of differentially methylated CpGs was significantly lower in CD4 of the mastitis group compared with the control group. RT-qRCR analysis showed that the mastitis group had significantly reduced expression of CD4 and TRAPPC9 genes compared to the control group (p < 0.05). Furthermore, Mac-T cells treated with lipopolysaccharide and lipoteichoic acid showed significant downregulation of the TRAPPC9 gene in the mastitis group compared with the control group. The identified epigenetic biomarkers provide theoretical reference for treating cow mastitis, breeding management, and the genetic improvement of mastitis resistance in Xinjiang brown cattle.

Transcriptomic Changes and Regulatory Networks Associated with Resistance to Mastitis in Xinjiang Brown Cattle

Xinjiang brown cattle are highly resistant to disease and tolerant of roughage feeding. The identification of genes regulating mastitis resistance in Xinjiang brown cattle is a novel means of genetic improvement. In this study, the blood levels of IL-1β, IL-6, IL-10, TNF-α, and TGF-β in Xinjiang brown cattle with high and low somatic cell counts (SCCs) were investigated, showing that cytokine levels were higher in cattle with high SCCs. The peripheral blood transcriptomic profiles of healthy and mastitis-affected cattle were constructed by RNA-seq. Differential expression analysis identified 1632 differentially expressed mRNAs (DE-mRNAs), 1757 differentially expressed lncRNAs (DE-lncRNAs), and 23 differentially expressed circRNAs (DE-circRNAs), which were found to be enriched in key pathways such as PI3K/Akt, focal adhesion, and ECM-receptor interactions. Finally, ceRNA interaction networks were constructed using the differentially expressed genes and ceRNAs. It was found that keynote genes or mRNAs were also enriched in pathways such as PI3K-Akt, cholinergic synapses, cell adhesion molecules, ion binding, cytokine receptor activity, and peptide receptor activity, suggesting that the key genes and ncRNAs in the network may play an important role in the regulation of bovine mastitis.

Identification of polymorphic loci in OSMR and GHR genes and analysis of their association with growth traits in sheep

The aim of this study was to analyze the effect of OSMR and GHR genes polymorphisms on growth traits in sheep. The single nucleotide polymorphisms of sheep OSMR and GHR genes were identified by DNA sequencing technology. A total of two intronic mutations g.2443 T > C and g.170196 A > G were identified in OSMR and GHR, respectively. Correlation analysis was carried out between the obtained genotypes and the growth traits of sheep. The results showed that at the OSMR g.2443 T > C locus, the body weight, chest circumference and cannon circumference of the TT genotype sheep were significantly higher than those of the CC genotype sheep (p < .05). At the GHR g.170196 A > G locus, the body weight, body length, chest circumference and cannon circumference of the AA genotype sheep were significantly higher than those of the AG genotype and GG genotype sheep (p < .05). Moreover, the body weight of sheep of combination TTOSMR/AAGHR genotype was significantly higher than that of other combination genotypes (p < .05). Therefore, we believe that the polymorphic sites identified in the OSMR and GHR genes can be used as candidate molecular markers for breeding good traits in sheep.

Methylation Levels in the Promoter Region of FHIT and PIAS1 Genes Associated with Mastitis Resistance in Xinjiang Brown Cattle

Mastitis causes serious economic losses in the dairy industry, but there are no effective treatments or preventive measures. In this study, the ZRANB3, PIAS1, ACTR3, LPCAT2, MGAT5, and SLC37A2 genes in Xinjiang brown cattle, which are associated with mastitis resistance, were identified using a GWAS. Pyrosequencing analysis showed that the promoter methylation levels of the FHIT and PIAS1 genes in the mastitis group were higher and lower, respectively, than those in the healthy group (65.97 ± 19.82% and 58.00 ± 23.52%). However, the methylation level of the PIAS1 gene promoter region in the mastitis group was lower than that in the healthy group (11.48 ± 4.12% and 12.17 ± 4.25%). Meanwhile, the methylation levels of CpG3, CpG5, CpG8, and CpG15 in the promoter region of the FHIT and PIAS1 genes in the mastitis group were significantly higher than those in the healthy group (p < 0.01), respectively. RT-qPCR showed that the expression levels of the FHIT and PIAS1 genes were significantly higher in the healthy group than those in the mastitis group (p < 0.01). Correlation analysis showed that the promoter methylation level of the FHIT gene was negatively correlated with its expression. Hence, increased methylation in the promoter of the FHIT gene reduces the mastitis resistance in Xinjiang brown cattle. Finally, this study provides a reference for the molecular-marker-assisted selection of mastitis resistance in dairy cattle.

A 67-bp variable duplication in the promoter region of the ADIPOQ is associated with milk traits in Xinjiang brown cattle

Adiponectin, also known as ADIPOQ, is a hormone protein secreted by adipocytes. The ADIPOQ gene is expressed primarily in adipose tissue, and the encoded protein circulates in the bloodstream and has the potential to regulate both animal fat metabolism and hormone production. Our previous work uncovered a 67-bp variable duplication in the promoter region of ADIPOQ, which reduced the basal transcriptional activity of ADIPOQ in the 3T3_L1 cell and also inhibits the ADIPOQ mRNA expression in adipose tissue. Accordingly, the present study aimed to identify the relationship between the 67-bp structural variations in ADIPOQ promoter region and the milk traits of Xinjiang brown cattle (XJBC). The results revealed two genotypes, DD and ID, in the XJBC, and minor allelic frequency (MAF) for the 'I' allele was more than 1%. Moreover, the association analysis revealed that the 67-bp duplication in the promoter region of the ADIPOQ gene was significantly correlated with the 305 days of milk production volume, fat yield, and milk fat percentage in the XJBC (p < 0.05). These results obtained in this study suggested that the identified variable duplication could be considered as the potential genetic marker for improving milk traits of XJBC.

Two Different Copy Number Variations of the SOX5 and SOX8 Genes in Yak and Their Association with Growth Traits

Copy number variation (CNV) is a structural variant with significant impact on genetic diversity. CNV has been widely used in breeding for growth traits, meat production or quality, and coat color. SRY-like box genes (SOXs) are a class of transcription factors that play a regulatory role in cell fate specification and differentiation. SOX5 and SOX8 belong to subgroups D and E of the SOXs, respectively. Previous studies have shown that SOX5 and SOX8 are essential in the development of bones. In this study, we explored the association between the growth traits and CNVs of SOX5 and SOX8 in 326 Ashidan yaks and detected mRNA expression levels in different tissues. Our results illustrated that CNVs of SOX5 and SOX8 were significantly associated with withers height at 18 months of age and chest girth at 30 months of age (p < 0.05). The CNV combination of SOX5 and SOX8 was significantly associated with withers height at 18 months of age (p < 0.01). SOX5 expression in the lung was significantly higher than in the heart, spleen, kidney, and muscle (p < 0.05). SOX8 expression in the lung was significantly higher than in the liver and muscle (p < 0.05). Our results provide evidence that the CNVs of SOX5 and SOX8 genes could be used as new markers for the selection of yak growth traits.

Combined transcriptome and proteome analyses reveal differences in the longissimus dorsi muscle between Kazakh cattle and Xinjiang brown cattle

Objective

With the rapid development of proteomics sequencing and RNA sequencing technology, multi-omics analysis has become a current research hotspot. Our previous study indicated that Xinjiang brown cattle have better meat quality than Kazakh cattle. In this study, Xinjiang brown cattle and Kazakh cattle were used as the research objects.

Methods

Proteome sequencing and RNA sequencing technology were used to analyze the proteome and transcriptome of the longissimus dorsi muscle of the two breeds of adult steers (n = 3).

Results

In this project, 22,677 transcripts and 1,874 proteins were identified through quantitative analysis of the transcriptome and proteome. By comparing the identified transcriptome and proteome, we found that 1,737 genes were identified at both the transcriptome and proteome levels. The results of the study revealed 12 differentially expressed genes and proteins: troponin I1, crystallin alpha B, cysteine, and glycine rich protein 3, phosphotriesterase-related, myosin-binding protein H, glutathione s-transferase mu 3, myosin light chain 3, nidogen 2, dihydropyrimidinase like 2, glutamate-oxaloacetic transaminase 1, receptor accessory protein 5, and aspartoacylase. We performed functional enrichment of these differentially expressed genes and proteins. The Kyoto encyclopedia of genes and genomes results showed that these differentially expressed genes and proteins are enriched in the fatty acid degradation and histidine metabolism signaling pathways. We performed parallel reaction monitoring (PRM) verification of the differentially expressed proteins, and the PRM results were consistent with the sequencing results.

Conclusion

Our study provided and identified the differentially expressed genes and proteins. In addition, identifying functional genes and proteins with important breeding value will provide genetic resources and technical support for the breeding and industrialization of new genetically modified beef cattle breeds.

Determination of genetic variation within the DYRK2 gene and its associations with milk traits in cattle

To speed up the progress of marker-assisted selection (MAS) in cattle breeding, the dual-specificity tyrosine phosphorylation-regulated kinase 2 (DYRK2), cadherin 2 (CDH2), and kinesin family member 1A (KIF1A) genes were chosen based on our pervious genome-wide association study (GWAS) analysis results. DYRK2 is a kinase that may participate in cell growth and/or development; it shows phosphorylation activity toward serine, threonine, and tyrosine fragments of proteins, and it is different from other protein kinases. The CDH2 gene encodes a classic cadherin, which is a member of the cadherin superfamily. The protein encoded by KIF1A is a member of the kinesin family and plays a role in the transportation of membrane organelles along axon microtubules. We detected insertion/deletion (InDel) variation in these three candidate genes in 438 individual cattle (Xinjiang Brown cattle and Wagyu × Luxi crossbreed cattle). Only DYRK2-P3-11 bp was polymorphic and genotyped. The polymorphism information content of DYRK2-P3-11 bp was 0.336. Correlation analyses showed that InDel polymorphism was significantly associated with six different milk traits. These findings may aid future analyses of InDel genotypes in cattle breeds, and speed up the progress of MAS in cattle breeding.