In silico identification of conserved microRNAs and their targets in bovine fat tissue

Transcriptome analysis of adipose tissue in grazing cattle: Identifying key regulators of fat metabolism

The taste and tenderness of meat are the main determinants of carcass quality in many countries. This study aimed to discuss the mechanisms of intramuscular fat deposition in grazing and house-breeding cattle. We performed transcriptome analysis to characterize messenger RNA and microRNA (miRNA) expression profiles. A total of 456 and 66 differentially expressed genes (DEGs) and differentially expressed (DE) miRNAs were identified in the adipose tissue of grazing and house-breeding cattle. Kyoto Encyclopedia of Genes and Genomes pathway analysis identified the association of DEGs with fatty acid metabolism, fatty acid degradation, peroxisome proliferator-activated receptors signaling pathway, adenosine monophosphate-activated protein kinase signaling pathway, adipocytokine signaling pathway, and the association of DE miRNAs with mitogen-activated protein kinase signaling pathway. Apolipoprotein L domain containing 1, pyruvate dehydrogenase kinase 4, and sphingosine-1-phosphate lyase 1 genes may be the key regulators of fat metabolism in grazing cattle. Finally, we found that miR-211 and miR-331-5p were negatively correlated with the elongation of very long-chain fatty acids protein 6 (ELOVL6), and miR-331-5p might be the new regulator involved in fat metabolism. The results indicated that ELOVL6 participated in various functions and pathways related to fat metabolism. Meanwhile, miR-331-5p, as a new regulator, might play an essential role in this process. Our findings laid a more in-depth and systematic research foundation for the formation mechanism and characteristics of adipose tissue in grazing cattle.

Identification and in silico characterization of structural and functional impacts of genetic variants in milk protein genes in the Zebu breeds Guzerat and Gyr

Whole genome sequencing of bovine breeds has allowed identification of genetic variants in milk protein genes. However, functional repercussion of such variants at a molecular level has seldom been investigated. Here, the results of a multistep Bioinformatic analysis for functional characterization of recently identified genetic variants in Brazilian Gyr and Guzerat breeds is described, including predicted effects on the following: (i) evolutionary conserved nucleotide positions/regions; (ii) protein function, stability, and interactions; (iii) splicing, branching, and miRNA binding sites; (iv) promoters and transcription factor binding sites; and (v) collocation with QTL. Seventy-one genetic variants were identified in the caseins (CSN1S1, CSN2, CSN1S2, and CSN3), LALBA, LGB, and LTF genes. Eleven potentially regulatory variants and two missense mutations were identified. LALBA Ile60Val was predicted to affect protein stability and flexibility, by reducing the number the disulfide bonds established. LTF Thr546Asn is predicted to generate steric clashes, which could mildly affect iron coordination. In addition, LALBA Ile60Val and LTF Thr546Asn affect exonic splicing enhancers and silencers. Consequently, both mutations have the potential of affecting immune response at individual level, not only in the mammary gland. Although laborious, this multistep procedure for classifying variants allowed the identification of potentially functional variants for milk protein genes.

In silico identification of variations in microRNAs with a potential impact on dairy traits using whole ruminant genome SNP datasets

MicroRNAs are small noncoding RNAs that have important roles in the lactation process and milk biosynthesis. Some polymorphisms have been studied in various livestock species from the perspective of pathology or production traits. To target variants that could be the causal variants of dairy traits, genetic variants of microRNAs expressed in the mammary gland or present in milk and localized in dairy quantitative trait loci (QTLs) were investigated in bovine, caprine, and ovine species. In this study, a total of 59,124 (out of 28 millions), 13,427 (out of 87 millions), and 4761 (out of 38 millions) genetic variants in microRNAs expressed in the mammary gland or present in milk were identified in bovine, caprine, and ovine species, respectively. A total of 4679 of these detected bovine genetic variants are located in dairy QTLs. In caprine species, 127 genetic variants are localized in dairy QTLs. In ovine species, no genetic variant was identified in dairy QTLs. This study leads to the detection of microRNA genetic variants of interest in the context of dairy production, taking advantage of whole genome data to identify microRNA genetic variants expressed in the mammary gland and localized in dairy QTLs.

Differential expression of miRNAs and related mRNAs during heat stress in buffalo heifers

