mRNA-Seq of testis and liver tissues reveals a testis-specific gene and alternative splicing associated with hybrid male sterility in dzo.

Alternative splicing (AS) plays an important role in the co-transcription and post-transcriptional regulation of gene expression during mammalian spermatogenesis. The dzo is the male F1 offspring of an interspecific hybrid between a domestic bull (Bos taurus ♂) and a yak (Bos grunniens ♀) which exhibits male sterility. This study aimed to identify the testis-specific genes and AS associated with hybrid male sterility in dzo. The iDEP90 program and rMATS software were used to identify the differentially expressed genes (DEG) and differential alternative splicing genes (DSG) based on RNA-seq data from the liver (n = 9) and testis (n = 6) tissues of domestic cattle, yak, and dzo. Splicing factors (SF) were obtained from the AmiGO2 and the NCBI databases, and Pearson correlation analysis was performed on the differentially expressed SFs and DSGs. We focused on the testis-specific DEGs and DSGs between dzo and cattle and yak. Among the top 3,000 genes with the most significant variations between these 15 samples, a large number of genes showed testis-specific expression involved with spermatogenesis. Cluster analysis showed that the expression levels of these testis-specific genes were dysregulated during mitosis with a burst downregulation during the pachynema spermatocyte stage. The occurrence of AS events in the testis was about 2.5 fold greater than in the liver, with exon skipping being the major AS event (81.89% to 82.73%). A total of 74 DSGs were specifically expressed in the testis and were significantly enriched during meiosis I, synapsis, and in the piRNA biosynthesis pathways. Notably, STAG3 and DDX4 were of the exon skipping type, and DMC1 was a mutually exclusive exon. A total of 36 SFs were significantly different in dzo testis, compared with cattle and yak. DDX4, SUGP1, and EFTUD2 were potential SFs leading to abnormal AS of testis-specific genes in dzo. These results show that AS of testis-specific genes can affect synapsis and the piRNA biosynthetic processes in dzo, which may be important factors associated with hybrid male sterility in dzo.

The interspecific hybrid offspring of a domestic bull (Bos taurus) and a yak (Bos grunniens) display heterosis in meat and milk production. The hybrid offspring are particularly adaptable to the harsh environments of the Qinghai-Tibet Plateau. However, the male F1 to F3 offspring of this interspecies hybrid are infertile, and spermatogenesis is arrested at meiosis preventing the prolonged utilization of the benefits of heterosis. This study aimed to identify the testis-specific genes and alternative splicing (AS) associated with hybrid male sterility using RNA-Seq data from the liver and testis tissues of domestic cattle, yaks, and their F1 offspring (dzo). The expression of the testis-specific genes became disordered during mitosis and meiosis in dzo. Their testis-specific genes with AS events were enriched during synapsis and in the piRNA biosynthetic processes. In addition, we identified the potential splicing factors associated with abnormal testis-specific AS gene expression in dzo. These results reveal the important role of AS in the meiotic arrest of dzo.

Identification and screening of circular RNAs during adipogenic differentiation of ovine preadipocyte by RNA-seq.

Circular RNAs (circRNAs) are a class of non-coding RNAs that play important roles in preadipocyte differentiation and adipogenesis. However, little is known about genome-wide identification, expression profile, and function of circRNAs in sheep. To investigate the role of circRNAs during ovine adipogenic differentiation, the subcutaneous adipose tissue of Tibetan rams was collected in June 2022. Subsequently, the preadipocytes were immediately isolated from collected adipose tissue and then induced to begin differentiation. The adipocytes samples cultured on days 0, 2, and 8 of preadipocytes differentiation were used to perform RNA sequencing (RNA-seq) analysis to construct the expression profiles of circRNAs. Subsequently, the function of differentially expressed circRNAs was investigated by performing the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of their parent genes. Finally, a circRNAs-miRNAs-mRNAs network involved in adipogenic differentiation was been analyzed. As a result, a total of 6,449 candidate circRNAs were identified in ovine preadipocytes. Of these circRNAs identified, 63 candidate circRNAs were differentially expressed among the three differentiation stages and their parent genes were mainly enriched in acetyl-CoA metabolic process, positive regulation of lipid biosynthetic process, positive regulation of steroid biosynthetic process, and focal adhesion pathway (P < 0.05). Based on a circRNAs-miRNAs-mRNAs regulatory network constructed, circ_004977, circ_006132 and circ_003788 were found to function as competing endogenous RNAs (ceRNAs) to regulate ovine preadipocyte differentiation and lipid metabolism. The results provide an improved understanding of functions and molecular mechanisms of circRNAs underlying ovine adipogenesis in sheep.

