Meta-analysis of RNA-Seq datasets highlights novel genes/pathways involved in fat deposition in fat-tail of sheep.

Introduction: Fat-tail in sheep is considered as an important energy reservoir to provide energy as a survival buffer during harsh challenges. However, fat-tail is losing its importance in modern sheep industry systems and thin-tailed breeds are more desirable. Using comparative transcriptome analysis to compare fat-tail tissue between fat- and thin-tailed sheep breeds provides a valuable approach to study the complex genetic factors associated with fat-tail development. However, transcriptomic studies often suffer from issues with reproducibility, which can be improved by integrating multiple studies based on a meta-analysis. Methods: Hence, for the first time, an RNA-Seq meta-analysis on sheep fat-tail transcriptomes was performed using six publicly available datasets. Results and discussion: A total of 500 genes (221 up-regulated, 279 down-regulated) were identified as differentially expressed genes (DEGs). A jackknife sensitivity analysis confirmed the robustness of the DEGs. Moreover, QTL and functional enrichment analysis reinforced the importance of the DEGs in the underlying molecular mechanisms of fat deposition. Protein-protein interactions (PPIs) network analysis revealed the functional interactions among the DEGs and the subsequent sub-network analysis led to identify six functional sub-networks. According to the results of the network analysis, down-regulated DEGs in green and pink sub-networks (like collagen subunits IV, V, and VI, integrins 1 and 2, SCD, SCD5, ELOVL6, ACLY, SLC27A2, and LPIN1) may impair lipolysis or fatty acid oxidation and cause fat accumulation in tail. On the other hand, up-regulated DEGs, especially those are presented in green and pink sub-networks (like IL6, RBP4, LEPR, PAI-1, EPHX1, HSD11B1, and FMO2), might contribute to a network controlling fat accumulation in the tail of sheep breed through mediating adipogenesis and fatty acid biosynthesis. Our results highlighted a set of known and novel genes/pathways associated with fat-tail development, which could improve the understanding of molecular mechanisms behind fat deposition in sheep fat-tail.

Investigation of the transability of dietary small non-coding RNAs to animals.

Our daily diet not only provides essential nutrients needed for survival and growth but also supplies bioactive ingredients to promote health and prevent disease. Recent studies have shown that exogenous microRNAs (miRNAs), xenomiRs, may enter the consumer's body through dietary intake and regulate gene expression. This fascinating phenomenon suggests that xenomiRs can act as a new class of bioactive substances associated with mammalian systems. In contrast, several studies have failed to detect xenomiRs in consumers and reported that the observed diet-derived miRNAs in the previous studies can be related to the false positive effects of experiments. This discrepancy can be attributed to the potential artifacts related to the process of experiments, small sample size, and inefficient bioinformatics pipeline. Since this hypothesis is not generally accepted yet, more studies are required. Here, a stringent and reliable bioinformatics pipeline was used to analyze 133 miRNA sequencing data from seven different studies to investigate this phenomenon. Generally, our results do not support the transfer of diet-derived miRNAs into the animal/human tissues in every situation. Briefly, xenomiRs were absent from most samples, and also, their expressions were very low in the samples where they were present, which is unlikely to be sufficient to regulate cell transcripts. Furthermore, this study showed that the possibility of miRNAs being absorbed through animals' diets and thus influencing gene expression during specific periods of biological development is not inconceivable. In this context, our results were in agreement with the theory of the transfer of small RNAs under certain conditions and periods as xenomiRs were found in colostrum which may modulate infants' immune systems via post-transcriptional regulation. These findings provide evidence for the selective absorption of diet-derived small RNAs, which need to be investigated in future studies to shed light on the mechanisms underlying the transference of diet-derived miRNAs.

Transcriptome analysis reveals the potential roles of long non-coding RNAs in feed efficiency of chicken.

