CFAP58 is involved in the sperm head shaping and flagellogenesis of cattle and mice.

CFAP58 is a testis-enriched gene that plays an important role in the sperm flagellogenesis of humans and mice. However, the effect of CFAP58 on bull semen quality and the underlying molecular mechanisms involved in spermatogenesis remain unknown. Here, we identified two single-nucleotide polymorphisms (rs110610797, A>G and rs133760846, G>T) and one indel (g.-1811_ g.-1810 ins147bp) in the promoter of CFAP58 that were significantly associated with semen quality of bulls, including sperm deformity rate and ejaculate volume. Moreover, by generating gene knockout mice, we found for the first time that the loss of Cfap58 not only causes severe defects in the sperm tail, but also affects the manchette structure, resulting in abnormal sperm head shaping. Cfap58 deficiency causes an increase in spermatozoa apoptosis. Further experiments confirmed that CFAP58 interacts with IFT88 and CCDC42. Moreover, it may be a transported cargo protein that plays a role in stabilizing other cargo proteins, such as CCDC42, in the intra-manchette transport/intra-flagellar transport pathway. Collectively, our findings reveal that CFAP58 is required for spermatogenesis and provide genetic markers for evaluating semen quality in cattle.

Coordinated alternation of DNA methylation and alternative splicing of PBRM1 affect bovine sperm structure and motility.

DNA methylation and gene alternative splicing drive spermatogenesis. In screening DNA methylation markers and transcripts related to sperm motility, semen from three pairs of full-sibling Holstein bulls with high and low motility was subjected to reduced representation bisulphite sequencing. A total of 948 DMRs were found in 874 genes (gDMRs). Approximately 89% of gDMR-related genes harboured alternative splicing events, including SMAD2, KIF17, and PBRM1. One DMR in exon 29 of PBRM1 with the highest 5mC ratio was found, and hypermethylation in this region was related to bull sperm motility. Furthermore, alternative splicing events at exon 29 of PBRM1 were found in bull testis, including PBRM1-complete, PBRM1-SV1 (exon 28 deletion), and PBRM1-SV2 (exons 28-29 deletion). PBRM1-SV2 exhibited significantly higher expression in adult bull testes than in newborn bull testes. In addition, PBRM1 was localized to the redundant nuclear membrane of bull sperm, which might be related to sperm motility caused by sperm tail breakage. Therefore, the hypermethylation of exon 29 may be associated with the production of PBRM1-SV2 in spermatogenesis. These findings indicated that DNA methylation alteration at specific loci could regulate gene splicing and expression and synergistically alter sperm structure and motility.

Signatures of selection in indigenous Chinese cattle genomes reveal adaptive genes and genetic variations to cold climate.

Cold climate shapes the genome of animals and drives them to carry sufficient genetic variations to adapt to changes in temperature. However, limited information is available about the genome-wide pattern of adaptations to cold environments in cattle. In the present study, we used 777K SNP genotyping (15 Chinese cattle breeds, 198 individuals) and whole genome resequencing data (54 cattle breeds of the world, 432 individuals) to disentangle divergent selection signatures, especially between the cold-adapted (annual average temperature of habitat, 6.24 °C to 10.3 °C) and heat-adapted (20.2 °C to 24.73 °C) Chinese native cattle breeds. Genomic analyses revealed a set of candidate genes (e.g., UQCR11, DNAJC18, EGR1, and STING1) were functionally associated with thermogenesis and energy metabolism. We also characterized the adaptive loci of cattle exposed to cold temperatures. Our study finds new candidate genes and provides new insights into adaptations to cold climates in cattle.

Cold climates can affect cattle performance, survival, and health. Local cattle breeds have been adapted to the local environments including extremely cold temperatures after a long period of natural and artificial selection. Selection and local adaptation are shaping populations. However, identifying loci associated with cold adaptation has been a major challenge. We used high-density SNP arrays and resequencing data to comprehensively analyze and compare the genomic selection signatures of Chinese northern and southern cattle, and elucidated several adaptive genes and alleles involved in cold adaptation. The complexity of genetic adaptation mechanism among different low-temperature adapted cattle breeds was also emphasized.

Selection Signature and CRISPR/Cas9-Mediated Gene Knockout Analyses Reveal ZC3H10 Involved in Cold Adaptation in Chinese Native Cattle.