The moderate fat deposition contributes to improve mutton quality, which is associated with the differentiation of preadipocytes. To investigate roles of circular RNAs (circRNAs) in preadipocyte differentiation, we identified circRNAs on days 0, 2, and 8 of preadipocytes differentiation and compared the expression profile of circRNAs at different adipogenic differentiation stages. A total of 6,449 candidate circRNAs were identified, among which 63 candidate circRNAs were differentially expressed among the three differentiation stages. The parent genes of differentially expressed circRNAs were enriched in several biological process and pathways related to lipid metabolism and synthesis. In addition, several circRNAs may regulate ovine preadipocyte differentiation by interacting with microRNAs (miRNAs). The results reveal the potential roles of circRNAs in adipogenic differentiation of sheep.

Combined ATAC-seq, RNA-seq, and GWAS analysis reveals glycogen metabolism regulatory network in Jinjiang oyster ( Crassostrea ariakensis).

Glycogen serves as the principal energy reserve for metabolic processes in aquatic shellfish and substantially contributes to the flavor and quality of oysters. The Jinjiang oyster ( Crassostrea ariakensis) is an economically and ecologically important species in China. In the present study, RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) were performed to investigate gene expression and chromatin accessibility variations in oysters with different glycogen contents. Analysis identified 9 483 differentially expressed genes (DEGs) and 7 215 genes with significantly differential chromatin accessibility (DCAGs) were obtained, with an overlap of 2 600 genes between them. Notably, a significant proportion of these genes were enriched in pathways related to glycogen metabolism, including "Glycogen metabolic process" and "Starch and sucrose metabolism". In addition, genome-wide association study (GWAS) identified 526 single nucleotide polymorphism (SNP) loci associated with glycogen content. These loci corresponded to 241 genes, 63 of which were categorized as both DEGs and DCAGs. This study enriches basic research data and provides insights into the molecular mechanisms underlying the regulation of glycogen metabolism in C. ariakensis.

Characterization and comparative transcriptomic analysis of skeletal muscle in female Pekin duck and Hanzhong Ma duck during different growth stages using RNA-seq.

Duck is an economically important poultry, and there is currently a major focus on improving its meat quality through breeding. There are wide variations in the growth regulation mechanisms of different duck breeds, that fundamental research on skeletal muscle growth is essential for understanding the regulation of unknown genes. The study aimed to broaden the understanding the duck skeletal muscle development and thereby to improve the performance of domestic ducks. In this study, RNA-seq data from skeletal muscles (breast muscle and leg muscle) of Pekin duck and Hanzhong Ma duck sampled at d 17, 21, and 27 of embryo (E17d, E21d, and E27d), as well as at 6-mo-old following birth (M6), to investigate and compare the mRNA temporal expression profiles and associated pathways that regulate skeletal myogenesis of different duck breeds. There were 331 to 1,440 annotated differentially expressed genes (DEGs) in breast muscle and 380 to 1,790 annotated DEGs in leg muscle from different databases between 2 duck breeds. Gene ontology (GO) enrichment in skeletal muscles indicated that these DEGs were mainly involved in biosynthetic process, developmental process, regulation of protein metabolic process and regulation of gene expression. KEGG analysis in skeletal muscles showed that a total of 41 DEGs were mapped to 7 KEGG pathways, including ECM-receptor interaction, focal adhesion, carbon metabolism, regulation of actin cytoskeleton, calcium signaling pathway, biosynthesis of amino acids and PPAR signaling pathway. The differential expression of 8 selected DEGs was verified by qRT-PCR, and the results were consistent with RNA-seq data. The identified DEGs, such as SDC, SPP1, PAK1, MYL9, PGK1, NOS1, PHGDH, TNNT2, FN1, and AQP4, were specially highlighted, indicating their associations with muscle development in the Pekin duck and Hanzhong Ma duck. This study provides a basis for revealing the differences in skeletal muscle development between Pekin duck and Hanzhong Ma duck.

RNA-seq reveals the role of miR-223 in alleviating inflammation of bovine mammary epithelial cells.