Feed efficiency is an important economic trait and reduces the production costs per unit of animal product. Up to now, few studies have conducted transcriptome profiling of liver tissue in feed efficiency-divergent chickens (Ross vs native breeds). Also, molecular mechanisms contributing to differences in feed efficiency are not fully understood, especially in terms of long non-coding RNAs (lncRNAs). Hence, transcriptome profiles of liver tissue in commercial and native chicken breeds were analyzed. RNA-Seq data along with bioinformatics approaches were applied and a series of lncRNAs and target genes were identified. Furthermore, protein-protein interaction network construction, co-expression analysis, co-localization analysis of QTLs and functional enrichment analysis were used to functionally annotate the identified lncRNAs. In total, 2,290 lncRNAs were found (including 1,110 annotated, 593 known and 587 novel), of which 53 (including 39 known and 14 novel), were identified as differentially expressed genes between two breeds. The expression profile of lncRNAs was validated by RT-qPCR. The identified novel lncRNAs showed a number of characteristics similar to those of known lncRNAs. Target prediction analysis showed that these lncRNAs have the potential to act in cis or trans mode. Functional enrichment analysis of the predicted target genes revealed that they might affect the differences in feed efficiency of chicken by modulating genes associated with lipid metabolism, carbohydrate metabolism, growth, energy homeostasis and glucose metabolism. Some gene members of significant modules in the constructed co-expression networks were reported as important genes related to feed efficiency. Co-localization analysis of QTLs related to feed efficiency and the identified lncRNAs suggested several candidates to be involved in residual feed intake. The findings of this study provided valuable resources to further clarify the genetic basis of regulation of feed efficiency in chicken from the perspective of lncRNAs.

An integrated analysis of mRNAs, lncRNAs, and miRNAs based on weighted gene co-expression network analysis involved in bovine endometritis.

In dairy cattle, endometritis is a severe infectious disease that occurs following parturition. It is clear that genetic factors are involved in the etiology of endometritis, however, the molecular pathogenesis of endometritis is not entirely understood. In this study, a system biology approach was used to better understand the molecular mechanisms underlying the development of endometritis. Forty transcriptomic datasets comprising of 20 RNA-Seq (GSE66825) and 20 miRNA-Seq (GSE66826) were obtained from the GEO database. Next, the co-expressed modules were constructed based on RNA-Seq (Rb-modules) and miRNA-Seq (mb-modules) data, separately, using a weighted gene co-expression network analysis (WGCNA) approach. Preservation analysis was used to find the non-preserved Rb-modules in endometritis samples. Afterward, the non-preserved Rb-modules were assigned to the mb-modules to construct the integrated regulatory networks. Just highly connected genes (hubs) in the networks were considered and functional enrichment analysis was used to identify the biological pathways associated with the development of the disease. Furthermore, additional bioinformatic analysis including protein-protein interactions network and miRNA target prediction were applied to enhance the reliability of the results. Thirty-five Rb-modules and 10 mb-modules were identified and 19 and 10 modules were non-preserved, respectively, which were enriched in biological pathways related to endometritis like inflammation and ciliogenesis. Two non-preserved Rb-modules were significantly assigned to three mb-modules and three and two important sub-networks in the Rb-modules were identified, respectively, including important mRNAs, lncRNAs and miRNAs genes like IRAK1, CASP3, CCDC40, CCDC39, ZMYND10, FOXJ1, TLR4, IL10, STAT3, FN1, AKT1, CD68, ENSBTAG00000049936, ENSBTAG00000050527, ENSBTAG00000051242, ENSBTAG00000049287, bta-miR-449, bta-miR-484, bta-miR-149, bta-miR-30b and bta-miR-423. The potential roles of these genes have been previously demonstrated in endometritis or related pathways, which reinforced putative functions of the suggested integrated regulatory networks in the endometritis pathogenesis. These findings may help further elucidate the underlying mechanisms of bovine endometritis.

