scholarly journals Changes in Transcriptomic Profiles in Different Reproductive Periods in Yaks

Biology ◽  
2021 ◽  
Vol 10 (12) ◽  
pp. 1229
Author(s):  
Shaoke Guo ◽  
Mengli Cao ◽  
Xingdong Wang ◽  
Lin Xiong ◽  
Xiaoyun Wu ◽  
...  

Yak reproductive characteristics have received extensive attention, though the molecular regulation mechanism of its ovarian activity remains to be explored. Therefore, this study initially conducted a comparative analysis of yak ovarian activities in anestrus, estrus, and pregnancy regarding their morphology and histology, followed by implementing RNA sequencing (RNA-seq) technology to detect the overall gene expression and biological mechanism in different reproductive stages. H&E staining showed that there were more growing follicles and mature follicles in ovarian tissue sections during estrus than ovarian tissues during non-estrus. The RNA-seq analysis of yak ovary tissues in three periods showed that DEGs related to follicular development and hormone metabolism were screened in the three comparison groups, such as COL1A2, NR4A1, THBS2, PTGS2, SCARB1, STAR, and WNT2B. Bioinformatics analysis showed that these DEGs are involved in ion binding, cell development, metabolic processes, enriched in ECM–receptor interactions, steroid biosynthesis, together with aldosterone generation/discharge and Wnt/PI3K-Akt signaling pathways. In addition, we speculate alternate splice development events to have important role/s in regulating ovarian functional genomic expression profiles. These results provide essential knowledge aimed at scrutinizing pivotal biomarkers for yak ovarian activity, together with paving the way for enhancing researchers’ focus on improving yak reproductive performance.

2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Samina Shabbir ◽  
Prerona Boruah ◽  
Lingli Xie ◽  
Muhammad Fakhar-e-Alam Kulyar ◽  
Mohsin Nawaz ◽  
...  

AbstractOvary development is an important determinant of the procreative capacity of female animals. Here, we performed genome-wide sequencing of long non-coding RNAs (lncRNAs) and mRNAs on ovaries of 1, 3 and 8 months old Hu sheep to assess their expression profiles and roles in ovarian development. We identified 37,309 lncRNAs, 45,404 messenger RNAs (mRNAs) and 330 novel micro RNAs (miRNAs) from the transcriptomic analysis. Six thousand, seven hundred and sixteen (6716) mRNAs and 1972 lncRNAs were significantly and differentially expressed in ovaries of 1 month and 3 months old Hu sheep (H1 vs H3). These mRNAs and target genes of lncRNAs were primarily enriched in the TGF-β and PI3K-Akt signalling pathways which are closely associated with ovarian follicular development and steroid hormone biosynthesis regulation. We identified MSTRG.162061.1, MSTRG.222844.7, MSTRG.335777.1, MSTRG.334059.16, MSTRG.188947.6 and MSTRG.24344.3 as vital genes in ovary development by regulating CTNNB1, CCNA2, CDK2, CDC20, CDK1 and EGFR expressions. A total of 2903 mRNAs and 636 lncRNAs were differentially expressed in 3 and 8 months old ovaries of Hu sheep (H3 vs H8); and were predominantly enriched in PI3K-Akt, progesterone-mediated oocyte maturation, estrogen metabolism, ovulation from the ovarian follicle and oogenesis pathways. These lncRNAs were also found to regulate FGF7, PRLR, PTK2, AMH and INHBA expressions during follicular development. Our result indicates the identified genes participate in the development of the final stages of follicles and ovary development in Hu sheep.


BMC Genomics ◽  
2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Yanlei Yue ◽  
Ze Jiang ◽  
Enoch Sapey ◽  
Tingting Wu ◽  
Shi Sun ◽  
...  