Bovine mammary epithelial cells (bMECs) are involved in the early defense against the invasion of intramammary pathogens and are essential for the health of bovine mammary gland. MicroRNA (MiRNA) is a key factor that regulates cell state and physiological function. In the present study, the transcriptome profiles of miR-223 inhibitor transfection group (miR-223_Inhibitor) and negative control inhibitor transfection group (NC_Inhibitor) within bMECs were detected via the RNA sequencing (RNA-seq) platform. Based on these experiments, the differentially expressed mRNAs (DE-mRNAs) of the miR-223_Inhibitor transfection group were screened, and the Gene Ontology and Kyoto Encyclopedia of Genes and Genomes functional analyses of DE-mRNAs were performed. The results revealed that compared with the NC_Inhibitor, 224 differentially expressed genes (DEGs) were identified in the miR-223_Inhibitor, including 184 upregulated and 40 downregulated genes. The functional annotation of the above DEGs indicated that some of these genes are involved in the immune response generated by extracellular substance stimulation, regulation of the activity of cytokines and chemokines, and the immune signaling pathways of NF-κB and TNF. Meanwhile, miR-223_inhibitor upregulated the immune key genes IRF1 and NFκBIA, cytokines IL-6 and IL-24, as well as chemokines CXCL3, CXCL5, and CCR6, triggering a signaling cascade response that exacerbated inflammation in bMECs. These results suggested that miR-223 plays an important role in inhibiting the inflammatory response and maintaining the stability of bMECs, and is a potential target for treating mastitis in dairy cows.

Research Note: Combined analysis of BSA-seq based mapping and RNA-seq reveals candidate genes associated with sub-Columbian plumage in H line chickens.

Columbian coloration patterns in plumage are widespread phenomena in several standard breeds of poultry, such as the Columbian Plymouth Rock chicken. H line chicken plumage is generally a pure white except in the hackle, wing, and tail plumage, which coloration is very similar to the Columbian plumage pattern, but with the barring substituting for the black vertical striping. Thus, we refer to this plumage coloration as "sub-Columbian" pattern. However, the genetic basis of this phenotype remains unknown. Here, a F3 cross population between yellow plumage roosters and sub-Columbian plumage hens was constructed, for verifying sub-Columbian plumage was sex-linked dominant inheritance. To identify the candidate regions, F2 generation sub-Columbian plumage hens and yellow plumage hens with their parental lines were used for BSA-seq, and sub-Columbian plumage genes were mapped to a 10.46 Mb interval on chromosome Z. Remarkably, by transcriptome analysis of the neck and wing tip follicle tissues of the 2 plumage colors, we demonstrated that within the interval, only 1 gene, SLC45A2 expressed significant differently (P < 0.05). Through KASP, we identified L347M and A10331272T in solute carrier family 45 member 2 (SLC45A2), and B2 haplotype of cyclin-dependent kinase inhibitor 2A (CDKN2A), showed near complete association with the phenotype. Eventually, we designed a hybridization experiment for verifying the locus of sub-Columbian plumage, which is inherited through Z-linked dominant inheritance and is controlled by SLC45A2 and CDKN2A.

Comprehensive RNA-seq profiling to evaluate the rabbit mammary gland transcriptome after mastitis.

Mastitis is a relatively common disease in rabbit does. The aim of this study was to investigate a relationship between the severity of clinical signs and pathological observations and to analyze differentially expressed genes (DEGs) in the mammary gland with mastitis versus healthy mammary gland. The result showed that rectal temperatures of the rabbits with both mild mastitis and severe mastitis were higher than that of control. Cell counting results showed that the somatic cell count (SCC) only in milk of the rabbit with severe mastitis was significantly higher than that in the control group. However, the number of heterophils in the histological sections of mammary glands with mild mastitis was significantly higher than that of control. A total of 1,096 DEGs between the control and mastitis mammary glands was identified by RNA-sequencing (RNA-seq). Gene ontology (GO) showed that most of up-regulated genes were enriched in terms such as response to stimulus, signal transduction, and cell communication. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that these genes were mostly enriched in the pathways such as Rap1 signaling pathway, proteoglycans in cancer, and PI3K-Akt signaling pathway. However, the downregulated genes were mainly enriched in metabolic processes and significantly involved in metabolic pathways. The data provides useful information to further dissect the molecular genetic mechanisms underlying rabbit mastitis, which is a prerequisite for designing effective intervention strategies.

Gene expression pattern and molecular mechanisms involved in Shal and Sangsari sheep fertility using RNA-Seq.

Ovulation rate and litter size are the main reproductive traits with high economic value in the sheep breeding industry. In this study, three Shal ewes (multiparous) and three Sangsari ewes (uniparous) at the age of 5 were used. The live weight was between 45 and 50 kg at an extremely body condition score of 3. These breeds are marked seasonal reproduction activity and are often bred in semi-closed breeding systems. Total RNAs were extracted from the ovarian tissues, and RNA sequencing was carried out. The DAVID (Database for Annotation, Visualization and Integrated Discovery) database was then used to annotate genes, and the string database and the Cytoscape software were used to investigate their interactions. Then path-act network analysis and gene-act network analysis were investigated. The results indicated that 19,932 genes were differentially expressed. The 5968 differentially expressed genes were identified in Shal ewe's ovarian tissue compared to Sangsari ewes (FDR < 0.05), of which 2921 genes were up-regulated and 3047 genes were down-regulated. Bioinformatics analysis exhibited that most of the biological processes and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways associated with significant DEGs (Differentially Expressed Genes) in the two studied breeds are associated with oocyte maturation and metabolism. MAPK signalling pathways and Ubiquitin-mediated proteolysis are the most important biological pathways associated with reproductive and fertility traits in the Shal breed. AKT3, MAPK8, MAPK9 and RELA genes are also important genes related to the fertility of multiparous sheep. Analysis of ovarian RNA-seq data identified that most of the differentially expressed genes were involved in various reproductive processes including folliculogenesis, ovulation, ovarian and embryonic development. The MAPK signalling pathway had the most interaction with other pathways, and the AKT3 gene could be a powerful candidate gene in the reproduction and fertility of Shal sheep. These results could pave the way for future efforts to address sheep prolificacy barriers.