Maternal supplementation with fish oil modulates inflammation-related MicroRNAs and genes in suckling lambs.

Dietary n-3 long-chain fatty acids (n-3 LCFA) have been shown to modify lipid metabolism and immune function. The objective of this study was to evaluate the effect of periparturient fish oil (FO) supplementation on the inflammation and metabolic health of ewes and their lambs at a molecular level. Prepartum ewes were fed control diet (CON, n = 12) or CON supplemented with 2% DM of calcium soap of FO (n = 12) from 28 days before until 21 days after parturition. The ewes were evaluated for plasma metabolites and milk composition. The experiment was followed by analyzing the relative transcript abundance of circulating microRNAs (miRNAs) in plasma and targeted miRNA/mRNA expression in peripheral blood mononuclear cells (PBMCs) in both ewes and lambs. FO treatment decreased prepartum feed intake (1812 ± 35 vs 1674 ± 33 g/day, P < 0.01), whereas the influence on plasma metabolites was negligible. Dietary FO supplementation decreased milk fat percentage (8.82 ± 0.49 vs 7.03 ± 0.45, P = 0.02) and reduced milk n-6/n-3 (P < 0.05). Also, it altered the expression of plasma-circulating miRNAs in both ewe and lamb (P < 0.05). Furthermore, maternal nutrition of FO downregulated the relative expression of miR-33a and miR-146b and transcript abundance of genes IL-1ß (0.41-fold) and NF-κB (0.25-fold) in lambs' PBMC. In conclusion, results showed that FO supplementation starting antepartum affects milk composition and circulating miRNA in dams and the inflammatory markers in lambs delivered by the supplemented ewes. These may provide a strategy to maintain immune balance during gestation and develop the immune system in lambs.

Deep transcriptome analysis using RNA-Seq suggests novel insights into molecular aspects of fat-tail metabolism in sheep.

Fat-tail content of sheep breeds is varied and the molecular mechanisms regulating fat-tail development have not been well characterized. Aiming at better identifying the important candidate genes and their functional pathways contributing to fat deposition in the tail, a comparative transcriptome analysis was performed between fat- (Lori-Bakhtiari) and thin-tailed (Zel) Iranian sheep breeds using RNA-seq. The experiment was conducted on six male lambs (three lambs per each breed) at seven months of age. Four different combinations of aligners and statistical methods including Hisat2 + edgeR, Hisat2 + DESeq2, STAR + edgeR and STAR + DESeq2 were used to identify the differentially expressed genes (DEGs). The DEGs were selected for functional enrichment analysis and protein-protein interaction (PPI) network construction. Module analysis was also conducted to mine the functional sub-networks from the PPI network. In total, 264 genes including 80 up- and 184 down-regulated genes were identified as DEGs. The RNA-Seq results were validated by Q-RT-PCR. Functional analysis of DEGs and the module analysis of PPI network demonstrated that in addition to pathways affecting lipid metabolism, a series of enriched functional terms related to "response to interleukin", "MAPK signaling pathways", "Wnt signaling pathway", "ECM-receptor interaction", "regulation of actin cytoskeleton", and "response to cAMP" might contribute to the deposition of fat in tails of sheep. Overall results using RNA-Seq analysis characterized important candidate genes involved in the fatty acid metabolism and regulation of fat deposition, suggesting novel insights into molecular aspects of fat-tail metabolism in sheep. Selected DEGs should be further investigated as potential markers associated with the fat-tail development in sheep breeds.

Genome-wide identification and analysis of A-to-I RNA editing events in bovine by transcriptome sequencing.

