scholarly journals Integrated analysis of lncRNA and mRNA reveals novel insights into cashmere fineness in Tibetan cashmere goats

PeerJ ◽  
2020 ◽  
Vol 8 ◽  
pp. e10217
Author(s):  
Xuefeng Fu ◽  
Bingru Zhao ◽  
Kechuan Tian ◽  
Yujiang Wu ◽  
Langda Suo ◽  
...  

Tibetan cashmere goats are famous for producing the finest, softest and lightest cashmere fiber in China. The growth and development of skin are closely related to fineness and are the key factors affecting the quality of cashmere. To investigate the specific role of long noncoding RNAs (lncRNAs) and messenger RNAs (mRNAs) in regulating cashmere fineness of Tibetan Cashmere goats in the anagen phase, we conducted high-throughput RNA sequencing of fine-type and coarse-type skin tissues. We identified 2,059 lncRNA candidates (1,589 lncRNAs annotated, 470 lncRNAs novel), and 80 differentially expressed (DE) lncRNAs and their potential targets were predicted. We also identified 384 DE messenger RNAs (mRNAs) out of 29,119 mRNAs. Several key genes in KRT26, KRT28, KRT39, IFT88, JAK3, NOTCH2 and NOTCH3 and a series of lncRNAs, including ENSCHIT00000009853, MSTRG.16794.17, MSTRG.17532.2, were shown to be potentially important for regulating cashmere fineness. GO and KEGG enrichment analyses of DE mRNAs and DE lncRNAs targets significantly enriched in positive regulation of the canonical Wnt signaling pathway, regulation of protein processing and metabolism processes. The mRNA-mRNA and lncRNA-mRNA regulatory networks further revealed potential transcripts involved in cashmere fineness. We further validated the expression patterns of DE mRNAs and DE lncRNAs by quantitative real-time PCR (qRT-PCR), and the results were consistent with the sequencing data. This study will shed new light on selective cashmere goat breeding, and these lncRNAs and mRNAs that were found to be enriched in Capra hircus RNA database.

2020 ◽  
Author(s):  
Silva Kasela ◽  
Victor E Ortega ◽  
Molly Martorella ◽  
Suresh Garudadri ◽  
Jenna Nguyen ◽  
...  

Introductory paragraphParticular host and environmental factors influence susceptibility to severe COVID-19. We analyzed RNA-sequencing data from bronchial epithelial brushings - a relevant tissue for SARS-CoV-2 infection - obtained from three cohorts of uninfected individuals, and investigated how non-genetic and genetic factors affect the regulation of host genes implicated in COVID-19. We found that ACE2 expression was higher in relation to active smoking, obesity, and hypertension that are known risk factors of COVID-19 severity, while an association with interferon-related inflammation was driven by the truncated, non-binding ACE2 isoform. We discovered that expression patterns of a suppressed airway immune response to early SARS-CoV-2 infection, compared to other viruses, are similar to patterns associated with obesity, hypertension, and cardiovascular disease, which may thus contribute to a COVID-19-susceptible airway environment. eQTL mapping identified regulatory variants for genes implicated in COVID-19, some of which had pheWAS evidence for their potential role in respiratory infections. These data provide evidence that clinically relevant variation in the expression of COVID-19-related genes is associated with host factors, environmental exposures, and likely host genetic variation.


2020 ◽  
Author(s):  
Liliang Yang ◽  
Kaizhen Wang ◽  
Wenjing Guo ◽  
Xian Chen ◽  
Qinglong Guo ◽  
...  

Abstract Background:RNA polymerase II subunit K (POLR2K) belongs to one of the multiple subunits of RNA polymerase II (Pol II), whose biological function is to synthesize mRNA. Aberrant POLR2K expression is related to carcinogenesis. However, POLR2K’s underlying role in bladder cancer has not been explored. In the current study, we intend to analyze the function of POLR2K and its regulatory network within bladder cancer.Methods: Public sequencing data was obtained from GEO and TCGA to investigate POLR2K expression and regulatory network within bladder cancer (BLCA) by using GEPIA and Oncomine as well as cBioPortal online tool. LinkedOmics was employed to identify genes displaying significantly differential expression patterns and to perform GO and KEGG analyses. After differential genes was assigned and ranked, GSEA analyses was performed to obtain target networks for transcription factors, miRNAs, and kinases that could regulate POLR2K–associated gene network. Subsequent functional webwork analyses were used to identify cancer-relevant pathways Moreover, POLR2K gene is verified, by ChIP-seq in MCF-7 cell line , with transcription factor binding evidence in the ENCODE Transcription Factor Binding Site Profiles dataset.Conclusions: The current study implies that POLR2K gene is overexpressed and often amplified in BLCA, providing the first evidence that POLR2K deregulation, in particular increased transcription, may promote BLCA. These findings uncover a unique expression patterns of POLR2K and its potential regulatory networks in BLCA, contributing greatly to study of the role of POLR2K in cancer development.