Exploring differentially expressed genes in hypothalamic, pars tuberalis and pineal gland transcriptomes in different sexual behavior phenotypes in rams using RNA-Seq.

Reproductive seasonality is a limiting factor in sheep production. Sexual behavior is a key element in reproductive efficiency, and this function is regulated by the hypothalamus-pituitary-gonadal (HPG) axis. To understand the mechanisms of sexual behavior, transcriptomic sequencing technology was used to identify differentially expressed genes (DEGs) in the hypothalamus (HT), pars tuberalis (PT) and pineal gland (PG) in Rasa Aragonesa rams with different sexual behavior. Bioinformatics analysis of the 16,401 identified genes by RNA-Seq revealed 103 and 12 DEGs in the HT and the PG, respectively, at a false discovery rate (FDR) of 5% with an absolute value of expression ≥ 1 (log2FC). However, no DEGs were found in the PT. Functional annotation and pathway enrichment analysis showed that DEGs of HT were enriched mainly in neuroactive ligand-receptor interactions and signaling pathways, including notable candidate genes such as MTNR1A, CHRNA2, FSHB, LHB, GNRHR, AVP, PRL, PDYN, CGA, GABRD, and TSHB, which play a crucial role in sexual behavior. The GnRH and cAMP signaling pathways were also highlighted. In addition, gene set enrichment analysis (GSEA) identified potential pathways, dominated mainly by biological process category, that could be responsible for the differences in sexual behavior observed in rams. The intracellular protein transport and pattern specification process were enriched within the PT and the transcription factor binding and protein ubiquitination pathways for the PG. Thus, these pathways together may play an important role in the regulation of the sexual behavior in Rasa Aragonesa rams through the hypothalamic-pituitary-gonadal axis. The validation of 5 DEGs using reverse transcription quantitative polymerase chain reaction (RT-qPCR) showed expression patterns like the found with RNA-Seq. Overall, these results contribute to understanding the genomic basis of sexual behavior in rams. Our study demonstrates that multiple networks and pathways orchestrate sexual behavior in sheep.

Male sexual behavior is a key factor in reproduction, especially in seasonal breeders such as sheep. The identification of differentially expressed genes (DEGs) in brain regions involved in male reproduction and sexual behavior between rams with different sexual activity by RNA high-throughput sequencing can provide useful information to the sheep meat industry. This work aimed to determine the possible molecular mechanisms underlying the sexual behavior of Rasa Aragonesa rams by investigating transcriptional changes in the hypothalamus (HT), pars tuberalis (PT) and pineal gland (PG) between active (A) and nonactive (NA) rams. Comparative analysis revealed 103 and 12 DEGs between the A vs. NA comparison in the HT and the PG, respectively, but no DEGs were found in the PT. Gene ontology (GO) enrichment analysis of DEGs in HT samples revealed significant pathways, associated mainly with neuroactive ligand-receptor interactions, and the GnRH and cAMP signaling pathways. Furthermore, gene set enrichment analysis (GSEA) detected many overrepresented pathways related to sexual behavior via an interaction network within the hypothalamic-pituitary-gonadal axis. These data will be helpful for further investigations to look for mutations or functional single nucleotide polymorphisms (SNPs) that may be used for genetic assisted selection to improve sexual behavior in sheep.

Exploration the role of INHBA in Hu sheep granulosa cells using RNA-Seq.