The present experiment was aimed to study differential expression of miRNAs and related mRNAs during heat stress (HS) in buffalo heifers. Twelve Murrah buffalo heifers aged between 1.5 and 2.0 years, weighting between 250 and 300 Kg were randomly assigned into two equal groups. The animals were kept in the psychrometric chamber under Thermo-neutral (TN; THI = 72) and HS (THI = 87-90) conditions for 6 h every day between 1000 and 1600 h for 21 days. The blood sampling was done at 1500 h on 15th day of the experiment and physiological parameters viz. pulse rate (PR), respiratory rate (RR) and rectal temperature (RT) were recorded at 1500 h on day -5, -3, -1, 0, +1, +3, +5 with respect to blood sampling. PBMCs were used for extraction of miRNAs and total RNA; and first strand cDNA was synthesized. qPCR was performed for relative gene expression studies. Physiological, hematological (erythrocytic indices), biochemical (triglycerides, urea, ALT, AST, LDH), redox (SOD, ROS) and endocrine parameters (T₄) altered significantly (P < 0.05) during HS as compared to TN. Out of eight targeted miRNAs only four were expressed in buffalo heifers. The relative expression of bta-mir-142, bta-mir-1248 and bta-mir-2332 was significantly (P < 0.05) up-regulated whereas expression of bta-mir-2478 was significantly (P < 0.05) down-regulated during HS as compared to TN. The relative expression of the predicted target genes i.e. HSF1, HSP60, HSP70, HSPA8 and HSP90 were significantly (P < 0.05) up-regulated whereas HSF4 expression was significantly (P < 0.05) down-regulated during HS as compared to TN. It can be concluded that a THI of 87-90 could lead to a moderate HS in buffalo heifers. Differential expression studies of miRNAs and related mRNAs in present study deciphers the role of miRNAs in the heat tolerance in buffalo heifers.

In silico identification of conserved miRNAs and their selective target gene prediction in indicine (Bos indicus) cattle

The modern cattle was domesticated from aurochs, sharing its physiological traits into two subspecies Bos taurus and Bos indicus. MicroRNAs (miRNAs) are a class of non-coding short RNAs of ~22nt which have a key role in the regulation of many cellular and physiological processes in the animal. The current study was aimed to predict and annotate the potential mutations in indicine miRNAs throughout the genome using de novo and homology-based in silico approaches. Genome-wide mapping was performed in available indicine assembly by the homology-based approach and 768 miRNAs were recovered out of 808 reported taurine miRNAs belonging to 521 unique mature miRNA families. While 42 precursors were dropped due to lack of secondary miRNA structure, increasing stringency or decreasing similarity between the two genomes' miRNA. Increasing tendency of miRNAs incidence was observed on chr5, chr7, chr8, chr12 and chr21 with 19 polycistronic miRNA within 1-kilobase distance throughout the indicine genome. Notably, 12 miRNAs showed copy number variation. Eighteen miRNAs showed a mutation in their mature sequences in which eight were found in their seed region. Whilst in de novo based approach, 12 novel potential miRNAs on Y chromosome in indicine cattle along with a new miRNA (bind-miR-1264) on chrX were found. The final data set is annotated and explains the impending target genes that are responsible for enhanced immunity, heat tolerance and disease tolerance regulation in indicine. The study conforms to better understanding and perceptive approach towards indicine genome.

TRIENNIAL GROWTH AND DEVELOPMENT SYMPOSIUM: Molecular mechanisms related to bovine intramuscular fat deposition in the longissimus muscle

The intramuscular fat (IMF) content of the LM, also known as marbling, is particularly important in determining the price of beef in Korea, Japan, and the United States. Deposition of IMF is influenced by both genetic (e.g., breed, gender, and genotype) and nongenetic factors (e.g., castration, nutrition, stressors, animal weight, and age). Castration of bulls markedly increases deposition of IMF, resulting in improved beef quality. Here, we present a comparative gene expression approach between bulls and steers. Transcriptomic and proteomic studies have demonstrated that the combined effects of increases in lipogenesis, fatty acid uptake, and fatty acid esterification and decreased lipolysis are associated with increased IMF deposition in the LM. Several peripheral tissues (LM, adipose tissues, and the liver) are involved in lipid metabolism. Therefore, understanding the significance of the tissue network in lipid metabolism is important. Here, we demonstrate that lipid metabolism in LM tissues is crucial for IMF deposition, whereas lipid metabolism in the liver plays only a minor role. Metabolism of body fat and IMF deposition in bovine species has similarities with these processes in metabolic diseases, such as obesity in humans and rodents. Extensive studies on metabolic diseases using epigenome modification (DNA methylation, histone modification, and microRNA), microbial metagenomics, and metabolomics have been performed in humans and rodents, and new findings have been reported using these technologies. The importance of applying "omics" fields (epigenomics, metagenomics, and metabolomics) to the study of IMF deposition in cattle is described. New information on the molecular mechanisms of IMF deposition may be used to design nutritional or genetic methods to manipulate IMF deposition and to modify fatty acid composition in beef cattle. Applying nutrigenomics could maximize the expression of genetic potential of economically important traits (e.g., marbling) in animals.