RNA editing increases the diversity of the transcriptome and proteome. Adenosine-to-inosine (A-to-I) editing is the predominant type of RNA editing in mammals and it is catalyzed by the adenosine deaminases acting on RNA (ADARs) family. Here, we used a largescale computational analysis of transcriptomic data from brain, heart, colon, lung, spleen, kidney, testes, skeletal muscle and liver, from three adult animals in order to identify RNA editing sites in bovine. We developed a computational pipeline and used a rigorous strategy to identify novel editing sites from RNA-Seq data in the absence of corresponding DNA sequence information. Our methods take into account sequencing errors, mapping bias, as well as biological replication to reduce the probability of obtaining a false-positive result. We conducted a detailed characterization of sequence and structural features related to novel candidate sites and found 1,600 novel canonical A-to-I editing sites in the nine bovine tissues analyzed. Results show that these sites 1) occur frequently in clusters and short interspersed nuclear elements (SINE) repeats, 2) have a preference for guanines depletion/enrichment in the flanking 5'/3' nucleotide, 3) occur less often in coding sequences than other regions of the genome, and 4) have low evolutionary conservation. Further, we found that a positive correlation exists between expression of ADAR family members and tissue-specific RNA editing. Most of the genes with predicted A-to-I editing in each tissue were significantly enriched in biological terms relevant to the function of the corresponding tissue. Lastly, the results highlight the importance of the RNA editome in nervous system regulation. The present study extends the list of RNA editing sites in bovine and provides pipelines that may be used to investigate the editome in other organisms.

Bioactive glass ceramic nanoparticles-coated poly(l-lactic acid) scaffold improved osteogenic differentiation of adipose stem cells in equine.

Horses with big bone fractures have low chance to live mainly due to the lake of a proper treatment strategy. We believe that further attempts in equine bone tissue engineering will probably be required to meet all the needs for the lesion therapies. Therefore in this study we aimed to investigate the osteogenic differentiation capacity of equine adipose-derived stem cells (e-ASCs) on nano-bioactive glass (nBGs) coated poly(l-lactic acid) (PLLA) nanofibers scaffold (nBG-PLLA). Using electrospinning technique, PLLA scaffold was prepared successfully and coated with nBGs. Fabricated nanofibers were characterized by MTT, SEM, and FTIR analyses, and then osteogenic differentiation potential of isolated e-ASCs was investigated by the most key osteogenic markers, namely Alizarin red-S, ALP, calcium content and bone related (RUNX2, Collagen I, Osteonectin, and ALP) gene markers. Our results indicated that nBGs was successfully coated on PLLA scaffold and this scaffold had no negative (p>0.05) effect on cell growth rate as indicated by MTT assay. Moreover, e-ASCs that differentiated on nBGs-PLLA scaffold showed a higher (p<0.05) ALP activity, more (p<0.05) calcium content, and higher (p<0.05) expression of bone-related genes than that on uncoated PLLA scaffold and TCPS. According to the results, a combination of bioceramics and biopolymeric nanofibers hold valuable promising potentials to use for bone tissue engineering application and regenerative medicine.

The Association of Bovine Osteopontin (OPN) Gene with Milk Production Traits in Iranian Holstein Bulls.

BACKGROUND: The Osteopontin (OPN) is a highly phosphorylated glycoprotein in numbers of bovine tissues and milk. OPN has been reported to be associated with milk production in cattle. OBJECTIVE: The genotype and allelic frequencies for OPN and its association with milk production will be evaluated in Iranian Holstein Bulls. MATERIALS AND METHODS: Bulls DNA (100) was isolated. Oligo was used for primer design. Polymerase Chain Reaction was implemented to amplify a 826 bp fragment and the amplicon was digested by BsrI. Restricted Maximum likelihood (REML) method based on average information algorithm using ASRMEL programs (version 3.1) was employed to estimate the genetic parameters and variance of components. The association of OPN genotypes with milk production traits were analysed by the least square method as applied in the general linear model (GLM) procedure of SAS. Allele substitution effects were performed by regression analyses. RESULTS: Allele frequencies of T and C were 0.59±0.03 and 0.41±0.03, respectively. Genotype frequencies of TT, CT and CC were 34.69, 48.62, and 16.69, respectively. The chi-square test showed the deviation from Hardy-Weinberg equilibrium. Estimated heritability for milk yield, fat yield and its percent, protein yield and its percent were 0.28±0.0061, 0.21±0.0064, 0.22±0.0086, 0.32±0.0065 and 0.34±0.0096 respectively. Allelic substitution effects and differences between genotypes were not significant for milk production traits. CONCLUSIONS: This study suggested that the C allele frequency of OPN was noticeable in Iranian proven bull Holstein population, but was not associated with milk production traits. However, before being practical for the breeding improvement of Iranian Holsteins a larger sample size is required.