Activin/inhibin is an important factor for the fecundity of Hu sheep, and it is involved in follicular development in ovaries. Inhibin subunit beta A (INHBA) participates in the synthesis of activin A and inhibin A. In this study, we also noted a positive correlation between INHBA level and the secretion of both activin A and inhibin A in culture medium. Nevertheless, both knockdown and overexpression of INHBA downregulated the expression of Inhibin Subunit Alpha (INHA). Based on RNA-Sequencing, we further examined the effect and molecular mechanism of INHBA knockdown in GCs on mRNA expression. A total of 1,687 differentially expressed genes (DEGs) were identified (Fold change ≥ 2; False-discovory-rates (FDR) ≤ 0.01), of which 602 genes were upregulated and 1,087 genes were downregulated in the INHBA interference group compared with the control groups. Gene Ontology (GO) enrichment indicated that these DEGs were mainly involved in the regulation of cell cycle, protein serine/threonine kinase activity, and actin cytoskeleton reorganization. Moreover, DEGs were significantly enriched in 40 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including P53, progesterone-mediated oocyte maturation, and PI3K-AKT signaling pathways. We also noted a positive correlation between INHBA level and many PI3K/Akt/mTOR pathway-related genes at the gene or/and protein expression. Overall, this study may contribute to a better understanding of the roles of INHBA on GCs of prolific sheep, as well as the molecular effect of low INHBA expression on GCs, clarifying some reproductive failures.

RNA-seq identifies differentially expressed genes involved in csal1 overexpression in granulosa cells of prehierarchical follicles in Chinese Dagu hens.

The transcription factor csal1 is an important molecule that plays a critical regulatory function in ovarian follicle development, as confirmed by our previous data. However, the candidate genes of csal1 and its regulatory mechanism remain poorly understood in the granulosa cells (GCs) of chicken prehierarchical follicles (PFs). Six transcriptomes of csal1 and empty vector were analyzed in Chinese Dagu hens by RNA sequencing. Six cDNA libraries were constructed, with more than 42 million clean reads and 16,779 unigenes. Of these 16,779 unigenes, 2,762 differentially expressed genes (DEGs) were found in GCs, including 1,605 upregulated and 1,157 downregulated unigenes. Fourteen genes, including BMP5, TACR2, AMH, PLAG1, MYOD1, BOP1, SIPA1, NOTCH1, BCL2L1, SOX9, ADGRA2, WNT5A, SLC7A11, and GATAD2B, were related to GC proliferation and differentiation, hormone production, ovarian follicular development, regulation of reproductive processes, and signaling pathways in the PFs. Further analysis demonstrated the DEGs in GCs of ovarian follicles were enriched in neuroactive ligand-receptor interaction, cell adhesion molecules, and pathways related to cytochrome P450, indicating a critical function for csal1 in the generation of egg-laying features by controlling ovarian follicle development. For the first time, the current study represents the transcriptome analysis with ectopic csal1 expression. These findings provide significant evidence for investigating the molecular mechanism by which csal1 controls PF development in the hen ovary.

Exploring differentially expressed genes related to metabolism by RNA-Seq in porcine embryonic fibroblast after insulin treatment.

BACKGROUND: Insulin regulates glucose homeostasis and has important effects on metabolism, cell growth, and differentiation. Depending on the cell type and physiological context, insulin signal has specific pathways and biological outcomes in different tissues and cells. For studying the signal pathway of insulin on glycolipid metabolism in porcine embryonic fibroblast (PEF), we used high-throughput sequencing to monitor gene expression patterns regulated by insulin. OBJECTIVES: The goal of our research was to see how insulin affected glucose and lipid metabolism in PEFs. METHODS: We cultured the PEFs with the addition of insulin and sampled them at 0, 48, and 72 h for RNA-Seq analysis in triplicate for each time point. RESULTS: At 48 and 72 h, 801 and 1,176 genes were differentially expressed, respectively. Of these, 272 up-regulated genes and 264 down-regulated genes were common to both time points. Gene Ontology analysis was used to annotate the functions of the differentially expressed genes (DEGs), the biological processes related to lipid metabolism and cell cycle were dominant. And the DEGs were significantly enriched in interleukin-17 signaling pathway, phosphatidylinositol-3-kinase-protein kinase B signaling pathway, pyruvate metabolism, and others pathways related to lipid metabolism by Kyoto Encyclopedia of Genes and Genomes enrichment analysis. CONCLUSIONS: These results elucidate the transcriptomic response to insulin in PEF. The genes and pathways involved in the transcriptome mechanisms provide useful information for further research into the complicated molecular processes of insulin in PEF.

RNA-seq and LC-MS/MS analysis of antiviral effects mediated by cold stress and stress hormone corticosterone in chicken DF-1 cells.

