Multi-Omics Integration and Network Analysis Reveal Potential Hub Genes and Genetic Mechanisms Regulating Bovine Mastitis.

Mastitis, inflammation of the mammary gland, is the most prevalent disease in dairy cattle that has a potential impact on profitability and animal welfare. Specifically designed multi-omics studies can be used to prioritize candidate genes and identify biomarkers and the molecular mechanisms underlying mastitis in dairy cattle. Hence, the present study aimed to explore the genetic basis of bovine mastitis by integrating microarray and RNA-Seq data containing healthy and mastitic samples in comparative transcriptome analysis with the results of published genome-wide association studies (GWAS) using a literature mining approach. The integration of different information sources resulted in the identification of 33 common and relevant genes associated with bovine mastitis. Among these, seven genes-CXCR1, HCK, IL1RN, MMP9, S100A9, GRO1, and SOCS3-were identified as the hub genes (highly connected genes) for mastitis susceptibility and resistance, and were subjected to protein-protein interaction (PPI) network and gene regulatory network construction. Gene ontology annotation and enrichment analysis revealed 23, 7, and 4 GO terms related to mastitis in the biological process, molecular function, and cellular component categories, respectively. Moreover, the main metabolic-signalling pathways responsible for the regulation of immune or inflammatory responses were significantly enriched in cytokine-cytokine-receptor interaction, the IL-17 signaling pathway, viral protein interaction with cytokines and cytokine receptors, and the chemokine signaling pathway. Consequently, the identification of these genes, pathways, and their respective functions could contribute to a better understanding of the genetics and mechanisms regulating mastitis and can be considered a starting point for future studies on bovine mastitis.

Diet-derived transmission of MicroRNAs from host plant into honey bee Midgut.

BACKGROUND: MicroRNA (miRNA) is a class of small noncoding RNAs, which targets on thousands of mRNA and thus plays important roles in many biological processes. It has been reported that miRNA has cross-species regulation functions between parasitoid-host, or plant-animal, etc. For example, several plant miRNAs enter into the honey bees and regulate gene expression. However, whether cross-species regulation function of miRNAs is a universal mechanism remains a debate question. RESULTS: We have evaluated transmission of miRNAs from sunflower and sedr plants into the midgut of honey bee using RNA-Seq analyses complemented with confirmation by RT-qPCR. The results showed that at least 11 plant miRNAs were found in the midgut of honey bee feeding by sunflower and sedr pollen. Among which, nine miRNAs, including miR-30d, miR-143, miR-148a, miR-21, let-7 g, miR-26a, miR-126, miR-27a, and miR-203, were shared between the sunflower- and sedr-fed honey bees, suggesting they might have essential roles in plant-insect interactions. Moreover, existence of these co-shared miRNAs presents a strong evidence to support the successful transmission of miRNAs into the midgut of the insect. In total, 121 honeybee mRNAs were predicted to be the target of these 11 plant-derived miRNAs. Interestingly, a sedr-derived miRNA, miR-206, targets on 53 honeybee genes. Kyoto Encyclopedia of Genes and Genome (KEGG) analyses showed that these target genes are significantly involved in hippo signaling pathway-fly, Wnt signaling pathway, and N-Glycan biosynthesis. CONCLUSIONS: In summary, these results provide evidence of cross-species regulation function of miRNA between honeybee and flowering host plants, extending our understanding of the molecular interactions between plants and animals.

ADAR expression and copy number variation in patients with advanced gastric cancer.