Abstract Background In soybean, some circadian clock genes have been identified as loci for maturity traits. However, the effects of these genes on soybean circadian rhythmicity and their impacts on maturity are unclear. Results We used two geographically, phenotypically and genetically distinct cultivars, conventional juvenile Zhonghuang 24 (with functional J/GmELF3a, a homolog of the circadian clock indispensable component EARLY FLOWERING 3) and long juvenile Huaxia 3 (with dysfunctional j/Gmelf3a) to dissect the soybean circadian clock with time-series transcriptomal RNA-Seq analysis of unifoliate leaves on a day scale. The results showed that several known circadian clock components, including RVE1, GI, LUX and TOC1, phase differently in soybean than in Arabidopsis, demonstrating that the soybean circadian clock is obviously different from the canonical model in Arabidopsis. In contrast to the observation that ELF3 dysfunction results in clock arrhythmia in Arabidopsis, the circadian clock is conserved in soybean regardless of the functional status of J/GmELF3a. Soybean exhibits a circadian rhythmicity in both gene expression and alternative splicing. Genes can be grouped into six clusters, C1-C6, with different expression profiles. Many more genes are grouped into the night clusters (C4-C6) than in the day cluster (C2), showing that night is essential for gene expression and regulation. Moreover, soybean chromosomes are activated with a circadian rhythmicity, indicating that high-order chromosome structure might impact circadian rhythmicity. Interestingly, night time points were clustered in one group, while day time points were separated into two groups, morning and afternoon, demonstrating that morning and afternoon are representative of different environments for soybean growth and development. However, no genes were consistently differentially expressed over different time-points, indicating that it is necessary to perform a circadian rhythmicity analysis to more thoroughly dissect the function of a gene. Moreover, the analysis of the circadian rhythmicity of the GmFT family showed that GmELF3a might phase- and amplitude-modulate the GmFT family to regulate the juvenility and maturity traits of soybean. Conclusions These results and the resultant RNA-seq data should be helpful in understanding the soybean circadian clock and elucidating the connection between the circadian clock and soybean maturity.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Xiaoqian Zhang ◽  
Chang Li ◽  
Bingzhou Zhang ◽  
Zhonghua Li ◽  
Wei Zeng ◽  
...  

AbstractThe variant virulent porcine epidemic diarrhea virus (PEDV) strain (YN15) can cause severe porcine epidemic diarrhea (PED); however, the attenuated vaccine-like PEDV strain (YN144) can induce immunity in piglets. To investigate the differences in pathogenesis and epigenetic mechanisms between the two strains, differential expression and correlation analyses of the microRNA (miRNA) and mRNA in swine testicular (ST) cells infected with YN15, YN144, and mock were performed on three comparison groups (YN15 vs Control, YN144 vs Control, and YN15 vs YN144). The mRNA and miRNA expression profiles were obtained using next-generation sequencing (NGS), and the differentially expressed (DE) (p-value < 0.05) mRNA and miRNA were obtained using DESeq R package. mRNAs targeted by DE miRNAs were predicted using the miRanda algortithm. 8039, 8631 and 3310 DE mRNAs, and 36, 36, and 22 DE miRNAs were identified in the three comparison groups, respectively. 14,140, 15,367 and 3771 DE miRNA–mRNA (targeted by DE miRNAs) interaction pairs with negatively correlated expression patterns were identified, and interaction networks were constructed using Cytoscape. Six DE miRNAs and six DE mRNAs were randomly selected to verify the sequencing data by real-time relative quantitative reverse transcription polymerase chain reaction (qRT-PCR). Based on bioinformatics analysis, we discovered the differences were mostly involved in host immune responses and viral pathogenicity, including NF-κB signaling pathway and bacterial invasion of epithelial cells, etc. This is the first comprehensive comparison of DE miRNA–mRNA pairs in YN15 and YN144 infection in vitro, which could provide novel strategies for the prevention and control of PED.


Cells ◽  
2021 ◽  
Vol 10 (2) ◽  
pp. 423
Author(s):  
Genxi Zhang ◽  
Mingliang He ◽  
Pengfei Wu ◽  
Xinchao Zhang ◽  
Kaizhi Zhou ◽  
...  