Cold stress in poultry is a global problem that causes huge economic losses and threatens the health and welfare of poultry. However, knowledge of chicken responses to virus infection under cold stress is limited. The purpose of this research was to investigate the effects of cold stress on gene expression and viral replication in chicken DF-1 cells in hypothermia. In addition, the characterization of circulating steroid hormone profiles in the plasma of chickens under cold stress was analyzed by liquid chromatography-tandem mass spectrometry. Herein, we performed RNA sequencing to obtain DF-1 cell transcriptional profiles under cold stress. A total of 9499 differentially expressed genes (DEGs) were identified in DF-1 cells. Overexpressed DEGs were related to the proteasome, cell cycle, spliceosome, ribosome biogenesis, and mammalian target of rapamycin (mTOR). Down-regulated DEGs were related to ribosomes, oxidative phosphorylation, apoptosis, and the p53 signaling pathway. Gene set enrichment analysis showed that the DEGs mainly affect host ribosome translation and mitochondrial respiratory electron transport. The principal steroid hormone alterations in chickens subjected to cold stress included dihydrotestosterone, testosterone, ß-sitosterol, androstenedione, 7a,27-dihydroxycholesterol,7-ketocholesterol, and desmosterol, which are associated with endocrine resistance, ovarian steroidogenesis, and steroid hormone biosynthesis. In addition, Infectious bronchitis virus (IBV), Newcastle disease virus (NDV), and Influenza A (H9N2) Virus replication in DF-1 cells is significantly inhibited by cold stress. Moreover, the plasma concentrations of corticosterone, an important stress hormone in poultry, were significantly elevated in chickens subjected to cold stress, and we found that IBV and vesicular stomatitis virus (VSV) replication were strongly inhibited in DF-1 cells pretreated with CORT, but NDV and H9N2 replication were unaffected. In conclusion, in response to cold stress, the translation efficiency and mitochondrial respiratory chain are temporarily weakened in DF-1 cells, which affects virus replication. Chickens may regulate aromatase deficiency, androstenedione metabolism, androgen and estrogen metabolism, and 17-beta hydroxysteroid dehydrogenase III deficiency through steroid hormones in response to cold stress. This study provides valuable insights into the molecular regulatory mechanisms of poultry under cold stress and may support further research on the intrinsic link between steroid hormones and virus replication under stress.

Delineating transcriptional crosstalk between Mycobacterium avium subsp. paratuberculosis and human THP-1 cells at the early stage of infection via dual RNA-seq analysis.

Mycobacterium avium subsp. paratuberculosis (MAP) is the causative agent of Johne's disease, a chronic debilitating disease in ruminants. To control this disease, it is crucial to understand immune evasion and the mechanism of persistence by analyzing the early phase interplays of the intracellular pathogens and their hosts. In the present study, host-pathogen interactions at the transcriptomic level were investigated in an in vitro macrophage infection model. When differentiated human THP-1 cells were infected with MAP, the expression of various genes associated with stress responses and metabolism was altered in both host and MAP at 3 h post-infection. MAP upregulates stress-responsive global gene regulators, such as two-component systems and sigma factors, in response to oxidative and cell wall stress. Downstream genes involved in type VII secretion systems, cell wall synthesis (polyketide biosynthesis proteins), and iron uptake were changed in response to the intracellular environment of macrophages. On the host side, upregulation of inflammatory cytokine genes was observed along with pattern recognition receptor genes. Notably, alterations in gene sets involved in arginine metabolism were observed in both the host and MAP, along with significant downregulation of NOS2 expression. These observations suggest that the utilization of metabolites such as arginine by intracellular MAP might affect host NO production. Our dual RNA-seq data can provide novel insights by capturing the global transcriptome with higher resolution, especially in MAP, thus enabling a more systematic understanding of host-pathogen interactions.

Identification of differentially expressed long noncoding RNAs in the ovarian tissue of Shal and Sangsari ewes using RNA-seq.

BACKGROUND: The ovary has an important role in reproductive function. Animal reproduction is dominated by numerous coding genes and noncoding elements. Although long noncoding RNAs (LncRNAs) are important in biological activity, little is known about their role in the ovary and fertility. METHODS: Three adult Shal ewes and three adult Sangsari ewes were used in this investigation. LncRNAs in ovarian tissue from two breeds were identified using bioinformatics analyses, and then target genes of LncRNAs were discovered. Target genes were annotated using the DAVID database, and their interactions were examined using the STRING database and Cytoscape software. The expression levels of seven LncRNAs with their target genes were assessed by real-time PCR to confirm the RNA-seq. RESULTS: Among all the identified LncRNAs, 124 LncRNAs were detected with different expression levels between the two breeds (FDR < 0.05). According to the DAVID database, target genes were discovered to be engaged in one biological process, one cellular component, and 21 KEGG pathways (FDR < 0.05). The PES1, RPS9, EF-1, Plectin, SURF6, CYC1, PRKACA MAPK1, ITGB2 and BRD2 genes were some of the most crucial target genes (hub genes) in the ovary. CONCLUSION: These results could pave the way for future efforts to address sheep prolificacy barriers.

RNA-seq reveals insights into molecular mechanisms of metabolic restoration via tryptophan supplementation in low birth weight piglet model.