BACKGROUND: Gastric cancer (GC) is a world health problem and it is the third leading cause of cancer deaths worldwide. The current practice for prognosis assessment in GC is based on radiological and pathological criteria and they may not result in an accurate prognosis. The aim of this study is to evaluate expression and copy number variation of the ADAR gene in advanced GC and clarify its correlation with survival and histopathological characteristics. METHODS: Forty two patients with stage III and IV GC were included in this study. ADAR gene expression and copy number variation were measured by real-time PCR and Quantitative multiplex fluorescent-PCR, respectively. Survival analysis performed based on the Kaplan-Meier method and Mantel-Cox test. RESULTS: ADAR mRNA was significantly overexpressed in the tumor tissues when compared to the adjacent normal tissues (p < 0.01). Also, ADAR expression level in stage IV was higher than stage III. 40% of patients showed amplification in ADAR gene and there was a positive correlation between ADAR copy number and expression. Increased ADAR expression was clearly correlated with poorer survival outcomes and Mantel-Cox test showed statistically significant differences between low and high expression groups (p < 0.0001). ADAR overexpression and amplification were significantly associated with metastasis, size and stage of tumor. CONCLUSIONS: Together, our data indicate that amplification leads to over expression of ADAR and it could be used as a prognostic biomarker for disease progression, especially for the metastatic process in GC.

Conditions of embryo culture from days 5 to 7 of development alter the DNA methylome of the bovine fetus at day 86 of gestation.

PURPOSE: We tested whether in vitro production (IVP) causes changes in DNA methylation in fetal liver and skeletal muscle and if exposure of cultured embryos to colony-stimulating factor 2 (CSF2) alters DNA methylation. METHODS: Female fetuses were produced by artificial insemination or transfer of an IVP embryo. Embryos were treated from days 5 to 7 after fertilization with CSF2 or vehicle. DNA methylation in fetal liver and skeletal muscle was determined by post-bisulfite adaptor tagging-based sequencing. The degree of DNA methylation for CpG sites in 50-bp windows of the promoter region 500 bp upstream of the transcriptional start site was compared between treatments. RESULTS: For liver, there were 12 genes (6% of those analyzed) in which DNA methylation was affected by treatment, with one 50-bp window per gene affected by treatment. For muscle, the degree of DNA methylation was affected by treatment for 32 windows (19% of the total windows analyzed) representing 28 distinct genes (23% of analyzed genes). For 19 of the 28 genes in muscle, the greatest deviation in DNA methylation was for the CSF2 group. CONCLUSION: Results are consistent with alterations in the methylome being one of the mechanisms by which IVP can result in altered fetal development and postnatal function in the resultant offspring. In addition, results indicate that maternally derived cell-signaling molecules can regulate the pattern of DNA methylation.

Construction of an integrated gene regulatory network link to stress-related immune system in cattle.

The immune system is an important biological system that is negatively impacted by stress. This study constructed an integrated regulatory network to enhance our understanding of the regulatory gene network used in the stress-related immune system. Module inference was used to construct modules of co-expressed genes with bovine leukocyte RNA-Seq data. Transcription factors (TFs) were then assigned to these modules using Lemon-Tree algorithms. In addition, the TFs assigned to each module were confirmed using the promoter analysis and protein-protein interactions data. Therefore, our integrated method identified three TFs which include one TF that is previously known to be involved in immune response (MYBL2) and two TFs (E2F8 and FOXS1) that had not been recognized previously and were identified for the first time in this study as novel regulatory candidates in immune response. This study provides valuable insights on the regulatory programs of genes involved in the stress-related immune system.

OOgenesis_Pred: A sequence-based method for predicting oogenesis proteins by six different modes of Chou's pseudo amino acid composition.

Regarding to critical roles of oogenesis in formation of ova or unfertilized eggs from the oogonia by mitotic division and subsequent differentiation, the identification of oogenesis-related proteins is of great interest. However, the experimental determination of proteins involved in oogenesis is expensive, time consuming and labor-intensive. Therefore, a new powerful discriminating model is indispensable for classifying oogenesis/non-oogenesis-related proteins with high accuracy and precision. Hereby, for the first time we developed a support vector machine based oogenesis protein prediction method which differentiates oogenesis from non-oogenesis proteins. By means of informative protein physicochemical properties and in addition parameter optimization scheme, our method yields a robust and consistent performance. Our model achieved 87.68% and 84.82% prediction accuracy by five-fold cross validation test for datasets with 90% and 50% identity, respectively. The prediction model was also assessed using the independent dataset and yielded 91.62% and 85.38% prediction accuracy for datasets with 90% and 50% identity, respectively, which further demonstrates the effectiveness of our method. Moreover, by applying 10 different feature weighting methods, the more important protein features for oogenesis/non-oogenesis-related proteins discrimination, including serine and glycine frequency, quasi-sequence-order, pseudo-amino acid composition, distribution and conjoint triad, were determined. The success rates revealed that our model can be considered as a new encouraging and strong model for predicting proteins involved in oogenesis with appropriate performance. To enhance the value of the practical applications of the proposed method, we developed a standalone software for predicting oogenesis candidate proteins called OOgenesis_Pred. This software is the first predictor ever established for identifying oogenesis proteins. We also showed the capability of OOgenesis_Pred by making oogenesis-related proteins prediction for some of the oogenesis candidate proteins. It is anticipated that OOgenesis_Pred will become a powerful tool for future proteomic studies related to oogenesis.