microRNAs play an important role in the growth and development of chicken embryos, including the regulation of skeletal muscle genesis, myoblast proliferation, differentiation, and apoptosis. Our previous RNA-seq studies showed that microRNA-27b-3p (miR-27b-3p) might play an important role in regulating the proliferation and differentiation of chicken primary myoblasts (CPMs). However, the mechanism of miR-27b-3p regulating the proliferation and differentiation of CPMs is still unclear. In this study, the results showed that miR-27b-3p significantly promoted the proliferation of CPMs and inhibited the differentiation of CPMs. Then, myostatin (MSTN) was confirmed to be the target gene of miR-27b-3p by double luciferase reporter assay, RT-qPCR, and Western blot. By overexpressing and interfering with MSTN expression in CPMs, the results showed that overexpression of MSTN significantly inhibited the proliferation and differentiation of CPMs. In contrast, interference of MSTN expression had the opposite effect. This study showed that miR-27b-3p could promote the proliferation of CPMs by targeting MSTN. Interestingly, both miR-27b-3p and MSTN can inhibit the differentiation of CPMs. These results provide a theoretical basis for further understanding the function of miR-27b-3p in chicken and revealing its regulation mechanism on chicken muscle growth.


BMC Genomics ◽  
2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Qi Wu ◽  
Yiming Luo ◽  
Xiaoyong Wu ◽  
Xue Bai ◽  
Xueling Ye ◽  
...  

Abstract Background Night-break (NB) has been proven to repress flowering of short-day plants (SDPs). Long-noncoding RNAs (lncRNAs) play key roles in plant flowering. However, investigation of the relationship between lncRNAs and NB responses is still limited, especially in Chenopodium quinoa, an important short-day coarse cereal. Results In this study, we performed strand-specific RNA-seq of leaf samples collected from quinoa seedlings treated by SD and NB. A total of 4914 high-confidence lncRNAs were identified, out of which 91 lncRNAs showed specific responses to SD and NB. Based on the expression profiles, we identified 17 positive- and 7 negative-flowering lncRNAs. Co-expression network analysis indicated that 1653 mRNAs were the common targets of both types of flowering lncRNAs. By mapping these targets to the known flowering pathways in model plants, we found some pivotal flowering homologs, including 2 florigen encoding genes (FT (FLOWERING LOCUS T) and TSF (TWIN SISTER of FT) homologs), 3 circadian clock related genes (EARLY FLOWERING 3 (ELF3), LATE ELONGATED HYPOCOTYL (LHY) and ELONGATED HYPOCOTYL 5 (HY5) homologs), 2 photoreceptor genes (PHYTOCHROME A (PHYA) and CRYPTOCHROME1 (CRY1) homologs), 1 B-BOX type CONSTANS (CO) homolog and 1 RELATED TO ABI3/VP1 (RAV1) homolog, were specifically affected by NB and competed by the positive and negative-flowering lncRNAs. We speculated that these potential flowering lncRNAs may mediate quinoa NB responses by modifying the expression of the floral homologous genes. Conclusions Together, the findings in this study will deepen our understanding of the roles of lncRNAs in NB responses, and provide valuable information for functional characterization in future.


Animals ◽  
2021 ◽  
Vol 11 (6) ◽  
pp. 1745
Author(s):  
Ben-Ben Miao ◽  
Su-Fang Niu ◽  
Ren-Xie Wu ◽  
Zhen-Bang Liang ◽  
Bao-Gui Tang ◽  
...  