Single nucleotide polymorphism of FGF2 gene in Iranian Holstein proven bulls.

Fibroblast growth factor 2 (FGF2) serves in the uterine endometrium during estrous presenting in the bovine mammary gland which is responsible to express interferon-T (IFNT), and is an important agent to encourage the continuation of pregnancy in the ruminants. Significant associations have been found between genes affected by IFNT and genes that are responsible for milk production traits. Semen samples from 101Iranian Holstein proven bulls were collected to extract the genomic DNA. Forward and reverse primers were designed and a 710-base-pair fragment in intron 1 was amplified using PCR technique. To detect single nucleotide polymorphism (SNP), all samples were sequenced. Three positions including 11474 (C/G), 11513 (C/G) and 11646 (A/G) were considered. The 11474C, 11513C and 11646A alleles are known as wild type alleles. In this study all animals were distinguished as the11474C, 11513C and 11646A alleles. Furthermore, amplified fragments were under consideration to detect new SNPs. Only one new SNP in one sample was observed at position 11863 resulting substitution of thymine to cytosine. This new mutation has been registered on the NCBI database with accession number HM597774.

Application of myostatin in sheep breeding programs: A review.

Plasma membrane H+-ATPase is a major integral membrane protein with a role in various physiological processes including abiotic stress response. To study the effect of NaCl on the expression pattern of a gene encoding the plasma membrane H+-ATPase, an experiment was carried out in a completely random design with three replications. A pair of specific primers was designed based on the sequence of the gene encoding plasma membrane H+-ATPase in Aeluropus littoralis to amplify a 259 bp fragment from the target gene by PCR. A gene encoding actin was used as reference gene to normalize the expression level of the target gene. A pair of specific primers was designed to amplify a 157 bp fragment from the actin gene by PCR. Plants were treated with different concentrations of NaCl, 0, 50, 100, 150, 200, 250, 500 and 1000 mM, for two days. Our results showed that the expression level of the plasma membrane H+-ATPase gene increased dramatically at 500 mM and then decreased with increasing concentrations of NaCl. The results also indicated that the leaves of plants, were treated with high concentrations of NaCl changed morphologically, but those grown under low concentrations of NaCl as well as the control plants did not show morphological changes in their leaves. Our results suggest a relation between morphological changes of treated plants and the expression level of the plasma membrane H+-ATPase gene in Aeluropus littoralis.

Evaluation of two DNA extraction methods from maternal plasma for using in non-invasive bovine fetus gender determination.

BACKGROUND: Fetal DNA in maternal plasma and serum has been shown to be a useful material for prenatal fetal sex determination during early gestational ages. Non-invasive prenatal diagnosis is now possible at 8(th) week of pregnancy, by maternal blood sample testing. OBJECTIVE: The purpose of this study was to evaluate two DNA extraction methods from mother plasma and its routine clinical application in bovine fetus gender determination with non-invasive method. MATERIALS AND METHODS: Maternal blood samples were taken from 40 pregnant cows during the 8(th)-38(th) weeks of gestation. DNA was extracted from 350 µl of maternal plasma with two salting-out and phenol-chloroform methods. The absorption in A260 and purity (A260/A280) of extracted DNA were detected by ultraviolet spectrophotometer. Three µl of the extracted DNA with phenol-chloroform method was used as a template. The PCR reaction was carried out to amplify the fragments of X and Y chromosomes of amelogenin, TSPY and BC1.2 genes. RESULTS: The difference between the mean absorption of DNA extracted by phenol-chloroform method and salting-out method was not significant in A260 (p>0.05, p=0.3549), but the difference between mean purity (A260/A280) of DNA extracted by phenol-chloroform method and salting-out method was significant (p<0.001). X chromosome fragment was detected in all 40 samples and Y chromosome fragments were detected in 25 plasma samples which were delivered a male calf. The sensitivity and specificity of test was 100% with no false negative and false positive results. CONCLUSION: The results showed that phenol-chloroform method is a simple and sensitive method for isolation of fetal DNA in maternal plasma.