In silico prediction of long intergenic non-coding RNAs in sheep.

Long non-coding RNAs (lncRNAs) are transcribed RNA molecules >200 nucleotides in length that do not encode proteins and serve as key regulators of diverse biological processes. Recently, thousands of long intergenic non-coding RNAs (lincRNAs), a type of lncRNAs, have been identified in mammalians using massive parallel large sequencing technologies. The availability of the genome sequence of sheep (Ovis aries) has allowed us genomic prediction of non-coding RNAs. This is the first study to identify lincRNAs using RNA-seq data of eight different tissues of sheep, including brain, heart, kidney, liver, lung, ovary, skin, and white adipose. A computational pipeline was employed to characterize 325 putative lincRNAs with high confidence from eight important tissues of sheep using different criteria such as GC content, exon number, gene length, co-expression analysis, stability, and tissue-specific scores. Sixty-four putative lincRNAs displayed tissues-specific expression. The highest number of tissues-specific lincRNAs was found in skin and brain. All novel lincRNAs that aligned to the human and mouse lincRNAs had conserved synteny. These closest protein-coding genes were enriched in 11 significant GO terms such as limb development, appendage development, striated muscle tissue development, and multicellular organismal development. The findings reported here have important implications for the study of sheep genome.

Neural network and SVM classifiers accurately predict lipid binding proteins, irrespective of sequence homology.

Due to the central roles of lipid binding proteins (LBPs) in many biological processes, sequence based identification of LBPs is of great interest. The major challenge is that LBPs are diverse in sequence, structure, and function which results in low accuracy of sequence homology based methods. Therefore, there is a need for developing alternative functional prediction methods irrespective of sequence similarity. To identify LBPs from non-LBPs, the performances of support vector machine (SVM) and neural network were compared in this study. Comprehensive protein features and various techniques were employed to create datasets. Five-fold cross-validation (CV) and independent evaluation (IE) tests were used to assess the validity of the two methods. The results indicated that SVM outperforms neural network. SVM achieved 89.28% (CV) and 89.55% (IE) overall accuracy in identification of LBPs from non-LBPs and 92.06% (CV) and 92.90% (IE) (in average) for classification of different LBPs classes. Increasing the number and the range of extracted protein features as well as optimization of the SVM parameters significantly increased the efficiency of LBPs class prediction in comparison to the only previous report in this field. Altogether, the results showed that the SVM algorithm can be run on broad, computationally calculated protein features and offers a promising tool in detection of LBPs classes. The proposed approach has the potential to integrate and improve the common sequence alignment based methods.

Deciphering the role of alternative splicing as a potential regulator in fat-tail development of sheep: a comprehensive RNA-seq based study.