Pearl gentian grouper (Epinephelus fuscoguttatus ♀ × Epinephelus lanceolatus ♂) is a fish of high commercial value in the aquaculture industry in Asia. However, this hybrid fish is not cold-tolerant, and its molecular regulation mechanism underlying cold stress remains largely elusive. This study thus investigated the liver transcriptomic responses of pearl gentian grouper by comparing the gene expression of cold stress groups (20, 15, 12, and 12 °C for 6 h) with that of control group (25 °C) using PacBio SMRT-Seq and Illumina RNA-Seq technologies. In SMRT-Seq analysis, a total of 11,033 full-length transcripts were generated and used as reference sequences for further RNA-Seq analysis. In RNA-Seq analysis, 3271 differentially expressed genes (DEGs), two low-temperature specific modules (tan and blue modules), and two significantly expressed gene sets (profiles 0 and 19) were screened by differential expression analysis, weighted gene co-expression networks analysis (WGCNA), and short time-series expression miner (STEM), respectively. The intersection of the above analyses further revealed some key genes, such as PCK, ALDOB, FBP, G6pC, CPT1A, PPARα, SOCS3, PPP1CC, CYP2J, HMGCR, CDKN1B, and GADD45Bc. These genes were significantly enriched in carbohydrate metabolism, lipid metabolism, signal transduction, and endocrine system pathways. All these pathways were linked to biological functions relevant to cold adaptation, such as energy metabolism, stress-induced cell membrane changes, and transduction of stress signals. Taken together, our study explores an overall and complex regulation network of the functional genes in the liver of pearl gentian grouper, which could benefit the species in preventing damage caused by cold stress.


BMC Genomics ◽  
2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Ruining Liu ◽  
Gang Li ◽  
Haoli Ma ◽  
Xianlong Zhou ◽  
Pengcheng Wang ◽  
...  

Abstract Background Ventilator-induced diaphragmatic dysfunction (VIDD) is associated with weaning difficulties, intensive care unit hospitalization (ICU), infant mortality, and poor long-term clinical outcomes. The expression patterns of long noncoding RNAs (lncRNAs) and mRNAs in the diaphragm in a rat controlled mechanical ventilation (CMV) model, however, remain to be investigated. Results The diaphragms of five male Wistar rats in a CMV group and five control Wistar rats were used to explore lncRNA and mRNA expression profiles by RNA-sequencing (RNA-seq). Muscle force measurements and immunofluorescence (IF) staining were used to verify the successful establishment of the CMV model. A total of 906 differentially expressed (DE) lncRNAs and 2,139 DE mRNAs were found in the CMV group. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to determine the biological functions or pathways of these DE mRNAs. Our results revealed that these DE mRNAs were related mainly related to complement and coagulation cascades, the PPAR signaling pathway, cholesterol metabolism, cytokine-cytokine receptor interaction, and the AMPK signaling pathway. Some DE lncRNAs and DE mRNAs determined by RNA-seq were validated by quantitative real-time polymerase chain reaction (qRT-PCR), which exhibited trends similar to those observed by RNA-sEq. Co-expression network analysis indicated that three selected muscle atrophy-related mRNAs (Myog, Trim63, and Fbxo32) were coexpressed with relatively newly discovered DE lncRNAs. Conclusions This study provides a novel perspective on the molecular mechanism of DE lncRNAs and mRNAs in a CMV model, and indicates that the inflammatory signaling pathway and lipid metabolism may play important roles in the pathophysiological mechanism and progression of VIDD.


Genes ◽  
2021 ◽  
Vol 12 (6) ◽  
pp. 896
Author(s):  
Yuenan Zhou ◽  
Pei Yang ◽  
Shuang Xie ◽  
Min Shi ◽  
Jianhua Huang ◽  
...  

