Identification and expression pattern of two novel alternative splicing variants of EEF1D gene of dairy cattle

The Experimental and In Silico-Based Evaluation of NRF2 Modulators, Sulforaphane and Brusatol, on the Transcriptome of Immortalized Bovine Mammary Alveolar Cells

Changes during the production cycle of dairy cattle can leave these animals susceptible to oxidative stress and reduced antioxidant health. In particular, the periparturient period, when dairy cows must rapidly adapt to the sudden metabolic demands of lactation, is a period when the production of damaging free radicals can overwhelm the natural antioxidant systems, potentially leading to tissue damage and reduced milk production. Central to the protection against free radical damage and antioxidant defense is the transcription factor NRF2, which activates an array of genes associated with antioxidant functions and cell survival. The objective of this study was to evaluate the effect that two natural NRF2 modulators, the NRF2 agonist sulforaphane (SFN) and the antagonist brusatol (BRU), have on the transcriptome of immortalized bovine mammary alveolar cells (MACT) using both the RT-qPCR of putative NRF2 target genes, as well as RNA sequencing approaches. The treatment of cells with SFN resulted in the activation of many putative NRF2 target genes and the upregulation of genes associated with pathways involved in cell survival, metabolism, and antioxidant function while suppressing the expression of genes related to cellular senescence and DNA repair. In contrast, the treatment of cells with BRU resulted in the upregulation of genes associated with inflammation, cellular stress, and apoptosis while suppressing the transcription of genes involved in various metabolic processes. The analysis also revealed several novel putative NRF2 target genes in bovine. In conclusion, these data indicate that the treatment of cells with SFN and BRU may be effective at modulating the NRF2 transcriptional network, but additional effects associated with cellular stress and metabolism may complicate the effectiveness of these compounds to improve antioxidant health in dairy cattle via nutrigenomic approaches.

Detection of candidate polymorphisms around the QTL for fat area ratio to rib eye area on BTA7 using whole-genome resequencing in Japanese Black cattle

In our previous study, we performed genome-wide association study (GWAS) to identify the genomic region associated with Fat area ratio to rib eye area (FAR) and detected a candidate in BTA7 at 10-30 Mbp. The present study aims to comprehensively detect all polymorphisms in the candidate region using whole-genome resequencing data. Based on whole-genome resequencing of eight animals, we detected 127,090 polymorphisms within the region. Of these, 31,945 were located within the genes. We further narrowed the polymorphisms to 6,044 with more than five allele differences between the high and low FAR groups that were located within 179 genes. We subsequently investigated the functions of these genes and selected 170 polymorphisms in eight genes as possible candidate polymorphisms. We focused on SLC27A6 K81M as a putative candidate polymorphism. We genotyped the SNP in a Japanese Black population (n = 904) to investigate the effect on FAR. Analysis of variance revealed that SLC27A6 K81M had a lower p-value (p = .0009) than the most significant SNP in GWAS (p = .0049). Although only SLC27A6 K81M was verified in the present study, subsequent verification of the remaining candidate genes and polymorphisms could lead to the identification of genes and polymorphisms responsible for FAR.

Significance of Angiogenic Growth Factors in Bovine Ocular Squamous Cell Carcinoma

Bovine ocular squamous cell carcinoma (BOSCC) is the most common and economically significant neoplasm of the eye in cattle. This study investigated the role of angiogenic growth factors in the pathogenesis of BOSCC. Eighteen cases of BOSCC were classified histopathologically according to the degree of differentiation. Normal upper and lower eyelids and third eyelids collected from the right and left eyes of six healthy cattle aged 1-3 years, that had been presented for slaughter to abattoirs, served as controls. Transcription of genes encoding the angiogenic growth factors vascular endothelial growth factor-C (VEGF-C), basic fibroblast growth factor (bFGF), platelet-derived growth factor-C (PDGF-C) and platelet-derived growth factor receptor-α (PDGFR-α) was determined by quantitative real-time polymerase chain reaction (RT-PCR) in tissue obtained from paraffin wax blocks. Immunohistochemistry (IHC) was utilized to detect intensity of expression and tissue distribution of these growth factors. IHC results revealed that bFGF and PDGF-C were elevated significantly (P >0.05) and VEGF-C expression was decreased in BOSCC compared with healthy control tissue. PDGR-α expression was elevated; however, the difference, compared with control tissues, was not significant. RT-PCR results showed an inverse relationship to the results of IHC; where protein levels were elevated their corresponding mRNA levels were decreased or vice-versa. Angiogenic regulators therefore appear to play a role in the pathogenesis of BOSCC.

Differential expression of mRNA isoforms in the skeletal muscle of pigs with distinct growth and fatness profiles

Background