Although research on alternative splicing (AS) has been widely conducted in mammals, no study has investigated the splicing profiles of genes involved in fat-tail formation in sheep. Here, for the first time, a comprehensive study was designed to investigate the profile of AS events and their involvement in fat-tail development of sheep. In total, 45 RNA-Seq samples related to seven different studies, which have compared the fat-tailed vs thin-tailed sheep breeds, were analyzed. Two independent tools, rMATS and Whippet, along with a set of stringent filters were applied to identify differential AS (DAS) events between the breeds per each study. Only DAS events that were detected by both tools as well as in at least three datasets with the same ΔPSI trend (percent spliced in), were considered as the final high-confidence set of DAS genes. Final results revealed 130 DAS skipped exon events (69 negative and 61 positive ΔPSI) belonged to 124 genes. Functional enrichment analysis highlighted the importance of the genes in the underlying molecular mechanisms of fat metabolism. Moreover, protein-protein interaction network analysis revealed that DAS genes are significantly connected. Of DAS genes, five transcription factors were found that were enriched in the biological process associated with lipid metabolism like "Fat Cell Differentiation". Further investigations of the findings along with a comprehensive literature review provided a reliable list of candidate genes that may potentially contribute to fat-tail formation including HSD11B1, SIRT2, STRN3 and TCF7L2. Based on the results, it can be stated that the AS patterns may have evolved, during the evolution of sheep breeds, as another layer of regulation to contribute to biological complexity by reprogramming the gene regulatory networks. This study provided the theoretical basis of the molecular mechanisms behind the sheep fat-tail development in terms of AS.

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.

Transcriptome alterations of radish shoots exposed to cadmium can be interpreted in the context of leaf senescence.

Till now few transcriptome studies have described shoot responses of heavy metal (HM)-sensitive plants to excess Cd and still a unifying model of Cd action is lacking. Using RNA-seq technique, the transcriptome responses of radish (Raphanus sativus L.) leaves to Cd stress were investigated in plants raised hydroponically under control and 5.0 mg L-1 Cd. The element was mainly accumulated in roots and led to declined biomass and photosynthetic pigments, increased H2O2 and lipid peroxidation, and the accumulation of sugars, protein thiols, and phytochelatins. Out of 524 differentially expressed genes (DEGs), 244 and 280 upregulated and downregulated ones were assigned to 82 and 115 GO terms, respectively. The upregulated DEGs were involved in osmotic regulation, protein metabolism, chelators, and carbohydrate metabolisms, whereas downregulated DEGs were related to photosynthesis, response to oxidative stress, glucosinolate, and secondary metabolite biosynthesis. Our transcriptome data suggest that Cd triggers ROS production and photosynthesis decline associated with increased proteolysis through ubiquitin-proteasome system (UPS)- and chloroplast-proteases and in this way brings about re-mobilization of N and C stores into amino acids and sugars. Meanwhile, declined glucosinolate metabolism in favor of chelator synthesis and upregulation of dehydrins as inferred from transcriptome analysis confers shoots some tolerance to the HM-derived ionic/osmotic imbalances. Thus, the induction of leaf senescence might be a major long-term response of HM-sensitive plants to Cd toxicity.

A nutrigenomics approach to study the effects of ω-3 fatty acids in laying hens under physiological stress.

Supplement of ω-3 fatty acids can decrease the harmful effects of stress. However, the potential molecular mechanisms that are modulated by dietary ω-3 fatty acids in laying hens under stress remain unknown. Hence, RNA-sequencing (RNA-Seq) technology was used to gain new insights into different gene expression profiles and potential pathways involved in response to stress in the liver of 35-week-old Lohmann LSL-Lite laying hens supplemented with ω-3. Three groups including control (non-stress), stress, and stress_ω-3 fatty acids (three layers per each group) were applied. A total of 1,321 genes were detected as differentially expressed genes of which 701, 1,049, and 86 DEGs belonged to stress vs. control, stress_ω-3 vs. control, and stress vs. stress_ω-3 pairwise comparisons, respectively. Gene ontology and KEGG pathway analysis indicated that the DEGs were enriched in particular regulation of steroid and cholesterol biosynthetic process, fatty acid degradation, AMPK signaling pathway, fatty acid biosynthesis, and immune response. Our data represented a promising approach regarding the importance of ω-3 as anxiolytic and anti-stress. In this context, UNC13B and ADRA1B genes were downregulated in the stress_ω-3 group compared to the stress group, which are associated with decreased activity of glutamatergic stimulatory neurons and probably play important role in facilitating the response to stress. This study extends the current understanding of the liver transcriptome response to physiological stress, and provides new insights into the molecular responses to stress in laying hens fed a diet supplemented with ω-3 fatty acids.