The endoparasitic wasp Cotesia vestalis is an important biological agent for controlling the population of Plutella xylostella, a major pest of cruciferous crops worldwide. Though the genome of C. vestalis has recently been reported, molecular mechanisms associated with sexual development have not been comprehensively studied. Here, we combined PacBio Iso-Seq and Illumina RNA-Seq to perform genome-wide profiling of pharate adult and adult development of male and female C. vestalis. Taking advantage of Iso-Seq full-length reads, we identified 14,466 novel transcripts as well as 8770 lncRNAs, with many lncRNAs showing a sex- and stage-specific expression pattern. The differentially expressed gene (DEG) analyses showed 2125 stage-specific and 326 sex-specific expressed genes. We also found that 4819 genes showed 11,856 alternative splicing events through combining the Iso-Seq and RNA-Seq data. The results of comparative analyses showed that most genes were alternatively spliced across developmental stages, and alternative splicing (AS) events were more prevalent in females than in males. Furthermore, we identified six sex-determining genes in this parasitic wasp and verified their sex-specific alternative splicing profiles. Specifically, the characterization of feminizer and doublesex splicing between male and female implies a conserved regulation mechanism of sexual development in parasitic wasps.


2021 ◽  
Vol 11 (1) ◽  
Author(s):  
Zefang Sun ◽  
Jia Tan ◽  
Minqiong Zhao ◽  
Qiyao Peng ◽  
Mingqing Zhou ◽  
...  

AbstracttRNAs and tRNA-derived RNA fragments (tRFs) play various roles in many cellular processes outside of protein synthesis. However, comprehensive investigations of tRNA/tRF regulation are rare. In this study, we used new algorithms to extensively analyze the publicly available data from 1332 ChIP-Seq and 42 small-RNA-Seq experiments in human cell lines and tissues to investigate the transcriptional and posttranscriptional regulatory mechanisms of tRNAs. We found that histone acetylation, cAMP, and pluripotency pathways play important roles in the regulation of the tRNA gene transcription in a cell-specific manner. Analysis of RNA-Seq data identified 950 high-confidence tRFs, and the results suggested that tRNA pools are dramatically distinct across the samples in terms of expression profiles and tRF composition. The mismatch analysis identified new potential modification sites and specific modification patterns in tRNA families. The results also show that RNA library preparation technologies have a considerable impact on tRNA profiling and need to be optimized in the future.


2021 ◽  
Vol 22 (1) ◽  
Author(s):  
Chunyan Li ◽  
Xiaoyun He ◽  
Zijun Zhang ◽  
Chunhuan Ren ◽  
Mingxing Chu

Abstract Background Long noncoding RNA (lncRNA) has been identified as important regulator in hypothalamic-pituitary-ovarian axis associated with sheep prolificacy. However, little is known of their expression pattern and potential roles in the pineal gland of sheep. Herein, RNA-Seq was used to detect transcriptome expression pattern in pineal gland between follicular phase (FP) and luteal phase (LP) in FecBBB (MM) and FecB++ (ww) STH sheep, respectively, and differentially expressed (DE) lncRNAs and mRNAs associated with reproduction were identified. Results Overall, 135 DE lncRNAs and 1360 DE mRNAs in pineal gland between MM and ww sheep were screened. Wherein, 39 DE lncRNAs and 764 DE mRNAs were identified (FP vs LP) in MM sheep, 96 DE lncRNAs and 596 DE mRNAs were identified (FP vs LP) in ww sheep. Moreover, GO and KEGG enrichment analysis indicated that the targets of DE lncRNAs and DE mRNAs were annotated to multiple biological processes such as phototransduction, circadian rhythm, melanogenesis, GSH metabolism and steroid biosynthesis, which directly or indirectly participate in hormone activities to affect sheep reproductive performance. Additionally, co-expression of lncRNAs-mRNAs and the network construction were performed based on correlation analysis, DE lncRNAs can modulate target genes involved in related pathways to affect sheep fecundity. Specifically, XLOC_466330, XLOC_532771, XLOC_028449 targeting RRM2B and GSTK1, XLOC_391199 targeting STMN1, XLOC_503926 targeting RAG2, XLOC_187711 targeting DLG4 were included. Conclusion All of these differential lncRNAs and mRNAs expression profiles in pineal gland provide a novel resource for elucidating regulatory mechanism underlying STH sheep prolificacy.


Sign in / Sign up

Export Citation Format

Share Document