bta-miR-23a involves in adipogenesis of progenitor cells derived from fetal bovine skeletal muscle

Intramuscular fat deposition or marbling is essential for high quality beef. The molecular mechanism of adipogenesis in skeletal muscle remains largely unknown. In this study, we isolated Platelet-derived growth factor receptor α (PDGFRα) positive progenitor cells from fetal bovine skeletal muscle and induced into adipocytes. Using miRNAome sequencing, we revealed that bta-miR-23a was an adipogenic miRNA mediating bovine adipogenesis in skeletal muscle. The expression of bta-miR-23a was down-regulated during differentiation of PDGFRα⁺ progenitor cells. Forced expression of bta-miR-23a mimics reduced lipid accumulation and inhibited the key adipogenic transcription factor peroxisome proliferative activated receptor gamma (PPARγ) and CCAAT/enhancer binding protein alpha (C/EBPα). Whereas down-regulation of bta-miR-23a by its inhibitors increased lipid accumulation and expression of C/EBPα, PPARγ and fatty acid-binding protein 4 (FABP4). Target prediction analysis revealed that ZNF423 was a potential target of bta-miR-23a. Dual-luciferase reporter assay revealed that bta-miR-23a directly targeted the 3′-UTR of ZNF423. Together, our data showed that bta-miR-23a orchestrates early intramuscular adipogeneic commitment as an anti-adipogenic regulator which acts by targeting ZNF423.

Effect of pre-miRNA-1658 gene polymorphism on chicken growth and carcass traits

Objective

Polymorphisms occurring in the precursor region of microRNAs (miRNAs) affect the target gene and alter the biogenesis of miRNAs, resulting in phenotypic variation. The purpose of the study was to investigate the genetic effects of rs16681031 (C>G) mutation in the precursor region of gga-miR-1658 on the economic traits of the Gushi-Anka chicken F2 resource population.

Methods

To explore the effect of miR-1658 polymorphisms on chicken economic traits, the SNP was genotyped by MassArray matrix-assisted laser desorption/ionization-time of flight mass spectrometry. The association between the SNP and chicken body size, growth and carcass traits was determined by linear mixed models.

Results

The SNP was not only significantly associated with body weight at the age of 6, 8, 10, 12 weeks, respectively, but also with the breadth of the chicken chest, body slanting length and pelvic breadth at 4 weeks, chest depth at 8 weeks of age, and body slanting length at 12 weeks (p<0.05), respectively.

Conclusion

Our data serve as a useful resource for further analysis of miRNA function, and represent a molecular genetic basis for poultry breeding.

Maternal Plane of Nutrition During Late-Gestation and Weaning Age Alter Steer Calf Longissimus Muscle Adipogenic MicroRNA and Target Gene Expression

The main objective was to evaluate if different planes of maternal nutrition during late gestation and weaning age alter microRNA (miRNA) and target gene expression in offspring longissimus muscle (LM). Early (EW) and normal weaned (NW) Angus × Simmental calves (n = 30) born to cows that were grazing endophyte-infected tall fescue and red clover pastures with no supplement [low plane of nutrition (LPN)], or supplemented with 2.3 and 9.1 kg of dried distiller's grains with solubles and soy hulls [medium and high plane of nutrition (MPN, HPN), respectively] during the last 105 ± 11 days of gestation were used. Biopsies of LM were harvested at 78 (early weaning), 187 (normal weaning) and 354 days of age. Results indicate a role of pro-adipogenic miRNA in the control of adipogenesis in LM of NW-MPN steers between 78 and 187 days of age through upregulation of (1) miR-103 which inhibits CAV1, a protein that destabilizes INSR and leads to insulin resistance; (2) miR-143 which inhibits DLK1, a protein that inhibits adipocyte differentiation; and (3) miR-21 which impairs TGFBR2-induced inhibition of adipocyte differentiation. Among the studied anti-adipogenic miRNA, cow plane of nutrition resulted in downregulation of miR-34a expression in MPN steers compared with HPN and LPN at 78 days of age. Data for miR-34a provided a potential sign of epigenetic regulation of LM in beef offspring due to the cow plane of nutrition during late gestation.