Comprehensive RNA-Seq-based study and metabolite profiling to identify genes involved in podophyllotoxin biosynthesis in Linum album Kotschy ex Boiss. (Linaceae).

Linum album is a well-known rich source of anticancer compounds, i.e., podophyllotoxin (PTOX) and other lignans. These compounds play an important role in the plant's defensive system. The RNA-Seq data of flax (L. usitatissimum) were analyzed under various biotic and abiotic stresses to comprehend better the importance of lignans in plant defense responses. Then, the association between the lignan contents and some related gene expressions was experimented with HPLC and qRT-PCR, respectively. Transcriptomic profiling showed a specific expression pattern in different organs, and just the commonly regulated gene EP3 was detected with a significant increase under all stresses. The in silico analysis of the PTOX biosynthesis pathway identified a list of genes, including laccase (LAC11), lactoperoxidase (POD), 4-coumarate-CoA ligase (4CL), and secoisolariciresinol dehydrogenase (SDH). These genes increased significantly under individual stresses. The HPLC analysis showed that the measured lignan contents generally increased under stress. In contrast, a quantitative expression of the genes involved in this pathway using qRT-PCR showed a different pattern that seems to contribute to regulating PTOX content in response to stress. Identified modifications of critical genes related to PTOX biosynthesis in response to multiple stresses can provide a baseline for improving PTOX content in L. album.

Integrated co-expression analysis of regulatory elements (miRNA, lncRNA, and TFs) in bovine monocytes induced by Str. uberis.

Non-coding RNAs, including long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), together with transcription factors, are critical pre-, co-, and post-transcriptional regulators. In addition to their criteria as ideal biomarkers, they have great potential in disease prognosis, diagnosis, and treatment of complex diseases. Investigation of regulatory mechanisms in the context of bovine mastitis, as most common and economic disease in the dairy industry, to identify elements influencing the expression of candidate genes as key regulators of the mammary immune response is not yet fully understood. Transcriptome profiles (50 RNA-Seq and 50 miRNA-Seq samples) of bovine monocytes induced by Str. uberis were used for co-expression module detection and preservation analysis using the weighted gene co-expression network analysis (WGCNA) approach. Assigned mi-, lnc-, and m-modules used to construct the integrated regulatory networks and miRNA-lncRNA-mRNA regulatory sub-networks. Remarkably, we have identified 18 miRNAs, five lncRNAs, and seven TFs as key regulators of str. uberis-induced mastitis. Most of the genes introduced here, mainly involved in immune response, inflammation, and apoptosis, were new to mastitis. These findings may help to further elucidate the underlying mechanisms of bovine mastitis, and the discovered genes may serve as signatures for early diagnosis and treatment of the disease.

Fatty Acids in Subcutaneous Adipose Tissue of Pregnant Women with and without Polycystic Ovary Syndrome Are Associated with Genes Related to Steroidogenesis: A Case-Control Study.

BACKGROUND: The qualitative analysis of adipose tissue (AT) is an exciting area for research and clinical applications in several diseases and it is emerging along with the quantitative approach to research on overweight and obese people. While the importance of steroid metabolism in women with polycystic ovary syndrome (PCOS) has been reported, limited data exists on the effective roles of AT in pregnant women suffering from PCOS. The aim of this study was to determine association of fatty acid (FA) profiles with expression of 14 steroid genes in abdominal subcutaneous AT of PCOS vs. non-PCOS pregnant women. MATERIALS AND METHODS: In this case-control study, the AT samples of 36 non-PCOS pregnant women and 12 pregnant women with PCOS (3:1 ratio control: case) who underwent cesarean section were collected. Relationship of expressing gene targets and different features were performed using Pearson correlation analysis on the R 3.6.2 software. The ggplot2 package in R tool was used to draw the plots. RESULTS: Age (31.4 and 31.5 years, P=0.99), body mass index (BMI) (prior pregnancy 26 and 26.5 kg.m-2, P=0.62) and at delivery day (30.1 and 31, P=0.94), gestational period (264 and 267 days, P=0.70) and parity (1.4 and 1.4, P=0.42) of non-PCOS and PCOS pregnant women were similar. Expression of steroidogenic acute regulator (STAR) and 11ß-Hydroxysteroid dehydrogenase (11BHSD2) in non-PCOS pregnant women showed the highest association with eicosapentaenoic acid (EPA, C20:5 n-3, r=0.59, P=0.001) and (r=0.66, P=0.001), respectively. In the all participants, STAR mRNA level showed the greatest association with the EPA fatty acid concentration (P=0.001, r=0.51). CONCLUSION: Our results showed a link between the genes involved in steroid metabolism and fatty acids in AT of pregnant women, especially for omega-3 FA and the gene involved in the first step of steroidogenesis in subcutaneous AT. These findings warrant further studies.