The identification of genes differentially expressed in the skeletal muscle of pigs displaying distinct growth and fatness profiles might contribute to identify the genetic factors that influence the phenotypic variation of such traits. So far, the majority of porcine transcriptomic studies have investigated differences in gene expression at a global scale rather than at the mRNA isoform level. In the current work, we have investigated the differential expression of mRNA isoforms in the gluteus medius (GM) muscle of 52 Duroc HIGH (increased backfat thickness, intramuscular fat and saturated and monounsaturated fatty acids contents) and LOW pigs (opposite phenotype, with an increased polyunsaturated fatty acids content).

Results

Our analysis revealed that 10.9% of genes expressed in the GM muscle generate alternative mRNA isoforms, with an average of 2.9 transcripts per gene. By using two different pipelines, one based on the CLC Genomics Workbench and another one on the STAR, RSEM and DESeq2 softwares, we have identified 10 mRNA isoforms that both pipelines categorize as differentially expressed in HIGH vs LOW pigs (P-value < 0.01 and ±0.6 log₂fold-change). Only five mRNA isoforms, produced by the ITGA5, SEMA4D, LITAF, TIMP1 and ANXA2 genes, remain significant after correction for multiple testing (q-value < 0.05 and ±0.6 log₂fold-change), being upregulated in HIGH pigs.

Conclusions

The increased levels of specific ITGA5, LITAF, TIMP1 and ANXA2 mRNA isoforms in HIGH pigs is consistent with reports indicating that the overexpression of these four genes is associated with obesity and metabolic disorders in humans. A broader knowledge about the functional attributes of these mRNA variants would be fundamental to elucidate the consequences of transcript diversity on the determinism of porcine phenotypes of economic interest.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4515-2) contains supplementary material, which is available to authorized users.

Assessment of runs of homozygosity islands and estimates of genomic inbreeding in Gyr (Bos indicus) dairy cattle

Background

Runs of homozygosity (ROH) are continuous homozygous segments of the DNA sequence. They have been applied to quantify individual autozygosity and used as a potential inbreeding measure in livestock species. The aim of the present study was (i) to investigate genome-wide autozygosity to identify and characterize ROH patterns in Gyr dairy cattle genome; (ii) identify ROH islands for gene content and enrichment in segments shared by more than 50% of the samples, and (iii) compare estimates of molecular inbreeding calculated from ROH (F_ROH), genomic relationship matrix approach (F_GRM) and based on the observed versus expected number of homozygous genotypes (F_HOM), and from pedigree-based coefficient (F_PED).

Results

ROH were identified in all animals, with an average number of 55.12 ± 10.37 segments and a mean length of 3.17 Mb. Short segments (ROH_{1–2 Mb}) were abundant through the genomes, which accounted for 60% of all segments identified, even though the proportion of the genome covered by them was relatively small. The findings obtained in this study suggest that on average 7.01% (175.28 Mb) of the genome of this population is autozygous. Overlapping ROH were evident across the genomes and 14 regions were identified with ROH frequencies exceeding 50% of the whole population. Genes associated with lactation (TRAPPC9), milk yield and composition (IRS2 and ANG), and heat adaptation (HSF1, HSPB1, and HSPE1), were identified. Inbreeding coefficients were estimated through the application of F_ROH, F_GRM, F_HOM, and F_PED approaches. F_PED estimates ranged from 0.00 to 0.327 and F_ROH from 0.001 to 0.201. Low to moderate correlations were observed between F_PED-F_ROH and F_GRM-F_ROH, with values ranging from −0.11 to 0.51. Low to high correlations were observed between F_ROH-F_HOM and moderate between F_PED-F_HOM and F_GRM-F_HOM. Correlations between F_ROH from different lengths and F_PED gradually increased with ROH length.

Conclusions

Genes inside ROH islands suggest a strong selection for dairy traits and enrichment for Gyr cattle environmental adaptation. Furthermore, low F_PED-F_ROH correlations for small segments indicate that F_PED estimates are not the most suitable method to capture ancient inbreeding. The existence of a moderate correlation between larger ROH indicates that F_ROH can be used as an alternative to inbreeding estimates in the absence of pedigree records.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4365-3) contains supplementary material, which is available to authorized users.

Genome-wide association for milk production and female fertility traits in Canadian dairy Holstein cattle

BACKGROUND

Genome-wide association studies (GWAS) are a powerful tool for detecting genomic regions explaining variation in phenotype. The objectives of the present study were to identify or refine the positions of genomic regions affecting milk production, milk components and fertility traits in Canadian Holstein cattle, and to use these positions to identify genes and pathways that may influence these traits.

RESULT