Cold stress is an important factor affecting cattle health, production performance, and reproductive efficiency. Understanding of the potential mechanism underlying genetic adaptation to local environments, particularly extreme cold environment, is limited. Here, by using FLK and hapFLK methods, we found that the Zinc finger CCCH-type containing 10 (ZC3H10) gene underwent positive selection in the Menggu, Fuzhou, Anxi, and Shigatse humped cattle breeds that are distributed in the cold areas of China. Furthermore, ZC3H10 expression significantly increased in bovine fetal fibroblast (BFF) cells at 28 °C for 4 h. ZC3H10 knockout BFFs were generated using CRISPR/Cas9. Wild and ZC3H10-deleted BFFs were treated at two temperatures and were divided into four groups (WT, wild and cultured at 38 °C; KO, ZC3H10-/- and 38 °C; WT_LT, wild, and 28 °C for 4 h; and KO_LT, ZC3H10-/- and 28 °C for 4 h. A total of 466, 598, 519, and 650 differently expressed genes (two-fold or more than two-fold changes) were identified by determining transcriptomic difference (KO_LT vs. KO, WT_LT vs. WT, KO vs. WT, and KO_LT vs. WT_LT, respectively). Loss of ZC3H10 dysregulated pathways involved in thermogenesis and immunity, and ZC3H10 participated in immunity-related pathways induced by cold stress and regulated genes involved in glucose and lipid metabolism and lipid transport (PLTP and APOA1), thereby facilitating adaptability to cold stress. Our findings provide a foundation for further studies on the function of ZC3H10 in cold stress and development of bovine breeding strategies for combatting the influences of cold climate.

Complement component 3 haplotypes influence serum complement activity and milk production traits in Chinese Holstein cattle.

Complement component 3 (C3) is the key molecule of the three pathways of complement activation (alternative, classical, and lectin pathways), which are involved in phagocytosis, inflammation, and immunoregulation processes to destroy infectious microorganisms. In this study, three novel single-nucleotide polymorphisms (SNPs) (g.-1293C>G located in the 5'-flanking region, g.56T>C in exon I, and g.7017C>T in exon XII) of the C3 gene were detected using created restriction site polymerase chain reaction, restriction fragment length polymorphism, and DNA sequencing in 952 cattle from three Chinese breeds. The genotypes and haplotypes were analyzed to investigate the polymorphisms and their possible implications, with particular investigative focus on their associations with serum C3 level, complement hemolytic activity (CH50 and ACH50), and milk production traits. The g.56T>C SNP in exon I affected the serum ACH50 (P<0.01) and the milk somatic cell score (SCS) (P<0.05), and the g.7017C>T SNP in exon XII significantly affected the serum ACH50 values (P<0.01). Moreover, statistical analyses revealed that individuals with genotypic combination CCC/GCC showed significantly lower SCS and the lowest C3 concentration in serum compared with cows with CCC/GTT (P = 0.0007) and CTT/CTT (P = 0.0021); the individuals with CCC/CCT had significantly higher ACH50 values than those with CCC/CTC (P = 0.0008) and CTC/GTC (P = 0.001); cows with CCT/CTT had higher values of CH50 and 305-day milk yield (P>0.05). The C3 expression levels were significantly increased in lung and mammary tissues (P<0.05), while significantly decreased in heart, spleen, liver, and kidney tissues in mastitis cows compared with those in healthy animals (P<0.01), respectively. Bacterial counts of serum antibacterial activities were also completed to verify the effect of SNPs on resistance to mastitis pathogens. Genetically resistant cows (CCC/GCC) had serum with noticeably higher antibacterial activity against S. aureus and E. coli in vitro than the genetically susceptible CCC/GTT cows (P<0.05). Results from this study imply that the C3 gene plays a role in resistance to bacterial infection and that it can be used as a molecular marker for complement activity and traits related to milk production.

Introgression, admixture, and selection facilitate genetic adaptation to high-altitude environments in cattle.

Domestication and subsequent selection of cattle to form breeds and biological types that can adapt to different environments partitioned ancestral genetic diversity into distinct modern lineages. Genome-wide selection particularly for adaptation to extreme environments left detectable signatures genome-wide. We used high-density genotype data for 42 cattle breeds and identified the influence of Bos grunniens and Bos javanicus on the formation of Chinese indicine breeds that led to their divergence from India-origin zebu. We also found evidence for introgression, admixture, and migration in most of the Chinese breeds. Selection signature analyses between high-altitude (≥1800 m) and low-altitude adapted breeds (<1500 m) revealed candidate genes (ACSS2, ALDOC, EPAS1, EGLN1, NUCB2) and pathways that are putatively involved in hypoxia adaptation. Immunohistochemical, real-time PCR and CRISPR/cas9 ACSS2-knockout analyses suggest that the up-regulation of ACSS2 expression in the liver promotes the metabolic adaptation of cells to hypoxia via the hypoxia-inducible factor pathway. High altitude adaptation involved the introgression of alleles from high-altitude adapted yaks into Chinese Bos taurus taurus prior to their formation into recognized breeds and followed by selection. In addition to selection, adaptation to high altitude environments has been facilitated by admixture and introgression with locally adapted cattle populations.

Effect of functional single nucleotide polymorphism g.-572 A > G of apolipoprotein A1 gene on resistance to ketosis in Chinese Holstein cows.

The ketosis has negative effects on the high-yielding dairy cows during early lactation. Apolipoprotein A1 (APOA1) is a component of high-density lipoprotein. However, the association of APOA1 gene with ketosis, and the molecular mechanisms of expression of APOA1 gene are not fully understood in dairy cows. In this study, expression of APOA1 in the liver and blood was investigated by RT-qPCR and immunohistochemistry, and genetic variation in the 5'-flanking region of the AOPA1 gene was also screened and identified. In addition, correlation of the single nucleotide polymorphisms (SNPs) of APOA1 gene with blood ketone characters, and activity of APOA1 promoter were analyzed in dairy cows. The results showed that ApoA1 protein was expressed in the liver, and the mRNA level of APOA1 was significantly higher in the cows with ketosis comparing to the healthy cows. In addition, a novel SNP (g.-572 A > G) in the core promoter of the APOA1 gene was identified between base g.-714 and g.-68 through transient transfection in both HepG2 cell and FFb cell, and luciferase report assay. Moreover, there was lower concentration of blood ß-hydroxybutyrate in cows with genotype GG comparing to the cows with genotypes AA and AG. This study reported for the first time that the genetic variant g.-572 A > G in the core promoter region of APOA1 gene was associated with the ketosis in Chinese Holstein cows, and g.-572 A > G may be used as a genetic marker for ketosis prevention.

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.

Genome-wide methylation and transcriptome of blood neutrophils reveal the roles of DNA methylation in affecting transcription of protein-coding genes and miRNAs in E. coli-infected mastitis cows.

BACKGROUND: Neutrophils are the first effectors of inflammatory response triggered by mastitis infection, and are important defense cells against pathogenic Escherichia coli (E. coli). DNA methylation, as a critical epigenetic mechanism for regulating gene function, is involved in bovine mastitis. RESULTS: In this study, we sequenced the blood neutrophils of healthy and E. coli-infected mastitic half-sib cows for the overall DNA methylation levels using transcriptome sequencing and reduced representation bisulfite sequencing. The methylation levels in the mastitis cows (MCs) were decreased compared with healthy cows (HCs). A total of 494 differentially methylated regions were identified, among which 61 were up-methylated and 433 were down-methylated (MCs vs. HCs). The expression levels of 1094 differentially expressed genes were up-regulated, and 245 genes were down-regulated. Twenty-nine genes were found in methylation and transcription data, among which seven genes' promoter methylation levels were negatively correlated with expression levels, and 11 genes were differentially methylated in the exon regions. The bisulfite sequencing PCR and quantitative real-time PCR validation results demonstrated that the promoter methylation of CITED2 and SLC40A1 genes affected differential expression. The methylation of LGR4 exon 5 regulated its own alternative splicing. The promoter methylation of bta-miR-15a has an indirect effect on the expression of its target gene CD163. The CITED2, SLC40A1, and LGR4 genes can be used as candidates for E. coli-induced mastitis resistance. CONCLUSIONS: This study explored the roles of DNA methylation in affecting transcription of protein-coding genes and miRNAs in E. coli-induced mastitis, thereby helping explain the function of DNA methylation in the pathogenesis of mastitis and provided new target genes and epigenetic markers for mastitis resistance breeding in dairy cattle.

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.

Production of an aberrant splice variant of CCL5 is not caused by genetic mutation in the mammary glands of mastitis­infected Holstein cows.

Genetic mutations, including single nucleotide polymorphisms (SNPs), result in aberrant alternatively splicing of gene and involves in susceptibility of inflammatory diseases, including bovine mastitis. C­C motif chemokine ligand 5 (CCL5) is an immune­associated gene, but its alternative splicing (AS) mechanism of gene expression has not yet been understood. The present study identified the splice variant of CCL5 and the compared differential expression of various transcripts between healthy and mastitic mammary gland tissue from cows. A novel transcript lacking exon 2 with a deletion of 112 bp (referred to as CCL5­AS) was identified in the mammary gland. The expression of CCL5­AS was lower compared with CCL5­reference in the healthy and mastitic mammary tissues. A total of two novel SNPs (g.1647 C>T and g.1804 G>A) were identified in exon 2 of CCL5. Using the splicing mini­gene reporter assay in bovine mammary epithelial cells (MAC­T) and 293T cells, it was confirmed that the production of CCL5­AS was not caused by the two SNPs. The present findings suggested that alternative splicing is one of the mechanisms of CCL5 expression regulation and is involved in mastitis infection, but that genetic mutation was not responsible for the generation of the abnormal transcript of CCL5 in cows.

Integrated analysis of mRNAs and long noncoding RNAs in the semen from Holstein bulls with high and low sperm motility.

Sperm motility is the main index used to assess the quality of bull semen. It may also be used to evaluate the fertility potential of bulls. Protein-coding mRNA and long noncoding RNA (lncRNA) participate in the regulation of spermatogenesis. Here, we employed strand-specific RNA sequencing to profile the semen transcriptome (mRNA and lncRNA) of six paired full-sibling Holstein bulls with divergent sperm motility and to determine the functions of mRNA and lncRNA in sperm motility. Among 20,875 protein-encoding genes detected in semen, 19 were differentially expressed between the high sperm motility group (H: H1, H2, and H3) and low sperm motility group (L: L1, L2, and L3). Of the 11,561 lncRNAs identified in sperm, 2,517 were differentially expressed between the H and L groups. We found that TCONS_00041733 lncRNA targets the node gene EFNA1 (ephrin A1), involved in male reproductive physiology. Our study provides a global mRNA and lncRNA transcriptome of bull semen, as well as novel insights into the regulation of neighboring protein coding by lncRNAs and the influence of mRNAs on sperm motility.

Population Structure, and Selection Signatures Underlying High-Altitude Adaptation Inferred From Genome-Wide Copy Number Variations in Chinese Indigenous Cattle.

Copy number variations (CNVs) have been demonstrated as crucial substrates for evolution, adaptation and breed formation. Chinese indigenous cattle breeds exhibit a broad geographical distribution and diverse environmental adaptability. Here, we analyzed the population structure and adaptation to high altitude of Chinese indigenous cattle based on genome-wide CNVs derived from the high-density BovineHD SNP array. We successfully detected the genome-wide CNVs of 318 individuals from 24 Chinese indigenous cattle breeds and 37 yaks as outgroups. A total of 5,818 autosomal CNV regions (683 bp-4,477,860 bp in size), covering ~14.34% of the bovine genome (UMD3.1), were identified, showing abundant CNV resources. Neighbor-joining clustering, principal component analysis (PCA), and population admixture analysis based on these CNVs support that most Chinese cattle breeds are hybrids of Bos taurus taurus (hereinafter to be referred as Bos taurus) and Bos taurus indicus (Bos indicus). The distribution patterns of the CNVs could to some extent be related to the geographical backgrounds of the habitat of the breeds, and admixture among cattle breeds from different districts. We analyzed the selective signatures of CNVs positively involved in high-altitude adaptation using pairwise Fst analysis within breeds with a strong Bos taurus background (taurine-type breeds) and within Bos taurus×Bos indicus hybrids, respectively. CNV-overlapping genes with strong selection signatures (at top 0.5% of Fst value), including LETM1 (Fst = 0.490), TXNRD2 (Fst = 0.440), and STUB1 (Fst = 0.420) within taurine-type breeds, and NOXA1 (Fst = 0.233), RUVBL1 (Fst = 0.222), and SLC4A3 (Fst=0.154) within hybrids, were potentially involved in the adaptation to hypoxia. Thus, we provide a new profile of population structure from the CNV aspects of Chinese indigenous cattle and new insights into high-altitude adaptation in cattle.

Functional haplotypes of ARID4A affect promoter activity and semen quality of bulls.

The AT-rich interaction domain 4 A (ARID4A) has an important role in regulating Sertoli cell function and male fertility. Its molecular mechanisms, however, remain largely unknown. In this study, two single nucleotide polymorphisms (SNPs) (g.53 G > T, ss 1966531596, and g.826 G > A, rs 210809648) were identified in the promoter region of ARID4A in 215 Chinese Holstein bulls using polymerase chain reaction (PCR)-restriction fragment length polymorphism and created restriction site-PCR. Results revealed that bulls with g.53 G > T-GG and g.826 G > A-G G genotype exhibited higher sperm deformity rate than those with g.53 G > T-TT and g.826 G > A-AA genotype (P < 0.01). Furthermore, three haplotypes (H1 (GG), H3 (TG), H4 (TA)) and six haplotype combinations (H1H1, H1H3, H1H4, H3H3, H3H4, H4H4) were obtained. The bulls with H4H4 exhibited lower sperm deformity rate than those with H1H1 and H1H3 (P < 0.05). In addition, results of bioinformatics analysis revealed that ARID4A has two promoters and that two SNPs of ARID4A are located in transcription factor binding sites. Compared with g.53 G > T-G and g.826 G > A-G allele, there was a greater fluorescence intensity in g.53 G > T-T and g.826 G > A-A allele by transient transfection in MLTC-1 cells and the luciferase report assay. qRT-PCR indicated the ARID4A expression was greater in bull spermatozoa with H4H4 haplotype combination than those with H1H1 haplotype combination (P < 0.05). Results of the present study indicate that g.53 G > T and g.826 G > A are functional mutations that are involved in regulation of ARID4A gene expression by affecting promoter activity and thus semen quality of Chinese Holstein bulls.

Identification of bta-miR-15a∼16a cluster expression, localization and regulated target in Holsteins.

Bovine mastitis is an inflammation response of the mammary gland tissues caused mainly by pathogenic bacteria in cows. Previous studies showed that bta-miR-15a and bta-miR-16a modulate immunity and inflammation responses. In this study, we investigated the expression pattern and tissue localization of bta-miR-15a and bta-miR-16a. The expression levels of bta-miR-15a and bta-miR-16a were significantly upregulated in mammary tissues and blood neutrophils of mastitis-infected cows, compared with those of healthy cows (P < 0.05). Through in situ hybridization, we examined the tissue localization of bta-miR-15a and bta-miR-16a and found that they were expressed in the ductal and acinar cells of mammary gland tissues, where they had a stronger expression signal in the mammary tissues of cows with mastitis than that in healthy cows' tissues. Moreover, we identified CD163 as the target gene of bta-miR-15a and bta-miR-16a. Luciferase assay indicated that bta-miR-15a, bta-miR-16a, and bta-miR-15a∼16a cluster led to the significant reduction in the luciferase activity of CD163 3'UTR vector (P < 0.05). Meanwhile, the luciferase activity had a significantly low value compared with that of single bta-miR-15a or bta-miR-16a plasmid (P < 0.05) in the presence of bta-miR-15a∼16a cluster. The bta-miR-15a∼16a cluster may function more effectively in inhibiting luciferase reporter gene activity of target gene CD163 than single miRNA. Our study provides an insight into the relationship between bovine mastitis and gene expression of bta-miR-15a/16a, which suggested that bta-miR-15a∼16a cluster may play a role against mastitis by binding to target CD163 gene in Holstein dairy cattle.

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.

Solexa sequencing and custom microRNA chip reveal repertoire of microRNAs in mammary gland of bovine suffering from natural infectious mastitis.

Identification of microRNAs (miRNAs), target genes and regulatory networks associated with innate immune and inflammatory responses and tissue damage is essential to elucidate the molecular and genetic mechanisms for resistance to mastitis. In this study, a combination of Solexa sequencing and custom miRNA chip approaches was used to profile the expression of miRNAs in bovine mammary gland at the late stage of natural infection with Staphylococcus aureus, a widespread mastitis pathogen. We found 383 loci corresponding to 277 known and 49 putative novel miRNAs, two potential mitrons and 266 differentially expressed miRNAs in the healthy and mastitic cows' mammary glands. Several interaction networks and regulators involved in mastitis susceptibility, such as ALCAM, COL1A1, APOP4, ITIH4, CRP and fibrinogen alpha (FGA), were highlighted. Significant down-regulation and location of bta-miR-26a, which targets FGA in the mastitic mammary glands, were validated using quantitative real-time PCR, in situ hybridization and dual-luciferase reporter assays. We propose that the observed miRNA variations in mammary glands of mastitic cows are related to the maintenance of immune and defense responses, cell proliferation and apoptosis, and tissue injury and healing during the late stage of infection. Furthermore, the effect of bta-miR-26a in mastitis, mediated at least in part by enhancing FGA expression, involves host defense, inflammation and tissue damage.

DNA methylation rather than single nucleotide polymorphisms regulates the production of an aberrant splice variant of IL6R in mastitic cows.

Interleukin-6 receptor-alpha (IL6R) interacts with IL6 and forms a ligand-receptor complex, which can stimulate various cellular responses, such as cell proliferation, cell differentiation, and activation of inflammatory processes. Both genetic mutation and epigenetic modification regulate gene transcription. We identified a novel splice variant of bovine IL6R, designated as IL6R-TV, which is characterized by the skipping of exon 2 of the NCBI-referenced IL6R gene (IL6R-reference). The expression levels of IL6R-TV and IL6R-reference transcripts were lower in normal mammary gland tissues. These transcripts play a potential role during inflammatory infection. We also detected two putative functional SNPs (g.19711 T > C and g.19731 G > C) located within the upstream 100 bp of exon 2. These SNPs formed two haplotypes (T-G and C-C). Two mutant pSPL3 exon-trapping plasmids (pSPL3-T-G and pSPL3-C-C) were transferred into the bovine mammary epithelial cells (MAC-T) and human embryonic kidney 293 T cells (HEK293T) to investigate the relationship between the two SNPs and the aberrant splicing of IL6R. DNA methylation levels of the alternatively spliced exon in normal and mastitis-infected mammary gland tissues were quantified through nested bisulfate sequencing PCR (BSP) and cloning sequencing. We found that DNA methylation regulated IL6R transcription. The DNA methylation level was high in mastitis-infected mammary gland tissues and stimulated IL6R expression, thereby promoting the inclusion of the alternatively spliced exon. The upregulated expression of the two transcripts was due to DNA methylation modification rather than genetic mutations.

Splicing-related single nucleotide polymorphism of RAB, member of RAS oncogene family like 2B (RABL2B) jeopardises semen quality in Chinese Holstein bulls.

RAB, member of RAS oncogene family like 2B (RABL2B) is a member of a poorly characterised clade of the RAS GTPase superfamily, which plays an essential role in male fertility, sperm intraflagellar transport and tail assembly. In the present study, we identified a novel RABL2B splice variant in bovine testis and spermatozoa. This splice variant, designated RABL2B-TV, is characterised by exon 2 skipping. Moreover, a single nucleotide polymorphism (SNP), namely c.125G>A, was found within the exonic splicing enhancer (ESE) motif, indicating that the SNP caused the production of the RABL2B-TV aberrant splice variant. This was demonstrated by constructing a pSPL3 exon capturing vector with different genotypes and transfecting these vectors into murine Leydig tumour cell line (MLTC-1) cells. Expression of the RABL2B-TV transcript was lower in semen from high- versus low-performance bulls. Association analysis showed that sperm deformity rate was significantly lower in Chinese Holstein bulls with the GG or GA genotype than in bulls with the AA genotype (P<0.05). In addition, initial sperm motility was significantly higher in individuals with the GG or GA genotype than in individuals with the AA genotype (P<0.05). The findings of the present study suggest that the difference in semen quality in bulls with different RABL2B genotypes is generated via an alternative splicing mechanism caused by a functional SNP within the ESE motif.