Genes ◽  
2019 ◽  
Vol 10 (4) ◽  
pp. 266 ◽  
Author(s):  
Yuan Y. Zheng ◽  
Sheng D. Sheng ◽  
Tai Y. Hui ◽  
Chang Yue ◽  
Jia M. Sun ◽  
...  

Animal growth and development are regulated by long non-coding RNAs (lncRNAs). However, the functions of lncRNAs in regulating cashmere fineness are poorly understood. To identify the key lncRNAs that are related to cashmere fineness in skin, we have collected skin samples of Liaoning cashmere goats (LCG) and Inner Mongolia cashmere goats (MCG) in the anagen phase, and have performed RNA sequencing (RNA-seq) approach on these samples. The high-throughput sequencing and bioinformatics analyses identified 437 novel lncRNAs, including 93 differentially expressed lncRNAs. We also identified 3,084 differentially expressed messenger RNAs (mRNAs) out of 27,947 mRNAs. Gene ontology (GO) analyses of lncRNAs and target genes in cis show a predominant enrichment of targets that are related to intermediate filament and intermediate filament cytoskeleton. According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, sphingolipid metabolism is a significant pathway for lncRNA targets. In addition, this is the first report to reveal the possible lncRNA–mRNA regulatory network for cashmere fineness in cashmere goats. We also found that lncRNA XLOC_008679 and its target gene, KRT35, may be related to cashmere fineness in the anagen phase. The characterization and expression analyses of lncRNAs will facilitate future studies on the potential value of fiber development in LCG.


2020 ◽  
Vol 2020 ◽  
pp. 1-27
Author(s):  
Xiaoqing Wang ◽  
Xiangyu Wu ◽  
Tianyou Li ◽  
Mingyu Cui ◽  
Lichao Zhu ◽  
...  

Rhabdoid tumor of the kidney (RTK) is a rare and severely malignant tumor occurring in infancy and early childhood, with the overall outcomes remain poor. Neither gene regulatory networks nor biomarkers to predict the prognostic outcomes have been elucidated in RTK. In this study, RNA sequencing data were obtained to identify differentially expressed messenger RNAs (mRNAs), long noncoding RNAs (lncRNAs), and microRNAs (miRNAs) between RTK samples and normal samples. A total of 4217 mRNAs, 284 lncRNAs, and 286 miRNAs were screened out. Of those, 103 mRNAs, 80 lncRNAs, and 45 miRNAs were identified for a competing endogenous RNA (ceRNA) regulatory network, in which three significant modules were identified. A protein-protein interaction (PPI) network was constructed, and the hub-gene cluster consisted of four core genes (EXOSC2, PAK1IP1, WDR43, and POLR1D) was selected. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were also performed to analyze the functional characteristics of differentially expressed mRNAs. Subsequently, among 211 mRNAs, 8 lncRNAs, and 12 miRNAs associated with overall survival (OS) obtained by univariate Cox analysis, 5 mRNAs, 7 lncRNAs, and 7 miRNAs were identified and the risk score formulas were constructed correspondingly using the least absolute shrinkage and selection operator (LASSO) Cox regression model analysis. The log-rank tests and Kaplan-Meier analyses were performed to confirm the predictive value of the risk scores for OS in RTK patients. A genomic-clinicopathologic nomogram integrating the stage and risk scores based on RNAs was established and demonstrated high predictive accuracy and clinical value, which was validated through calibration curves, time-dependent receiver operating characteristic (ROC) curve analyses, and decision curve analysis (DCA). In conclusion, this study not only provided potential insights into the mechanisms underlying RTK, but also presented a practicable tool for predicting the prognosis in children with RTK.


2019 ◽  
Vol 19 (1) ◽  
Author(s):  
Li Xue ◽  
Jian Wang ◽  
Jun Zhao ◽  
Yang Zheng ◽  
Hai-Feng Wang ◽  
...  

Abstract Background Pink-flowered strawberry is a promising new ornamental flower derived from intergeneric hybridization (Fragaria × Potentilla) with bright color, a prolonged flowering period and edible fruits. Its flower color ranges from light pink to red. Pigment compounds accumulated in its fruits were the same as in cultivated strawberry fruits, but different from that in its flowers. However, the transcriptional events underlying the anthocyanin biosynthetic pathway have not been fully characterized in petal coloration. To gain insights into the regulatory networks related to anthocyanin biosynthesis and identify the key genes, we performed an integrated analysis of the transcriptome and metabolome in petals of pink-flowered strawberry. Results The main pigments of red and dark pink petals were anthocyanins, among which cyanidins were the main compound. There were no anthocyanins detected in the white-flowered hybrids. A total of 50,285 non-redundant unigenes were obtained from the transcriptome databases involved in red petals of pink-flowered strawberry cultivar Sijihong at three development stages. Amongst the unigenes found to show significant differential expression, 57 were associated with anthocyanin or other flavonoid biosynthesis, in which they were regulated by 241 differentially expressed members of transcription factor families, such as 40 MYBs, 47 bHLHs, and 41 NACs. Based on a comprehensive analysis relating pigment compounds to gene expression profiles, the mechanism of flower coloration was examined in pink-flowered strawberry. A new hypothesis was proposed to explain the lack of color phenotype of the white-flowered strawberry hybrids based on the transcriptome analysis. The expression patterns of FpDFR and FpANS genes corresponded to the accumulation patterns of cyanidin contents in pink-flowered strawberry hybrids with different shades of pink. Moreover, FpANS, FpBZ1 and FpUGT75C1 genes were the major factors that led to the absence of anthocyanins in the white petals of pink-flowered strawberry hybrids. Meanwhile, the competitive effect of FpFLS and FpDFR genes might further inhibit anthocyanin synthesis. Conclusions The data presented herein are important for understanding the molecular mechanisms underlying the petal pigmentation and will be powerful for integrating novel potential target genes to breed valuable pink-flowered strawberry cultivars.