Low birth weight (LBW) is associated with metabolic disorders in early life. While dietary l-tryptophan (Trp) can ameliorate postprandial plasma triglycerides (TG) disposal in LBW piglets, the genetic and biological basis underlying Trp-caused alterations in lipid metabolism is poorly understood. In this study, we collected 24 liver samples from 1-mo-old LBW and normal birth weight (NBW) piglets supplemented with different concentrations of dietary Trp (NBW with 0% Trp, N0; LBW with 0% Trp, L0; LBW with 0.4% Trp, L4; LBW with 0.8% Trp, L8; N = 6 in each group.) and conducted systematic, transcriptome-wide analysis using RNA sequencing (RNA-seq). We identified 39 differentially expressed genes (DEG) between N0 and L0, and genes within "increased dose effect" clusters based on dose-series expression profile analysis, enriched in fatty acid response of gene ontology (GO) biological process (BP). We then identified RNA-binding proteins including SRSF1, DAZAP1, PUM2, PCBP3, IGF2BP2, and IGF2BP3 significantly (P < 0.05) enriched in alternative splicing events (ASE) in comparison with L0 as control. There were significant positive and negative relationships between candidate genes from co-expression networks (including PID1, ANKRD44, RUSC1, and CYP2J34) and postprandial plasma TG concentration. Further, we determined whether these candidate hub genes were also significantly associated with metabolic and cardiovascular traits in humans via human phenome-wide association study (Phe-WAS), and analysis of mammalian orthologs suggests a functional conservation between human and pig. Our work demonstrates that transcriptomic changes during dietary Trp supplementation in LBW piglets. We detected candidate genes and related BP that may play roles on lipid metabolism restoration. These findings will help to better understand the amino acid support in LBW metabolic complications.

Low birth weight (LBW) has been associated with higher rate of mortality and morbidity and the development of metabolic complications, leaving burdens on livestock production and human health care. The feasibility of LBW metabolic restoration via postnatal nutrition compensation has been verified and the role of one of essential amino acids, l-tryptophan (Trp), on rescuing lipid metabolism in LBW was determined, while the underlying molecular mechanism and key gene regulation is little known. Our study was conducted to identify the unique molecular mechanisms between LBW and normal birth weight (NBW), and to identify the metabolic restoration related genes and biological processes after dietary Trp supplementation in LBW piglet model. We found that differentially expressed genes (DEG) between LBW and NBW were related to fatty acid response based on gene ontology enrichment analysis, and LBW piglets supplemented with Trp showed lower postprandial plasma triglycerides (TG) level as NBW, with similar expression feature of lipid metabolism related genes.

RNA-Seq reveals the functional specificity of epididymal caput, corpus, and cauda genes of cattleyak.

The first filial generation of the cattleyaks demonstrates hybrid vigor; however, the male cattleyaks are infertile and restrict productivity and breeding. The discovery of genes in a segment-specific approach offers valuable information and understanding concerning fertility status, yet the biology of cattleyak epididymis is still progressing. Comparative transcriptome analysis was performed on segment pairs of cattleyak epididymis. The caput versus corpus epididymis provided the highest (57.8%) differentially expressed genes (DEGs), corpus versus cauda (25.1%) followed, whereas caput versus cauda pair (17.1%) had the least DEGs. The expression levels of genes coding EPHB6, TLR1, MUC20, MT3, INHBB, TRPV5, EI24, PAOX, KIF12, DEPDC5, and KRT25, which might have the potentials to regulate the homeostasis, innate immunity, differentiation, motility, transport, and sperm maturation-related function in epididymal cells, were downregulated in the distal segment of epididymis. Top enriched KEGG pathways included mTOR, axon guidance, and taste transduction signaling pathways. EIF4B, EPHB6, and TAS2R42 were enriched in the pathways, respectively. Identifying key, new, and unexplored DEGs among the epididymal segments and further analyzing them could boost cattleyak fertility by maximizing sperm quality from genetically better sires and also facilitate better understanding of the epididymal biology.

Integrated analysis of lncRNAs and mRNAs by RNA-Seq in secondary hair follicle development and cycling (anagen, catagen and telogen) of Jiangnan cashmere goat (Capra hircus).

