The posttranscriptional gene regulation mediated by microRNA plays an important role in the development and function of male and female reproductive organs and germ cells in mammals, including cattle. In the present study, we identified novel and differentially expressed miRNAs in the testis and ovary in Holstein cattle by combining the Solexa sequencing with bioinformatics. In total 100 and 104 novel pre-miRNAs were identified in testicular and ovarian tissues, encoding 122 and 136 mature miRNAs, respectively. Of these, 6 miRNAs appear to be bovine-specific. A total of 246 known miRNAs were co-expressed in the testicular and ovarian tissues. Of the known miRNAs, twenty-one testis-specific and nine ovary-specific (1-23 reads) were found. Approximately 30.5% of the known bovine miRNAs in this study were found to have >2-fold differential expression within the two respective reproductive organ systems. The putative miRNA target genes of miRNAs were involved in pathways associated with reproductive physiology. Both known and novel tissue-specific miRNAs are expressed by Real-time quantitative PCR analysis in dairy cattle. This study expands the number of miRNAs known to be expressed in cattle. The patterns of miRNAs expression differed significantly between the bovine testicular and ovarian tissues, which provide important information on sex differences in miRNA expression. Diverse miRNAs may play an important regulatory role in the development of the reproductive organs in Holstein cattle.

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.

Proteomics and bioinformatics may help us better understand the biological adaptations occurring during bovine mastitis. This systems approach also could help identify biomarkers for monitoring clinical and subclinical mastitis. The aim of the present study was to use isobaric tags for relative and absolute quantification (iTRAQ) to screen potential proteins associated with mastitis at late infectious stage. Healthy and mastitic cows’ mammary gland tissues were analyzed using iTRAQ combined with two-dimensional liquid chromatography-tandem mass spectrometry (2D-LC-MS/MS). Bioinformatics analyses of differentially expressed proteins were performed by means of Gene Ontology, metabolic pathways, transcriptional regulation networks using Blast2GO software, the Dynamic Impact Approach and Ingenuity Pathway Analysis. At a false discovery rate of 5%, a total of 768 proteins were identified from 6,499 peptides, which were matched with 15,879 spectra. Compared with healthy mammary gland tissue, 36 proteins were significantly up-regulated (>1.5-fold) while 19 were significantly down-regulated (<0.67-fold) in response to mastitis due to natural infections with Staphylococci aureus. Up-regulation of collagen, type I, alpha 1 (COL1A1) and inter-alpha (Globulin) inhibitor H4 (ITIH4) in the mastitis-infected tissue was confirmed by Western blotting and Immunohistochemistry. This paper is the first to show the protein expression in the late response to a mastitic pathogen, thus, revealing mechanisms associated with host tissue damage. The bioinformatics analyses highlighted the effects of mastitis on proteins such as collagen, fibrinogen, fibronectin, casein alpha and heparan sulfate proteoglycan 2. Our findings provide additional clues for further studies of candidate genes for mastitis susceptibility. The up-regulated expression of COL1A1 and ITIH4 in the mastitic mammary gland may be associated with tissue damage and repair during late stages of infection.

The membrane protein CD46, a ubiquitous cell surface pathogen receptor, can bind Streptococcus to trigger cell autophagy, which is a critical step in the control of infection. In this study, we found a new splice variant designated CD46 transcript variant (CD46-TV). The splice variant is characterized by the retention of a 48 bp sequence from intron 8 of the bovine CD46 gene, which encodes a putative protein enlarged by 16 amino acids. CD46-TV mRNA was found to be over expressed in mastitis-infected mammary gland tissues relative to healthy tissues. A single nucleotide polymorphism (c. 1033 + 2184 C > T) in the exonic splicing enhancer (ESE) motif region was shown to result in the CD46-TV aberrant splice variant through constructing alternative alleles using the pSPL3 exon capturing vector and transfecting these into 293 T cells. Allelic frequency in 56,682 individuals belonging to 112 Bos taurus, Bos indicus, Bos javanicus, Bos grunniens and Bos mutus, etc. suggests that the C allele (80.09%) is the ancestral allele. Association analysis found that the mean genomic estimated breeding values (gEBV) for milk somatic cell score and the occurrence of clinical mastitis, as well as the milk somatic cell score of Chinese Holsteins with the CT genotype was lower than those of individuals with either the CC or TT genotypes. The mean gEBV for udder health synthesis for the TT genotype was greater than those for the CC or CT genotypes. Our findings suggest that the CD46 gene likely plays a critical role in the risk of mastitis caused by Streptococcus in dairy cows via an alternative splicing mechanism caused by a functional mutation in intron 8. Our data also underline the importance of variation within ESEs in regulating transcript processing.

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.

Alternative splicing (AS) contributes to the complexity of the mammalian proteome and plays an important role in diseases, including infectious diseases. The differential AS patterns of these transcript sequences between the healthy (HS3A) and mastitic (HS8A) cows naturally infected by Staphylococcus aureus were compared to understand the molecular mechanisms underlying mastitis resistance and susceptibility. In this study, using the Illumina paired-end RNA sequencing method, 1352 differentially expressed genes (DEGs) with higher than twofold changes were found in the HS3A and HS8A mammary gland tissues. Gene ontology and KEGG pathway analyses revealed that the cytokine–cytokine receptor interaction pathway is the most significantly enriched pathway. Approximately 16k annotated unigenes were respectively identified in two libraries, based on the bovine Bos taurus UMD3.1 sequence assembly and search. A total of 52.62% and 51.24% annotated unigenes were alternatively spliced in term of exon skipping, intron retention, alternative 5′ splicing and alternative 3ʹ splicing. Additionally, 1,317 AS unigenes were HS3A-specific, whereas 1,093 AS unigenes were HS8A-specific. Some immune-related genes, such as ITGB6, MYD88, ADA, ACKR1, and TNFRSF1B, and their potential relationships with mastitis were highlighted. From Chromosome 2, 4, 6, 7, 10, 13, 14, 17, and 20, 3.66% (HS3A) and 5.4% (HS8A) novel transcripts, which harbor known quantitative trait locus associated with clinical mastitis, were identified. Many DEGs in the healthy and mastitic mammary glands are involved in immune, defense, and inflammation responses. These DEGs, which exhibit diverse and specific splicing patterns and events, can endow dairy cattle with the potential complex genetic resistance against mastitis.