2021 ◽  
Vol 5 (1) ◽  
Author(s):  
Robert L. Hollis ◽  
Barbara Stanley ◽  
John P. Thomson ◽  
Michael Churchman ◽  
Ian Croy ◽  
...  

AbstractEndometrioid ovarian carcinoma (EnOC) is an under-investigated ovarian cancer type. Recent studies have described disease subtypes defined by genomics and hormone receptor expression patterns; here, we determine the relationship between these subtyping layers to define the molecular landscape of EnOC with high granularity and identify therapeutic vulnerabilities in high-risk cases. Whole exome sequencing data were integrated with progesterone and oestrogen receptor (PR and ER) expression-defined subtypes in 90 EnOC cases following robust pathological assessment, revealing dominant clinical and molecular features in the resulting integrated subtypes. We demonstrate significant correlation between subtyping approaches: PR-high (PR + /ER + , PR + /ER−) cases were predominantly CTNNB1-mutant (73.2% vs 18.4%, P < 0.001), while PR-low (PR−/ER + , PR−/ER−) cases displayed higher TP53 mutation frequency (38.8% vs 7.3%, P = 0.001), greater genomic complexity (P = 0.007) and more frequent copy number alterations (P = 0.001). PR-high EnOC patients experience favourable disease-specific survival independent of clinicopathological and genomic features (HR = 0.16, 95% CI 0.04–0.71). TP53 mutation further delineates the outcome of patients with PR-low tumours (HR = 2.56, 95% CI 1.14–5.75). A simple, routinely applicable, classification algorithm utilising immunohistochemistry for PR and p53 recapitulated these subtypes and their survival profiles. The genomic profile of high-risk EnOC subtypes suggests that inhibitors of the MAPK and PI3K-AKT pathways, alongside PARP inhibitors, represent promising candidate agents for improving patient survival. Patients with PR-low TP53-mutant EnOC have the greatest unmet clinical need, while PR-high tumours—which are typically CTNNB1-mutant and TP53 wild-type—experience excellent survival and may represent candidates for trials investigating de-escalation of adjuvant chemotherapy to agents such as endocrine therapy.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Yu Qiao ◽  
Qiming Cheng ◽  
Yutong Zhang ◽  
Wei Yan ◽  
Fengyan Yi ◽  
...  

Abstract Background Sainfoin (Onobrychis viciifolia Scop) is not only a high-quality legume forage, but also a nectar-producing plant. Therefore, the flower color of sainfoin is an important agronomic trait, but the factors affecting its flower phenotype are still unclear. To gain insights into the regulatory networks associated with metabolic pathways of coloration compounds (flavonoids or anthocyanins) and identify the key genes, we conducted a comprehensive analysis of the phenotype, metabolome and transcriptome of WF and AF of sainfoin. Results Delphinidin, petunidin and malvidin derivatives were the main anthocyanin compounds in the AF of sainfoin. These substances were not detected in the WF of sainfoin. The transcriptomes of WF and AF in sainfoin at the S1 and S3 stages were obtained using the Illumina HiSeq4000 platform. Overall, 10,166 (4273 upregulated and 5893 downregulated) and 15,334 (8174 upregulated and 7160 downregulated) DEGs were identified in flowers at S1 and S3 stages, respectively (WF-VS-AF). KEGG pathway annotations showed that 6396 unigenes were annotated to 120 pathways and contained 866 DEGs at S1 stages, and 6396 unigenes were annotated to 131 pathways and included 1546 DEGs at the S3 stage. Nine DEGs belonging to the “flavonoid biosynthesis”and “phenylpropanoid biosynthesis” pathways involved in flower color formation were identified and verified by RT-qPCR analyses. Among these DEGs, 4CL3, FLS, ANS, CHS, DFR and CHI2 exhibited downregulated expression, and F3H exhibited upregulated expression in the WF compared to the AF, resulting in a decrease in anthocyanin synthesis and the formation of WF in sainfoin. Conclusions This study is the first to use transcriptome technology to study the mechanism of white flower formation in sainfoin. Our transcriptome data will be a great enrichment of the genetic information for sainfoin. In addition, the data presented herein will provide valuable molecular information for genetic breeding and provide insight into the future study of flower color polymorphisms in sainfoin.


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.


Sign in / Sign up

Export Citation Format

Share Document