Several QTL regions were detected for milk production (MILK), fat production (FAT), protein production (PROT) and fat and protein deviation (FATD, PROTD respectively). The identified QTL regions for production traits (including milk production) support previous findings and some overlap with genes with known relevant biological functions identified in earlier studies such as DGAT1 and CPSF1. A significant region on chromosome 21 overlapping with the gene FAM181A and not previous linked to fertility in dairy cattle was identified for the calving to first service interval and days open. A functional enrichment analysis of the GWAS results yielded GO terms consistent with the specific phenotypes tested, for example GO terms GO:0007595 (lactation) and GO:0043627 (response to estrogen) for milk production (MILK), GO:0051057 (positive regulation of small GTPase mediated signal transduction) for fat production (FAT), GO:0040019 (positive regulation of embryonic development) for first service to calving interval (CTFS) and GO:0043268 (positive regulation of potassium ion transport) for days open (DO). In other cases the connection between the enriched GO terms and the traits were less clear, for example GO:0003279 (cardiac septum development) for FAT and GO:0030903 (notochord development) for DO trait.

CONCLUSION

The chromosomal regions and enriched pathways identified in this study confirm several previous findings and highlight new regions and pathways that may contribute to variation in production or fertility traits in dairy cattle.

Identification of novel isoforms of dairy goat EEF1D and their mRNA expression characterization

Eukaryotic translation elongation factor 1 delta (EEF1D) gene encodes guanine nucleotide exchange protein eEF1Bδ, which participates in the eukaryotic protein synthesis, and plays important roles in regulating cell cycling and milk production. This study firstly focused on detecting the isoforms of dairy goat EEF1D gene and their mRNA expression characterization. Herein, two novel isoforms, EEF1Da and EEF1Dc, were identified in dairy goat. The entire coding sequences of EEF1Da and EEF1Dc isoforms were 843bp and 267bp in length, respectively. Goat EEF1Da had complete conserved domains of elongation factor 1 (EF1) family, and the evolution of goat EEF1Da isoform was agreed with the evolution of species. Expression pattern analysis of different isoforms revealed relatively ubiquitous expression of EEF1D and EEF1Da. While EEF1Dc only expressed in heart, lung, kidney, adipose and muscle. Combining with the analysis results of cloning, qRT-PCR and bioinformatics, EEF1Da is the major alternative splicing form of EEF1D gene. Interestingly, qRT-PCR result showed that the highest expression of EEF1D was in adipose, which is the major component of mammary. This result was consistent with the early research that EEF1D expressed highly in the mammary, which indicated that EEF1D played a potential key role in regulating adipose development and milk production. All these findings would provide a foundation for the further research of EEF1D gene and development of dairy goat industry.

Identification of novel alternative splicing transcript and expression analysis of bovine TMEM95 gene

Transmembrane protein 95 (TMEM95) is closely related to male reproductive performance in cattle, but does not affect semen quality. Alternative splicing plays an important role in regulating biological function as well as in generating proteomic and functional diversity in metazoan organisms. Thus, the aim of this study was to clone and identify transcripts of the TMEM95 gene in cattle using RT-PCR, characterize them via bioinformatics analysis, and detect their expression patterns using qRT-PCR. Two transcripts of TMEM95 were identified in cattle, including TMEM95-SV1 and TMEM95-SV2. Bioinformatics predicted that TMEM95-SV1 has a leucine-rich repeat C-terminal domain and a Pfam: IZUMO. These regions are closely related to protein interactions and the acrosome reaction, respectively. Interestingly, the two transcripts were exclusively expressed in the testes and brain in male fetus cattle, and TMEM95-SV1 was expressed in the brain at significantly higher levels than in the testis (P<0.05, 4.06-fold) and TMEM95-SV2 in the brain (P<0.05, 4.95-fold). These findings enrich the understanding of the TMEM95 gene function and benefit for enhancing male reproduction in cattle industry.

Comparison of significant single nucleotide polymorphisms selections in GWAS for complex traits

The goal of this study was to compare significant SNP selection approaches in the context of complex traits based on SNP estimates obtained by models: a model fitting a single SNP (M1), a model fitting a single SNP and a random polygenic effect (M2), the nonparametric CAR score (M3), a SNP-BLUP model with random effects of all SNPs fitted simultaneously (M4). There were 46,267 SNPs tested in a population of 2601 Holstein Friesian bulls, four traits (milk and fat yields, somatic cell score, non-return rate for heifers) were considered. The numbers of SNPs selected as significant differed among models. M1 selected a very large number of SNPs, except for a NRH in which no SNPs were significant. M2 and M3 both selected similar and low number of SNPs for each trait. M4 selected more SNPs than M2 and M3. Considering linkage disequilibrium between SNPs, for MY M2 and M3 selected SNPs more highly correlated with each other than in the case of M4, while for FY M3 selection contained more correlated SNPs than M2 and M4. In conclusion, if the research interest is to identify SNPs not only with strong, but also with moderate effects on a complex trait a multiple-SNP model is recommended. Such models are capable of accounting for at least a part of linkage disequilibrium between SNPs through the design matrix of SNP effects. Functional annotation of SNPs significant in M4 reveals good correspondence between selected polymorphisms and functional information as well as with QTL mapping results.