Artificial neural network (ANN)-based algorithms for high light stress phenotyping of tomato genotypes using chlorophyll fluorescence features.

High light (HL) is a common environmental stress directly imposes photoinhibition on the photosynthesis apparatus. Breeding plants for tolerance against HL is therefore highly demanded. Chlorophyll fluorescence (ChlF) is a sensitive indicator of stress in plants and can be evaluated using OJIP transients. In this study, we compared the ChlF features of plants exposed to HL (1200 µmol m-2 s-1) with that of control plants (300 µmol m-2 s-1). To extract the most reliable ChlF features for discrimination between HL-stressed and non-stressed plants, we applied three artificial neural network (ANN)-based algorithms, namely, Boruta, Support Vector Machine (SVM), and Recursive Feature Elimination (RFE). Feature selection algorithms identified multiple features but only two features, namely the maximal quantum yield of PSII photochemistry (FV/FM) and quantum yield of energy dissipation (ɸD0), remained consistent across all genotypes in control conditions, while exhibited variation in HL. Therefore, considered reliable features for HL stress screening. The selected features were then used for screening 14 tomato genotypes for HL. Genotypes were categorized into three groups, tolerant, semi-tolerant, and sensitive genotypes. Foliar hydrogen peroxide (H2O2) and malondialdehyde (MDA) contents were measured as independent proxies for benchmarking selected features. Tolerant genotypes were attributed with the lowest change in H2O2 and MDA contents, while the sensitive genotypes displayed the highest magnitude of increase in H2O2 and MDA by HL treatment compared to the control. Finally, a FV/FM higher than 0.77 and ɸD0 lower than 0.24 indicates a healthy electron transfer chain (ETC) when tomato plants are exposed to HL.

Alpha-linolenic acid alleviates the detrimental effects of lipopolysaccharide during in vitro ovine oocyte development.

During the transition period and early lactation of ruminants with higher production, the reproductive organs are exposed to various stressors, like inflammation stimulators such as lipopolysaccharides (LPS), as a consequence of high concentrate consumption. In this study, we aimed to determine the probable potential of α-linolenic acid (ALA) in alleviating LPS-induced effects in ovine oocytes in vitro as well as the underlying controlling mechanisms. Different concentrations of LPS (0, 0.01, 0.1, 1, and 10 µg/mL) were added to the oocyte maturation medium to evaluate its effect on oocyte developmental competence. Likewise, different concentrations of ALA (0, 10, 50, 100, and 200 µM/mL) were added to the maturation medium to define its effects on oocyte developmental competence. Accordingly, a combination of ALA and LPS in a dose-dependent manner was added to the maturation medium to elucidate their effect on oocyte developmental competence and uncover any possible potential of ALA to alleviate the detrimental effect induced by the presence of LPS. The expressions of candidate genes were measured in mature oocytes treated either with ALA, LPS, or ALA plus LPS. Adding LPS to the maturation medium decreased the cleavage rate of the treated oocytes, and those oocytes reached the blastocyst stage at a lower rate. Adding ALA to the maturation medium in the presence of LPS alleviated the detrimental effects of LPS in a dose-dependent manner, which ultimately led to higher cleavage and blastocyst formation. A higher expression of Trim26, GRHPR, NDUFA, PGC-1α, SOD, CS, SDH, p53, and CAT was observed in LPS-treated oocytes compared with the ALA and control groups. Additionally, CS and CAT transcripts were down-regulated in oocytes in LPS plus ALA-treated group compared to that of the LPS-treated group. These findings revealed that ALA has the potential to alleviate the detrimental effects induced by LPS on in ovine oocytes during maturation in vitro. Thus, LPS-detrimental effect and ALA-preventing mechanisms seem to be regulated through the expression of genes involved in mitochondrial biogenesis and function, oxidative stress, and antioxidant systems.