A low-cost efficient multiplex PCR for prenatal sex determination in bovine fetus using free fetal DNA in maternal plasma.

BACKGROUND: In order to establish a reliable non-invasive method for sex determination in a bovine fetus in a routine setting, the possibility of identifying specific sequence in the fetal X and Y-chromosomes has been evaluated in maternal plasma using conventional multiplex polymerase chain reaction (PCR) analysis. The aim of this study was to provide a rapid and reliable method for sexing bovine fetuses. MATERIALS AND METHODS: In this experimental study, peripheral blood samples were taken from 38 pregnant heifers with 8 to 38 weeks of gestation. DNA template was extracted by phenol-chloroform method from 350 µl maternal plasma. Two primer pairs for bovine amelogenin gene (bAML) and BC1.2 were used to amplify fragments from X and Y chromosomes. A multiplex PCR reaction has been optimized for amplification of 467 bp and 341 bp fragments from X and Y bAML gene and a 190 bp fragment from BC1.2 related to Y chromosome. RESULTS: The 467 bp fragment was observed in all 38 samples. Both 341 and 190 bp fragments were detected only in 24 plasma samples from male calves. The sensitivity and specificity of test were 100% with no false negative or false positive results. CONCLUSION: The results showed that phenol-chloroform method is a simple and suitable method for isolation of fetal DNA in maternal plasma. The multiplex PCR method is an available non-invasive approach which is cost efficient and reliable for sexing bovine fetuses.

PCR based RFLP genotyping of bovine lymphocyte antigen DRB3.2 in Iranian Holstein population.

Major Histocompatibility Complex (MHC) class II locus DRB3 was investigated by PCR based restriction fragment length polymorphism (PCR-RFLP) assay. A total of 262 Holstein cows participating in the national recording system were sampled from 10 herds. A two-step polymerase chain reaction was carried out in order to amplify a 284 base-pair fragment of exon 2 of the target gene. Second PCR products were treated with three restriction endonucleas enzymes RsaI, BstYI and HaeIII. Digested fragments were analyzed by polyacrylamid gel electrophoresis. Twenty-eight BoLA-DRB3 alleles were identified. Identified alleles are: BoLA-DRB3.2 *3, *6, *7, *8, *9, *10, *11, *12, *13, *14, *15, *16,20, *21, *22, *23, *24, *25, *26, *27, *28, *32, *36, *37, *40, *51, *iaa and *ibb. The BoLA-DRB3.2*40 allele that was observed in this study has not been reported previously. The calculated frequencies were as follows: 2.29, 1.34, 0.19, 14.5, 0.38, 3.05, 12.21, 1.34, 2.29, 1.34, 2.48, 9.16, 0.95, 0.77, 6.68, 9.16, 17.94, 1.15, 0.57, 1.15, 0.95, 0.57, 0.38, 1.91, 0.38, 5.73, 0.19 and 0.95% respectively. The six most frequently observed alleles (BoLA-DRB3.2 *8, *11, *16, *22, *23 and *24) accounted for 69.65% of the alleles in these 10 herds. The results of this study confirm the allelic distribution of six most frequent alleles in Holstein population's worldwide.