BACKGROUND: Among the world's finest natural fiber composites is derived from the secondary hair follicles (SHFs) of cashmere goats yield one of the world's best natural fibres. Their development and cycling are characterized by photoperiodism with diverse, well-orchestrated stimulatory and inhibitory signals. Long non-coding RNA (lncRNAs) and mRNAs play important roles in hair follicle (HF) development. However, not many studies have explored their specific functions in cashmere development and cycling. This study detected mRNAs and lncRNAs with their candidate genes and related pathways in SHF development and cycling of cashmere goat. We utilized RNA sequencing (RNA-Seq) and bioinformatics analysis on lncRNA and mRNA expressions in goat hair follicles to discover candidate genes and metabolic pathways that could affect development and cycling (anagen, catagen, and telogen). RESULTS: We identified 228 differentially expressed (DE) mRNAs and 256 DE lncRNA. For mRNAs, catagen and anagen had 16 upregulated and 35 downregulated DEGs, catagen and telogen had 18 upregulated and 9 downregulated DEGs and telogen and anagen had 52 upregulated and 98 downregulated DEGs. LncRNA witnessed 22 upregulated and 39 downregulated DEGs for catagen and anagen, 36 upregulated and 29 downregulated DEGs for catagen and telogen as well as 66 upregulated and 97 downregulated DEGs for telogen and anagen. Several key genes, including MSTRG.5451.2, MSTRG.45465.3, MSTRG.11609.2, CHST1, SH3BP4, CDKN1A, GAREM1, GSK-3ß, DEFB103A KRTAP9-2, YAP1, S100A7A, FA2H, LOC102190037, LOC102179090, LOC102173866, KRT2, KRT39, FAM167A, FAT4 and EGFL6 were shown to be potentially important in hair follicle development and cycling. They were related to, WNT/ß-catenin, mTORC1, ERK/MAPK, Hedgehog, TGFß, NFkB/p38MAPK, caspase-1, and interleukin (IL)-1a signaling pathways. CONCLUSION: This work adds to existing understanding of the regulation of HF development and cycling in cashmere goats via lncRNAs and mRNAs. It also serves as theoretical foundation for future SHF research in cashmere goats.

RNA-seq analysis of murine peyer's patches at 6 and 18 h post infection with Fasciola hepatica metacecariae.

Fasciola hepatica is a zoonotic parasite that not only economically burdens the agribusiness sector, but also infects up to 1 million people worldwide, with no commercial vaccine yet available. An ideal vaccine would induce protection in the gut, curtailing the extensive tissue damage associated with parasite's migration from the gut to the bile ducts. The design of such a vaccine requires greater knowledge of gut mucosal responses during the early stage of infection. We examined total mRNA expression of the peyer's patches at 6 and 18 h post F. hepatica infection using RNA sequencing. Differential expression analysis revealed 1341 genes upregulated and 61 genes downregulated at 6 h post infection, while 1562 genes were upregulated and 10 genes downregulated after 18 h. Gene-set enrichment analysis demonstrated that immune specific biological processes were amongst the most downregulated. The Toll-like receptor pathway in particular was significantly affected, the suppression of which is a well-documented immune evasive strategy employed by F. hepatica. In general, the genes identified were associated with suppression of inflammatory responses, helminth induced immune responses and tissue repair/homeostasis. This study provides a rich catalogue of the genes expressed in the early stages of F. hepatica infection, adding to the understanding of early host-parasite interactions and assisting in the design of future studies that look to advance the development of a novel F. hepatica vaccine.

Identifying the heat resistant genes by multi-tissue transcriptome sequencing analysis in Turpan Black sheep.

Heat stress not only affects the physical condition but also affects reproductive performance in sheep. A thorough understanding of the molecular and physiological mechanisms underlying heat stress would certainly improve livestock productivity and provide genetic evaluation ways for heat resistant breeds selection. In this study, 85 Turpan Black sheep, a breed exhibited excellent heat resistance from long-term artificial selection, and 85 heat sensitive Kazakh sheep in Turpan basin were tested for physiological and reproductive performance from July to August in summer. The results showed that the estrus rate was significantly higher in Turpan Black sheep (P < 0.05), while the heart rate and respiratory rate of Turpan Black sheep are significantly lower than that of Kazakh sheep (P < 0.05). Furthermore, to clarify genes participated in heat stress response, the pituitary, ovarian and hepatic tissues from three Turpan Black sheep and three Kazakh sheep were subjected to RNA-seq. The results indicated that 32, 49 and 69 genes were up-regulated, and 39, 60 and 145 genes were down-regulated in pituitary, ovarian and hepatic tissues in Turpan Black sheep compared with that of the Kazakh sheep, respectively. KEGG and gene set enrichment analysis showed that the differentially expressed genes were mainly involved in signal transduction pathways. In particular, the differentially expressed genes in hepar were enriched in the energy metabolism pathway, while the differentially expressed genes in ovarian tissue were enriched in the ovarium steroidogenesis pathway. In conclusion, our results implied that the pituitary-ovary axis might include hepar as downstream targeted organism in heat resistant regulation. Under heat stress, the signals released from pituitary would impact steroidogenesis in ovary, and further alter energy metabolism in hepar. As we know, this is the first comparative study to investigate the gene expression in multi-tissue in sheep under heat stress.