A framework for non-preserved consensus gene module detection in Johne's disease.

Johne's disease caused by Mycobacterium avium subsp. paratuberculosis (MAP) is a major concern in dairy industry. Since, the pathogenesis of the disease is not clearly known, it is necessary to develop an approach to discover molecular mechanisms behind this disease with high confidence. Biological studies often suffer from issues with reproducibility. Lack of a method to find stable modules in co-expression networks from different datasets related to Johne's disease motivated us to present a computational pipeline to identify non-preserved consensus modules. Two RNA-Seq datasets related to MAP infection were analyzed, and consensus modules were detected and were subjected to the preservation analysis. The non-preserved consensus modules in both datasets were determined as they are modules whose connectivity and density are affected by the disease. Long non-coding RNAs (lncRNAs) and TF genes in the non-preserved consensus modules were identified to construct integrated networks of lncRNA-mRNA-TF. These networks were confirmed by protein-protein interactions (PPIs) networks. Also, the overlapped hub genes between two datasets were considered hub genes of the consensus modules. Out of 66 consensus modules, 21 modules were non-preserved consensus modules, which were common in both datasets and 619 hub genes were members of these modules. Moreover, 34 lncRNA and 152 TF genes were identified in 12 and 19 non-preserved consensus modules, respectively. The predicted PPIs in 17 non-preserved consensus modules were significant, and 283 hub genes were commonly identified in both co-expression and PPIs networks. Functional enrichment analysis revealed that eight out of 21 modules were significantly enriched for biological processes associated with Johne's disease including "inflammatory response," "interleukin-1-mediated signaling pathway", "type I interferon signaling pathway," "cytokine-mediated signaling pathway," "regulation of interferon-beta production," and "response to interferon-gamma." Moreover, some genes (hub mRNA, TF, and lncRNA) were introduced as potential candidates for Johne's disease pathogenesis such as TLR2, NFKB1, IRF1, ATF3, TREM1, CDH26, HMGB1, STAT1, ISG15, CASP3. This study expanded our knowledge of molecular mechanisms involved in Johne's disease, and the presented pipeline enabled us to achieve more valid results.

RNA-Seq Analysis of Magnaporthe grisea Transcriptome Reveals the High Potential of ZnO Nanoparticles as a Nanofungicide.

Magnaporthe grisea is one of the most destructive pathogen that encounters a challenge to rice production around the worldwide. The unique properties of ZnO nanoparticles (NPs), have high attractiveness as nanofungicide. In the present study, the response of fungi to ZnO NPs was evaluated using RNA sequencing (RNA-seq). Two different aligners (STAR and Hisat2) were used for aligning the clean reads, and the DEseq2 package was used to identify the differentially expressed genes (DEGs). In total, 1,438 and 761 fungal genes were significantly up- and down-regulated in response to ZnO NPs, respectively. The DEGs were subjected to functional enrichment analysis to identify significantly enriched biological pathways. Functional enrichment analysis revealed that "cell membrane components," "ion (calcium) transmembrane transporter activity," "steroid biosynthesis pathway" and "catalytic activity" were the contributed terms to fungal response mechanisms. The genes involved in aflatoxin efflux pumps and ribosome maturation were among the genes showing significant up- and down-regulation after ZnO NPs application. To confirm the obtained RNA-seq results, the expression of six randomly selected genes were evaluated using q-RT-PCR. Overall, the RNA-seq results suggest that ZnO NPs primarily act on the fungal cell membrane, but accumulation of ROS inside the cell induces oxidative stress, the fungal catalytic system is disrupted, resulting into the inhibition of ROS scavenging and eventually, to the death of fungal cells. Our findings provide novel insights into the effect of ZnO NPs as a promising nanofungicide for effective control of rice blast disease.

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.