Identification of unique transcriptomic signatures through integrated multispecies comparative analysis and WGCNA in bovine oocyte development.
Ontology highlight
ABSTRACT: Background
Cattle (Bos taurus) are a major large livestock, however, compared with other species, the transcriptional specificity of bovine oocyte development has not been emphasised.Results
To reveal the unique transcriptional signatures of bovine oocyte development, we used integrated multispecies comparative analysis and weighted gene co-expression network analysis (WGCNA) to perform bioinformatic analysis of the germinal follicle (GV) and second meiosis (MII) gene expression profile from cattle, sheep, pigs and mice. We found that the expression levels of most genes were down-regulated from GV to MII in all species. Next, the multispecies comparative analysis showed more genes involved in the regulation of cAMP signalling during bovine oocyte development. Moreover, the green module identified by WGCNA was closely related to bovine oocyte development. Finally, integrated multispecies comparative analysis and WGCNA picked up 61 bovine-specific signature genes that participate in metabolic regulation and steroid hormone biosynthesis.Conclusion
In a short, this study provides new insights into the regulation of cattle oocyte development from a cross-species comparison.
SUBMITTER: Zhang FL
PROVIDER: S-EPMC10193610 | biostudies-literature | 2023 May
REPOSITORIES: biostudies-literature
Publications
Identification of unique transcriptomic signatures through integrated multispecies comparative analysis and WGCNA in bovine oocyte development.
Zhang Fa-Li FL Li Wei-Dong WD Zhang Geng G Zhang Min M Liu Zhao-Jun ZJ Zhu Ke-Xin KX Liu Qing-Chun QC Zhang Shu-Er SE Shen Wei W Zhang Xi-Feng XF
BMC genomics 20230518 1
<h4>Background</h4>Cattle (Bos taurus) are a major large livestock, however, compared with other species, the transcriptional specificity of bovine oocyte development has not been emphasised.<h4>Results</h4>To reveal the unique transcriptional signatures of bovine oocyte development, we used integrated multispecies comparative analysis and weighted gene co-expression network analysis (WGCNA) to perform bioinformatic analysis of the germinal follicle (GV) and second meiosis (MII) gene expression ...[more]
Similar Datasets
Project description:PurposeTranscriptomic studies on gastroesophageal reflux disease are scarce, and gene expression signatures in nonerosive reflux disease (NERD) remain elusive. The aim of the study was to identify gene expression profiles and potential hub genes in NERD.Patients and methodsWe performed RNA sequencing on biopsy samples from nine consecutive patients with NERD and six healthy controls. Differentially expressed genes (DEGs) were analysed with the DESeq2 R package. A DEG-based protein-protein interaction (PPI) network was constructed to filter hub genes using Cytoscape. Weighted gene coexpression network analysis (WGCNA) was conducted to identify the coexpression relationships of all modules and explore the relationship between gene sets and clinical traits.ResultsIn total, 1195 DEGs were identified, including 649 upregulated and 546 downregulated genes involved in regulating the inflammatory response and epithelial cell differentiation. Overlap of the PPI and WGCNA networks identified five shared genes, namely, THY1, BMP2, LOX, KDR and MMP9, as candidate hub genes in NERD. Quantitative PCR analysis of the expression of these five genes confirmed the sequencing results. Receiver operating characteristic analyses indicated that these hub genes had diagnostic potential for NERD patients. Gene set enrichment analysis confirmed that each hub gene was closely associated with the pathophysiological processes of NERD. In addition, a regulatory network comprising 42 transcription factors (TFs), 28 miRNAs and 5 hub genes was established.ConclusionThe five core genes may be promising biomarkers of NERD. The TF/miRNA/hub gene network can improve the understanding of the molecular mechanisms underlying disease progression.
| S-EPMC8627320 | biostudies-literature
Project description:BackgroundBovine papillomavirus (BPV) belongs to the Papillomaviridae family and infects epithelial cells of bovines and closely related animals, causing hyperproliferative lesions known as warts or papillomas, which may regress or progress to form benign or malignant tumors. The virus enters the host cell and interacts with it by altering the regulation of genes that are responsible for controlling the cell cycle, thus triggering lesion formation. It is not yet known which host genes are regulated by viral infection. Therefore, the objective of this study was to make use of next-generation RNA sequencing methods to identify differentially expressed genes associated with BPV infection, which might elucidate possible marker genes that could be used to control the disease.ResultsTranscriptome analysis revealed that 1343 genes were differentially regulated (FDR < 0.05). A comparison of gene expression in infected and noninfected cows indicated that 655 genes were significantly upregulated, and 688 genes were significantly downregulated. Most differentially expressed genes were associated with BPV infection pathways, which supports the hypothesis that viral infection was the mechanism associated with this regulation.ConclusionsThis is the first study that focused on a large-scale evaluation of gene expression associated with BPV infection, which is important to identify possible metabolic pathways regulated by host genes for lesion development. In addition, novel targets could be identified in order to find ligands that interact with BPV, with the aim of interrupting the infection cycle.
| S-EPMC6300001 | biostudies-literature
Project description:BackgroundA risk assessment model for prognostic prediction of colon adenocarcinoma (COAD) was established based on weighted gene co-expression network analysis (WGCNA).MethodsFrom the Cancer Genome Atlas (TCGA) database, RNA-seq data and clinical data of COAD patients were retrieved. After screening of differentially expressed genes (DEGs), WGCNA was performed to identify gene modules and screen those associated with COAD progression. Then, via protein-protein interaction (PPI) network construction of module genes, hub genes were obtained, which were then subjected to the least absolute shrinkage and selection operator (LASSO) and Cox regression to build a hub gene-based prognostic scoring model. The receiver operating characteristic curve (ROC curve) was plotted for the optimal cutoff (OCO) of the risk score, based on which, patients were assigned to high or low-risk groups. Areas under the ROC curve (AUCs) were calculated, and model performance was visualized using Kaplan-Meier (KM) survival curves and verified in the external dataset GSE29621. Finally, the model's independent prognostic value was evaluated by univariate and multivariate Cox regression analyses, and a nomogram was built.ResultsTotally 2840 DEGs were screened from COAD dataset of TCGA, including 1401 upregulated ones and 1439 downregulated ones, which were divided into 10 modules by WGCNA. The eigenvalue of the black module was found to have a high correlation with COAD progression. PPI interaction networks were constructed for genes in the black module, and 34 hub genes were obtained by using the MCODE plug-in. A LASSO-Cox regression approach was utilized to analyze the hub genes, and a prognostic risk score model based on the signatures of 9 genes (CHEK1, DEPDC1B, FANCI, MCM10, NCAPG, PARPBP, PLK4, RAD51AP1, and RFC4) was constructed. KM analysis identified shorter overall lower survival in the high-risk group. The model was verified to have favorable predictive ability through training set and validation set. The nomogram, composed of tumor node metastasis (TNM) staging and risk score, was of good predictability.ConclusionsThe COAD prognostic risk model constructed upon the signatures of 9 genes (CHEK1, DEPDC1B, FANCI, MCM10, NCAPG, PARPBP, PLK4, RAD51AP1, and RFC4) can effectively predict the survival status of COAD patients.
| S-EPMC8983226 | biostudies-literature
Project description:Esophageal adenocarcinoma (EAC) is a deadly cancer with high mortality rate, especially in economically advanced countries, while Barrett's esophagus (BE) is reported to be a precursor that strongly increases the risk of EAC. Due to the complexity of these diseases, their molecular mechanisms have not been revealed clearly. This study aims to explore the gene signatures shared between BE and EAC based on integrated network analysis. We obtained EAC- and BE-associated microarray datasets GSE26886, GSE1420, GSE37200, and GSE37203 from the Gene Expression Omnibus and ArrayExpress using systematic meta-analysis. These data were accompanied by clinical data and RNAseq data from The Cancer Genome Atlas (TCGA). Weighted gene co-expression network analysis (WGCNA) and differentially expressed gene (DEG) analysis were conducted to explore the relationship between gene sets and clinical traits as well as to discover the key relationships behind the co-expression modules. A differentially expressed gene-based protein-protein interaction (PPI) complex was used to extract hub genes through Cytoscape plugins. As a result, 403 DEGs were excavated, comprising 236 upregulated and 167 downregulated genes, which are involved in the cell cycle and replication pathways. Forty key genes were identified using modules of MCODE, CytoHubba, and CytoNCA with different algorithms. A dark-gray module with 207 genes was identified which having a high correlation with phenotype (gender) in the WGCNA. Furthermore, five shared hub gene signatures (SHGS), namely, pre-mRNA processing factor 4 (PRPF4), serine and arginine-rich splicing factor 1 (SRSF1), heterogeneous nuclear ribonucleoprotein M (HNRNPM), DExH-Box Helicase 9 (DHX9), and origin recognition complex subunit 2 (ORC2), were identified between BE and EAC. SHGS enrichment denotes that RNA metabolism and splicosomes play a key role in esophageal cancer development and progress. We conclude that the PPI complex and WGCNA co-expression network highlight the importance of phenotypic identifying hub gene signatures for BE and EAC.
| S-EPMC7438937 | biostudies-literature
Project description:Litsea coreana Levl. var. sinensis (Allen) Yang et P. H. Huang is a popular ethnic herb and beverage plant known for its high flavonoid content, which has been linked to a variety of pharmacological benefits and crucial health-promoting impacts in humans. The progress in understanding the molecular mechanisms of flavonoid accumulation in this plant has been hindered due to the deficiency of genomic and transcriptomic resources. We utilized a combination of Illumina and Oxford Nanopore Technology (ONT) sequencing to generate a de novo hybrid transcriptome assembly. In total, 126,977 unigenes were characterized, out of which 107,977 were successfully annotated in seven public databases. Within the annotated unigenes, 3,781 were categorized into 58 transcription factor families. Furthermore, we investigated the presence of four valuable flavonoids-quercetin-3-O-β-D-galactoside, quercetin-3-O-β-D-glucoside, kaempferol-3-O-β-D-galactoside, and kaempferol-3-O-β-D-glucoside in 98 samples, using high-performance liquid chromatography. A weighted gene co-expression network analysis identified two co-expression modules, MEpink and MEturquoise, that showed strong positive correlation with flavonoid content. Within these modules, four transcription factor genes (R2R3-MYB, NAC, WD40, and ARF) and four key enzyme-encoding genes (CHI, F3H, PAL, and C4H) emerged as potential hub genes. Among them, the R2R3-MYB (LcsMYB123) as a homologous gene to AtMYB123/TT2, was speculated to play a significant role in flavonol biosynthesis based on phylogenetic analysis. Our findings provided a theoretical foundation for further research into the molecular mechanisms of flavonoid biosynthesis. Additionally, The hybrid transcriptome sequences will serve as a valuable molecular resource for the transcriptional annotation of L. coreana var. sinensis, which will contribute to the improvement of high-flavonoid materials.
| S-EPMC10985888 | biostudies-literature
Project description:Despite great improvements in assisted reproductive technology (ARTs), the success of in vitro embryo production remains relatively low. Most of the oocytes used to produce in vitro embryos are recovered from smaller follicles, forming a heterogeneous population that must be matured in vitro. Since the original in vitro maturation (IVM ) experiments (Heilbrunn et al. , 1939) the process by which the most competent oocytes are selected to produce blastocysts remains similar and is still based on morphological aspects of the oocyte cytoplasm and the number of layers and compaction of cumulus cells attached to the oocyte surface (Armstrong, 2001, Coticchio et al. , 2004, Krisher, 2003, Lonergan et al. , 2003). Ovaries from crossbred cows (Bos indicus x Bos taurus) were collected immediately after slaughter and transported to the laboratory in saline solution (0.9% NaCl) supplemented with penicillin G (100 IU/mL) and streptomycin sulfate (100 mg/mL) at 35-37C. The follicles were dissected from the ovarian cortex using scissors, scalpels and tweezers in TCM-199 medium supplemented with Hank’s salts and 10% fetal calf serum (FCS; GIBCO BRL) at 36C. Follicles were measured using a graduated eyepiece (OSM-4; Olympus) and then classified morphologically into two groups according to their diameter: 1.0-3.0 mm or ≥8.0 mm. The criteria used for follicle selection included the presence of extensive and fine vascularization and a shiny and translucent appearance. After follicular rupture, the presence of granulosa cells with a regular and healthy appearance and no free-floating particles in the follicular fluid were also used as selection criteria. Only COCs with a homogeneous granulated cytoplasm and at least three layers of compact of cumulus cells were used in the present study. The selected COCs were transferred to a 50 μL droplet of phosphate-buffered saline (PBS), and the cumulus cells were mechanically removed by repeated pipetting. After cumulus collection, they were transferred to a 0.2-mL tube and centrifuged twice for 2 min at 700 g. The supernatant was removed, and 2μL of RNAlater (Applied Biosystems, Foster City, CA, USA) was added to the pellet, which was then stored at -20C until RNA extraction. For each follicle size group, three replicas of pooled cumulus cells corresponding to 30 oocytes were stored for RNA extraction and subsequent microarray assays, and four independent replicates corresponding to 20 oocytes were stored for RNA extraction for subsequent qPCR assays.
2015-01-27 | GSE65269 | GEO
Project description:Primordial follicle activation (PFA) is a pivotal event in female reproductive biology, coordinating the transition from quiescent to growing follicles. This study employed comprehensive single-cell RNA sequencing to gain insights into the detailed regulatory mechanisms governing the synchronized dormancy and activation between granulosa cells (GCs) and oocytes with the progression of the PFA process. Wntless (Wls) conditional knockout (cKO) mice served as a unique model, suppressing the transition from pre-GCs to GCs, and disrupting somatic cell-derived WNT signaling in the ovary. Our data revealed immediate transcriptomic changes in GCs post-PFA in Wls cKO mice, leading to a divergent trajectory, while oocytes exhibited modest transcriptomic alterations. Subpopulation analysis identified the molecular pathways affected by WNT signaling on GC maturation, along with specific gene signatures linked to dormant and activated oocytes. Despite minimal evidence of continuous up-regulation of dormancy-related genes in oocytes, the loss of WNT signaling in (pre-)GCs impacted gene expression in oocytes even before PFA, subsequently influencing them globally. The infertility observed in Wls cKO mice was attributed to compromised GC-oocyte molecular crosstalk and the microenvironment for oocytes. Our study highlights the pivotal role of the WNT-signaling pathway and its molecular signature, emphasizing the importance of intercellular crosstalk between (pre-)GCs and oocytes in orchestrating folliculogenesis.
| S-EPMC11498688 | biostudies-literature
Project description:Despite great improvements in assisted reproductive technology (ARTs), the success of in vitro embryo production remains relatively low. Most of the oocytes used to produce in vitro embryos are recovered from smaller follicles, forming a heterogeneous population that must be matured in vitro. Since the original in vitro maturation (IVM ) experiments (Heilbrunn et al. , 1939) the process by which the most competent oocytes are selected to produce blastocysts remains similar and is still based on morphological aspects of the oocyte cytoplasm and the number of layers and compaction of cumulus cells attached to the oocyte surface (Armstrong, 2001, Coticchio et al. , 2004, Krisher, 2003, Lonergan et al. , 2003). Ovaries from crossbred cows (Bos indicus x Bos taurus) were collected immediately after slaughter and transported to the laboratory in saline solution (0.9% NaCl) supplemented with penicillin G (100 IU/mL) and streptomycin sulfate (100 mg/mL) at 35-37C. The follicles were dissected from the ovarian cortex using scissors, scalpels and tweezers in TCM-199 medium supplemented with Hank’s salts and 10% fetal calf serum (FCS; GIBCO BRL) at 36C. Follicles were measured using a graduated eyepiece (OSM-4; Olympus) and then classified morphologically into two groups according to their diameter: 1.0-3.0 mm or ≥8.0 mm. The criteria used for follicle selection included the presence of extensive and fine vascularization and a shiny and translucent appearance. After follicular rupture, the presence of granulosa cells with a regular and healthy appearance and no free-floating particles in the follicular fluid were also used as selection criteria. Only COCs with a homogeneous granulated cytoplasm and at least three layers of compact of cumulus cells were used in the present study. The selected COCs were transferred to a 50 μL droplet of phosphate-buffered saline (PBS), and the cumulus cells were mechanically removed by repeated pipetting. After cumulus collection, they were transferred to a 0.2-mL tube and centrifuged twice for 2 min at 700 g. The supernatant was removed, and 2μL of RNAlater (Applied Biosystems, Foster City, CA, USA) was added to the pellet, which was then stored at -20C until RNA extraction. For each follicle size group, three replicas of pooled cumulus cells corresponding to 30 oocytes were stored for RNA extraction and subsequent microarray assays, and four independent replicates corresponding to 20 oocytes were stored for RNA extraction for subsequent qPCR assays. During oogenesis, the somatic cells that surround the oocyte proliferate and differentiate into cumulus cells (CCs), which are metabolically coupled to the oocyte, forming the cumulus-oocyte complexes (COCs) (van den Hurk and Zhao, 2005). The CCs remain in strong contact with the oocyte, maintained by cytoplasmic bridges called gap-junctions, as well as through justacrine and paracrine signaling networks (Albertini et al. , 2001, Gilchrist et al. , 2004, Tanghe et al. , 2002, Vozzi et al. , 2001, Webb et al. , 2002). This intense bidirectional communication between the CCs and the oocyte is maintained during all phases of folliculogenesis and is essential for the acquisition of oocyte competence that is required for blastocyst development (Assidi et al., 2008, Fair, 2003, Gilchrist et al., 2004). An informative model to access the level of oocyte competence that has been used by many research groups is follicle size (Donnison and Pfeffer, 2004, Franco et al. , 2013, Lequarre et al. , 2005, Mourot et al., 2006). In a previous study, we showed that oocytes obtained from follicles 1-3 mm in diameter are significantly less competent in producing blastocysts than oocytes obtained from follicles ≥8 mm in diameter (Franco et al., 2013). In the present study, we use the same model to compare the gene expression profile of more than 23,000 bovine transcripts between the CCs obtained from incompetent COCs (1-3 mm) and competent COCs (≥8 mm). We have found substantial differences in the expression of several gene clusters representing distinct metabolic pathways such as energy metabolism, cell signaling, cell cycle, DNA repair, meiosis, and inflammation.
2015-01-27 | E-GEOD-65269 | biostudies-arrayexpress
Project description:Background & aimsWe used patient-derived organoids (PDOs) to study the epithelial-specific transcriptional and secretome signatures of the ileum during Crohn's disease (CD) with varying phenotypes to screen for disease profiles and potential druggable targets.MethodsRNA sequencing was performed on isolated intestinal crypts and 3-week-old PDOs derived from ileal biopsies of CD patients (n = 8 B1, inflammatory; n = 8 B2, stricturing disease) and non-inflammatory bowel disease (IBD) controls (n = 13). Differentially expressed (DE) genes were identified by comparing CD vs control, B1 vs B2, and inflamed vs non-inflamed. DE genes were used for computational screening to find candidate small molecules that could potentially reverse B1and B2 gene signatures. The secretome of a second cohort (n = 6 non-IBD controls, n = 7 CD, 5 non-inflamed, 2 inflamed) was tested by Luminex using cultured organoid conditioned medium.ResultsWe found 90% similarity in both the identity and abundance of protein coding genes between PDOs and intestinal crypts (15,554 transcripts of 19,900 genes). DE analysis identified 814 genes among disease group (CD vs non-IBD control), 470 genes different between the CD phenotypes, and 5 false discovery rate correction significant genes between inflamed and non-inflamed CD. The PDOs showed both similarity and diversity in the levels and types of soluble cytokines and growth factors they released. Perturbagen analysis revealed potential candidate compounds to reverse B2 disease phenotype to B1 in PDOs.ConclusionsPDOs are similar at the transcriptome level with the in vivo epithelium and retain disease-specific gene expression for which we have identified secretome products, druggable targets, and corresponding pharmacologic agents. Targeting the epithelium could reverse a stricturing phenotype and improve outcomes.
| S-EPMC8455365 | biostudies-literature
Project description:BackgroundEndothelial cells (ECs) display considerable functional heterogeneity depending on the vessel and tissue in which they are located. Whereas these functional differences are presumably imprinted in the transcriptome, the pathways and networks that sustain EC heterogeneity have not been fully delineated.MethodsTo investigate the transcriptomic basis of EC specificity, we analyzed single-cell RNA sequencing data from tissue-specific mouse ECs generated by the Tabula Muris consortium. We used a number of bioinformatics tools to uncover markers and sources of EC heterogeneity from single-cell RNA sequencing data.ResultsWe found a strong correlation between tissue-specific EC transcriptomic measurements generated by either single-cell RNA sequencing or bulk RNA sequencing, thus validating the approach. Using a graph-based clustering algorithm, we found that certain tissue-specific ECs cluster strongly by tissue (eg, liver, brain), whereas others (ie, adipose, heart) have considerable transcriptomic overlap with ECs from other tissues. We identified novel markers of tissue-specific ECs and signaling pathways that may be involved in maintaining their identity. Sex was a considerable source of heterogeneity in the endothelial transcriptome and we discovered Lars2 to be a gene that is highly enriched in ECs from male mice. We found that markers of heart and lung ECs in mice were conserved in human fetal heart and lung ECs. We identified potential angiocrine interactions between tissue-specific ECs and other cell types by analyzing ligand and receptor expression patterns.ConclusionsWe used single-cell RNA sequencing data generated by the Tabula Muris consortium to uncover transcriptional networks that maintain tissue-specific EC identity and to identify novel angiocrine and functional relationships between tissue-specific ECs.
| S-EPMC7658053 | biostudies-literature