Computational Prediction of Alzheimer's and Parkinson's Disease MicroRNAs in Domestic Animals

As the most common neurodegenerative diseases, Alzheimer’s disease (AD) and Parkinson’s disease (PD) are two of the main health concerns for the elderly population. Recently, microRNAs (miRNAs) have been used as biomarkers of infectious, genetic, and metabolic diseases in humans but they have not been well studied in domestic animals. Here we describe a computational biology study in which human AD- and PD-associated miRNAs (ADM and PDM) were utilized to predict orthologous miRNAs in the following domestic animal species: dog, cow, pig, horse, and chicken. In this study, a total of 121 and 70 published human ADM and PDM were identified, respectively. Thirty-seven miRNAs were co-regulated in AD and PD. We identified a total of 105 unrepeated human ADM and PDM that had at least one 100% identical animal homolog, among which 81 and 54 showed 100% sequence identity with 241 and 161 domestic animal miRNAs, respectively. Over 20% of the total mature horse miRNAs (92) showed perfect matches to AD/PD-associated miRNAs. Pigs, dogs, and cows have similar numbers of AD/PD-associated miRNAs (63, 62, and 59). Chickens had the least number of perfect matches (34). Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses suggested that humans and dogs are relatively similar in the functional pathways of the five selected highly conserved miRNAs. Taken together, our study provides the first evidence for better understanding the miRNA-AD/PD associations in domestic animals, and provides guidance to generate domestic animal models of AD/PD to replace the current rodent models.

Bioinformatic identification and expression analysis of Nelumbo nucifera microRNA and their targets

PREMISE OF THE STUDY

Sacred lotus (Nelumbo nucifera) is a perennial aquatic herbaceous plant of ecological, ornamental, and economic importance. MicroRNAs (miRNAs) play an important role in plant development. However, reports of miRNAs and their role in sacred lotus have been limited.

METHODS

Using the homology search of known miRNAs with genome and transcriptome contig sequences, we employed a pipeline to identify miRNAs in N. nucifera. We also predicted the targets of these miRNAs.

RESULTS

We found 106 conserved miRNAs in N. nucifera, and 456 of their miRNA targets were annotated. Quantitative real-time PCR (qRT-PCR) analysis revealed the different expression levels of the 10 selected conserved miRNAs in tissues of young leaves, stems, and flowers of N. nucifera. Negative correlation of expression level between five miRNAs and their target genes was also revealed.

DISCUSSION

Combining bioinformatics and experiment analysis, we identified the miRNAs in N. nucifera. The results can be used as a workbench for further investigation of the roles of miRNAs in N. nucifera.

Dietary MicroRNA Database (DMD): An Archive Database and Analytic Tool for Food-Borne microRNAs

With the advent of high throughput technology, a huge amount of microRNA information has been added to the growing body of knowledge for non-coding RNAs. Here we present the Dietary MicroRNA Databases (DMD), the first repository for archiving and analyzing the published and novel microRNAs discovered in dietary resources. Currently there are fifteen types of dietary species, such as apple, grape, cow milk, and cow fat, included in the database originating from 9 plant and 5 animal species. Annotation for each entry, a mature microRNA indexed as DM0000*, covers information of the mature sequences, genome locations, hairpin structures of parental pre-microRNAs, cross-species sequence comparison, disease relevance, and the experimentally validated gene targets. Furthermore, a few functional analyses including target prediction, pathway enrichment and gene network construction have been integrated into the system, which enable users to generate functional insights through viewing the functional pathways and building protein-protein interaction networks associated with each microRNA. Another unique feature of DMD is that it provides a feature generator where a total of 411 descriptive attributes can be calculated for any given microRNAs based on their sequences and structures. DMD would be particularly useful for research groups studying microRNA regulation from a nutrition point of view. The database can be accessed at http://sbbi.unl.edu/dmd/.